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ECOLOGICAL  STUDIES  OF  LARVAL  Glyptotendipes  paripes  (CHIRONOMIDAE: 

DIPTERA)  IN  SELECTED  CENTRAL  FLORIDA  LAKES  FOR  CREATING  AN 

EXPLORATORY  TEMPORAL  AND  SPATIAL  MODEL  OF  NUISANCE 
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By 
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Chairman:  Arshad  Ali 

Major  Department:  Entomology  and  Nematology 

Ecological  and  physiological  influences  on  the  chironomid  midge,  Glyptotendipes 

paripes,  were  examined  to  develop  a  qualitative  spatial  and  temporal  computer  model  of 

larval  distributions  in  central  Florida  lakes.  Larval  population  distributions  of  G.  paripes 

were  examined  over  a  two-year  period  in  three  central  Florida  lakes  (Dora,  Yale,  and 

Wauburg)  in  relation  to  selected  water  and  sediment  physico-chemical  properties.  Lake 

bathymetry  and  sediment  physical  compositions  were  extensively  mapped  in  each  lake. 

Composition  of  midge  larvae  in  the  gut  contents  offish  in  relation  to  the  standing  crop  of 

midge  larvae  in  Lakes  Dora  and  Yale  was  studied.  Development  times  for  immature 

G.  paripes  were  determined  in  the  laboratory  by  rearing  under  a  series  of  constant 

temperatures.  These  data  were  used  to  construct  an  exploratory  spatio-temporal 

computer  model  of  G.  paripes  larval  distributions  in  central  Florida  lakes. 
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Canonical  correspondence  analysis  indicated  that  the  primary,  significant 
influences  on  distributions  of  G.  paripes  larvae  in  the  study  lakes  were  water  depth  and 
the  following  physico-chemical  sediment  characteristics:  dry  weight,  total  organic  carbon 
and  presence  of  sand,  muck  or  vegetation.  While  G  paripes  larval  distributions  tended  to 
be  associated  with  shallower,  sandy  substrates  in  Lakes  Dora  and  Wauburg,  these  larvae 
were  distributed  over/in  muck  bottom  areas  of  Lake  Yale.  The  development 
requirements  of  this  species  were  calculated  to  be  716  degree-days  above  9.0  °C  from 
egg  to  emergence.  Gut  content  analysis  indicated  that  bluegill  (Lepomis  macrochirus) 
were  indiscriminate  feeders  on  midge  larvae  in  shallower  portions  of  the  lakes  and  larvae 
in  deeper  portions  of  the  lakes  perhaps  were  not  within  the  grazing  range  of  the  fish. 
When  validated  against  previously  collected  data  from  Lakes  Jesup  and  Monroe,  the 
computer  model  successfully  estimated  G.  paripes  larval  distributions.  The  final  model 
was  effective  at  estimating  larval  distributions  in  Lakes  Dora,  Yale  and  Wauburg  with  the 
exception  of  the  field  data  collected  from  Lakes  Dora  and  Yale  during  the  latter  part  of 
the  field  study.  Further  work  is  needed  to  determine  and  understand  the  environmental 
conditions  that  produced  the  unusual  changes  in  G  paripes  distributions  noted  in  these 
lakes  that  resulted  in  model  predictions  deviating  from  field  data.  Theoretical 
manipulations  of  conditions  within  the  model  produced  results  that  indicated  the  model  to 
be  useful  in  preliminary  ecological  studies  by  providing  initial  hypotheses  on  G  paripes 
larval  population  responses  to  selected  environmental  conditions  studies. 
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CHAPTER  1 
GENERAL  INTRODUCTION 


Non-biting  aquatic  insects  of  the  family  Chironomidae  (Order:  Diptera)  can  cause 
nuisance  because  of  their  large  and  frequently  emanating  adult  swarms  (Ali  1995)  and 
associated  medical  problems  (Cranston  1995)  for  waterfront  residents,  visitors  and 
business  owners.  This  is  especially  true  in  Florida,  which  is  laced  with  numerous 
widespread  aquatic  ecosystems  in  urban  and  suburban  areas  and  significant  waterfront 
developments.  An  estimated  >1 0,000  species  of  Chironomidae  are  distributed 
worldwide,  with  >4,000  species  identified  globally  and  thus  far  280  species  in  1 19  genera 
identified  from  Florida  (Epler  1995).  Worldwide,  100  or  more  species  of  the  subfamilies 
Chironominae,  Orthocladiinae,  and  Tanypodinae  have  been  reported  to  emerge  at 
nuisance  levels  (Ali  1995).  In  Florida,  Beck  and  Beck  (1969)  identified  Glyptotendipes 
paripes,  Chironomus  crassicaudatus,  Chironomus  decorus,  Chironomus  fulvipilus, 
Chironomus  carus,  Chironomus  stigmaterus,  and  Glyptotendipes  lobiferus  occurring  at 
nuisance  levels.  In  central  Florida,  Ali  and  Fowler  (1985)  identified  G.  paripes  and 
C.  crassicaudatus  as  the  dominant  pest  species.  The  chironomid  nuisance  results  from 
the  rather  high  productivity  of  these  organisms  in  some  habitats.  Most  nuisance  species 
generally  have  short  adult  lifespans,  so  accumulations  of  large  numbers  of  dead  and 
odorous  midges  add  to  the  midge  nuisance  problems.  Adult  chironomids  can  soil 
buildings  and  equipment  with  egg  masses  and  meconium,  foul  sensitive  outdoor 


equipment,  and  in  extreme  situations  accumulate  on  roads  in  quantities  that  produce 
slippery  driving  conditions  (Ali  1995).  An  economic  impact  study  commissioned  by  the 
Greater  Sanford  Chamber  of  Commerce,  Sanford,  FL  (Anonymous  1977)  revealed  that 
midges  emerging  from  Lake  Monroe  (a  part  of  the  St.  Johns  River  in  Seminole  and 
Volusia  Counties,  FL)  and  other  nearby  aquatic  systems  resulted  in  a  loss  of  3-4  million 
dollars  annually  for  lakefront  businesses  of  the  city  of  Sanford.  This  study  also  reported 
that  at  least  10  counties  in  Florida  suffered  economic  losses  because  of  midge  swarms. 

Medical  problems  associated  with  midges  primarily  involve  allergies 
(conjunctivitis,  rhinitis,  hay  fever,  and  asthma)  triggered  by  inhalation  of  adult  midge 
fragments  or  skin  contact  with  the  larval  hemoglobins  (Cranston  1995).  Presently, 
midges  are  not  known  to  be  vectors  of  any  disease  organism.  However,  a  recent  report 
has  identified  Chironomus  sp.  egg  masses  from  waste-stabilization  lagoons  as  a  reservoir 
for  non-pathogenic  biotypes  of  Vibrio  cholerae,  the  cholera  pathogen,  and  possibly  a 
reservoir  for  pathogenic  biotypes  (Broza  and  Halpern  2001). 

Other  problems  associated  with  larval  midges  include  contamination  of  drinking 
water,  and  their  presence  as  agricultural  pests  on  rice,  Japanese  horseradish,  and  several 
ornamental  plants  (Ali  1995,  Cranston  1987,  Stevens  and  Warren  1994).  The  various 
problems  are  often  exacerbated  by  cultural  eutrophication  of  adjacent  aquatic  habitats, 
creating  favorable  conditions  for  midge  production  in  large  numbers.  In  many  situations, 
midge  problems  intensify  due  to  creation  of  new  midge  habitats,  such  as  stormwater 
holding  ponds,  effluent  discharge  channels,  and  residential-recreational  lakes  amid 
human  population  centers. 


Despite  the  problems  associated  with  adult  and  larval  chironomids,  midge  larvae 
and  pupae  are  an  important  element  of  most  aquatic  food  webs,  often  being  the  dominant 
or  near  dominant  organisms  in  the  benthos,  providing  food  for  many  fish  and 
macroinvertebrates  (Armitage  1995,  Tokeshi  1995).  Therefore,  while  developing  any 
midge  control  strategies,  care  must  be  taken  to  ameliorate  the  problems  for  humans 
without  having  any  serious  or  irreversible  effects  on  the  food-chain  and  overall  ecology 
of  the  aquatic  habitat  from  which  they  emerge. 

Within  the  framework  of  developing  an  integrated  population  management 
strategy  for  midges,  a  spatio-temporal  model  for  chironomid  larval  populations  seems 
necessary.  Although  a  few  lakes  in  central  Florida  (e.g.,  eutrophic  Lakes  Monroe  and 
Jesup)  have  been  thoroughly  investigated  for  chironomid  ecology  and  productivity  in  the 
past  two  decades  (Ali  and  Baggs  1982,  AH  and  Alam  1996,  Ali  et  al.  1996),  more 
information  beyond  the  present  database,  including  hypereutrophic  and  mesotrophic 
habitats  in  central  Florida,  was  needed.  The  objective  of  the  present  study  was  to  collect 
relevant  data  from  additional  central  Florida  lakes  to  assemble  a  spatial  and  temporal 
model  of  G.  paripes  larval  population  distributions.  Such  a  model  would  enable  lake 
managers  to  reduce  sampling  efforts  and  selectively  target  control  measures  to  specific 
areas  of  a  habitat  with  high  larval  densities,  and  only  at  the  times  when  they  are  abundant. 
This  would  result  in  control  measures  interfering  in  only  a  relatively  small,  partial  area  of 
a  waterbody.  Successful  management  efforts  that  reduce  adult  populations  will  have 
considerable  benefits,  such  as  reduced  nuisance  problems  with  associated  costs  and 
efforts  to  clean  and  maintain  lakefront  facilities,  reduced  economic  losses  to  water-related 
recreation  and  tourism  industries,  and  overall  a  better  quality  of  life  for  Florida  citizens. 


Based  upon  the  Sanford  Chamber  of  Commerce  estimates  (Anonymous  1977),  these 
savings  could  be  in  the  millions  of  dollars  statewide. 

To  develop  a  spatio-temporal  model  of  G.  paripes  larval  populations  in  central 
Florida  lakes,  the  following  specific  objectives  were  required.  Detailed  bathymetric  and 
sediment  characteristic  maps  of  the  three  study  lakes  (Lakes  Dora,  Yale  and  Wauburg) 
had  to  be  developed.  Habitat  specific  physico-chemical  characteristics  influencing 
chironomid  larval  populations  were  elucidated.  Predation  of  chironomid  larvae  by  fish 
was  analyzed  by  examining  fish  gut  contents  in  relation  to  concurrent  standing  crop  of 
midge  larvae.  Temperature-related  effects  on  G  paripes  development  were  determined 
and  degree-day  requirements  for  this  species  were  estimated  in  the  laboratory.  These 
data  were  used  to  develop  and  calibrate  a  model  using  the  matrix  techniques  first 
described  by  Allen  et  al.  (1996).  This  model  was  qualitatively  validated  using  G.  paripes 
larval  population  and  physico-chemical  environmental  data  from  Lakes  Jesup  and 
Monroe.  The  utility  of  the  model  as  a  theoretical  tool  in  developing  new  hypotheses  was 
tested  by  varying  different  parameters  within  the  model. 


CHAPTER  2 
REVIEW  OF  RELEVANT  LITERATURE 


In  the  past  several  decades,  numerous  investigations  concerning  bionomics  and 
management  of  chironomid  populations  have  been  undertaken  in  many  countries  of  the 
world  as  reviewed  by  Ali  (1995).  For  population  management  purposes,  a  variety  of 
physical,  cultural,  biological  and  chemical  means  have  been  employed,  with  a  majority  of 
attempts  focused  on  chemical  control.  Predators  of  midges  have  been  investigated,  and 
many  parasites  and  pathogens  have  been  identified.  Commercially-produced  bacteria, 
Bacillus  thuringiensis  serovar.  israelensis  and  Bacillus  sphaericus,  have  been  used  in 
attempts  to  reduce  midge  larvae  in  some  habitats.  Physical  and  cultural  control 
techniques  have  included  mechanical  means,  habitat  management  and  ecological 
manipulations,  and  behavioral  manipulations  of  adult  midge  populations.  Details  of 
midge  bionomics  and  management  studies  on  a  global  basis,  and  in  Florida,  were 
summarized  by  Ali  (1991,  1996  a,  b). 

In  Florida,  Provost  (1957,  1958,  1959)  provided  information  on  productivity,  food 
habits,  and  nutrient  interactions  of  chironomids  in  13  lakes  in  the  Winter  Haven/Polk 
County  area.  Glyptotendipes  paripes  was  the  predominant  chironomid  in  most  of  these 
lakes,  with  a  preference  for  sand  (75%  of  collected  larvae)  or  peat/sand  (23%  of  collected 
larvae)  bottoms.  Larval  densities  of  G.  paripes  up  to  4,730/m2  were  reported,  and  such 
high  densities  generally  existed  under  low  water  level  conditions.  This  species  was 


absent  from  lakes  with  minimal  human  influences,  as  well  as  from  pure  muck  bottoms. 
Gut  analysis  of  G  paripes  larvae  revealed  that  96%  of  the  total  diet  consisted  of  blue- 
green  algae,  with  the  balance  consisting  of  green  algae  (2%),  diatoms  and  other  materials. 
It  was  concluded  that  feeding  of  G.  paripes  larvae  was  non-selective,  since  blue-green 
algae  were  by  far  the  most  common  algae  in  the  feeding  zone  of  the  larvae.  Nielsen 
(1962)  studied  the  life  history  of  imagines  (adults)  of  G  paripes  in  Lake  Cannon,  Winter 
Haven,  FL  and  reported  the  unusual  dawn  mating  swarms  of  this  species,  and  the  peak 
period  of  adult  emergence  within  45-120  minutes  after  sunset.  More  than  3  decades  ago, 
midge  larval  and  adult  control  by  chemical  means  (primarily  organophosphorus 
insecticides)  in  some  central  Florida  situations  was  studied  by  Patterson  (1964,  1965)  and 
Patterson  and  Wilson  (1966)  who  reported  an  excellent  control  of  G  paripes  larvae  with 
temephos  used  at  rates  as  low  as  0.06  kg  active  ingredient  (AI)/ha. 

In  another  situation  in  Florida,  Lake  Thonotosassa,  Hillsborough  County,  was 
reported  to  be  highly  eutrophic,  with  chironomid  fauna  dominated  by  C.  crassicaudatus 
(31%),  G  paripes  (29%),  and  Polypedilum  halter  ale  (13%).  The  latter  two  species 
preferred  sand  bottom  areas,  with  P.  halterale  occurrence  limited  to  mostly  shallow 
water.  Chironomus  crassicaudatus  was  distributed  throughout  the  lake,  except  during 
summer  when  it  was  absent  from  the  deep  areas  of  the  lake  (Cowell  et  al.  1975).  Later, 
Cowell  and  Vodopich  (1981)  reported  the  distribution  and  seasonal  abundance  of 
Chironomidae  and  other  invertebrates  in  Lake  Thonotosassa.  Chironomids  comprised 
37.1%  of  total  invertebrates  collected.  The  dominant  chironomid  species  was 
P.  halterale,  followed  by  G  paripes  and  C.  crassicaudatus.  Seasonal  peaks  of 
abundance  for  P.  halterale  were  spring  and  summer,  for  G  paripes  summer  and  fall,  and 


spring  for  C.  crassicaudatus.  Life  cycles  of  these  species  ranged  between  14  and  22  days 
and  were  multivoltine.  Fuller  and  Cowell  (1985)  studied  the  invertebrate  community 
recovery  with  special  reference  to  Chironomidae  in  Lake  Thonotosassa  after  disturbances 
caused  by  nesting  of  the  cichlid  fish,  Saratherodon  aurea. 

Benthos  of  three  eutrophic  lakes  (Apopka,  Dora,  and  Griffin)  in  the  Ocklawaha 
chain  of  lakes,  Florida,  were  investigated  for  three  consecutive  years  by  Holcomb  et  al. 
(1974,  1975,  1976).  These  authors,  limiting  their  identification  of  chironomids  to  family 
level,  reported  annual  mean  larval  densities  of  641,  364,  and  764/m2  in  Lake  Apopka; 
2,094,  9,249,  and  6,012/m2  in  Lake  Dora;  and  1,793,  7,896,  and  5,546/m2  in  Lake  Griffin 
during  1973,  1974,  and  1975,  respectively.  In  Lake  Apopka,  midge  density  declined  with 
water  depth,  whereas  in  Lake  Dora,  low  densities  occurred  at  shallow  contours, 
increasing  at  the  3-ft  (0.91  m)  depth,  remaining  fairly  constant  down  to  1 1 -ft  (3.35  m), 
and  declining  significantly  below  1 1 -ft  (3.35  m)  water  depths.  In  Lake  Griffin,  however, 
the  midge  larval  densities  gradually  increased  with  water  depth.  The  authors  attributed 
these  midge  larval  distribution  trends  to  what  they  described  as  "poor  quality  of  deep 
water"  in  Lake  Apopka,  and  "gradually  improving  quality"  further  downstream  (deep 
water)  in  the  chain  of  lakes.  In  addition,  heavy  predation  by  fish  in  the  shallow  areas  of 
Lake  Griffin  was  speculated  to  limit  midge  densities  at  these  depths.  A  U.  S. 
Environmental  Protection  Agency  survey  (EPA  1978)  included  Lakes  Apopka,  Dora, 
Griffin,  and  Yale  in  the  National  Eutrophication  Survey  and  classified  the  status  of  these 
lakes  as  eutrophic  or  hypereutrophic.  Canfield  (1981)  studied  chemistry  and  trophic  state 
of  lakes  in  different  geographical  regions  of  Florida  and  included  six  lakes  of  the 
Ocklawaha  chain  in  his  investigation.  He  reported  that  overall  the  lakes  in  the 


Ocklawaha  chain  were  eutrophic  hard-water  lakes,  except  for  Lake  Yale,  which  was 
classified  as  mesotrophic,  a  classification  that  differed  from  that  of  the  earlier 
classification  by  the  EPA  for  Lake  Yale  (EPA  1978).  Preston  (1983)  presented  a  large 
amount  of  historical  data  and  limnological  data  based  on  a  4-year  study  of  five  lakes 
(Apopka,  Beauclair,  Dora,  Eustis,  and  Griffin)  in  the  Ocklawaha  chain  reported  by 
Brezonik  et  al.  (1978,  1981),  Tuschall  et  al.  (1979),  and  Pollman  et  al.  (1980)  in  four 
annual  reports.  These  authors  reported  the  relatively  poorest  (compared  to  other  lakes  in 
the  study)  water  quality  was  in  Lakes  Apopka  and  Griffin,  which  received  considerable 
farm  runoff  from  nearby  agricultural  areas,  followed  by  Lakes  Beauclair,  Dora,  and 
Eustis.  These  reports  described  declining  overall  water  quality  and  increasing 
eutrophication  in  all  lakes  during  1977-1980  [Brezonik  et  al.  (1978,  1981),  Tuschall  et  al. 
(1979),  and  Pollman  et  al.  (1980)]. 

In  1987,  Callahan  and  Morris,  though  surveying  primarily  for  Culicidae  in  some 
Polk  County  lakes,  reported  that  G  paripes  was  the  predominant  midge  species  in  these 
lakes,  preferring  sand  bottom.  These  authors  reported  G.  paripes  mean  emergence 
densities  of  108,625  and  150,012  adults/ha/day  from  lakes  rated  "expected  good"  and 
"expected  fair",  respectively,  for  production  of  the  Coquillettidia  mosquitoes.  Butler  et 
al.  (1992)  studied  benthic  invertebrates  in  Lake  Tohopekaliga,  Osceola  County,  for  two 
years  after  an  extreme  drawdown  and  organic  substrate  removal  program.  Among 
chironomids,  Glyptotendipes  sp.  I,  C.  stigmaterus,  Procladius  sp.,  and  Cryptochironomus 
fulvus,  in  that  order,  were  abundant  during  the  first  year.  In  the  second  year,  a 
considerable  numerical  decline  of  most  benthic  invertebrates  collected  was  noted  in  the 
restored  areas.  Chironomus  stigmaterus  and  Procladius  sp.  were  the  most  abundant 
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chironomids  at  that  time,  though  these  numbers  were  over  ten  fold  lower  when  compared 
to  the  first  year  after  restoration.  The  authors  suggested  that  the  improved  aquatic 
macrophyte  diversity  in  the  restored  areas  attracted  a  proportionately  greater  population 
of  predatory  fish  that  preyed  upon  the  benthic  macroinvertebrates,  resulting  in  the  overall 
lower  densities  during  the  second  year  after  restoration.  Warren  et  al.  (1995)  studied  the 
sublittoral  benthic  communities  of  highly  eutrophic  Lake  Okeechobee.  Chironomidae 
were  found  to  compose  25.9%  of  the  total  benthic  organisms,  with  Cladotany  tarsus  sp. 
being  the  most  abundant  (7.9%  of  total  benthic  organisms),  followed  by  C. 
crassicaudatus,  P.  halterale,  and  Djalmabatista  pulcher  forming  7.3,  3.0,  and  2.2%  of  the 
total  benthic  organisms,  respectively.  Seasonal  mean  larval  densities  during  February, 
May,  August,  and  November,  respectively,  were:  889.  616,  3,101,  and  600/m2  for 
Cladotany  tarsus  sp.,  10,  35,  1,621,  and  157/m2  for  C  crassicaudatus,  and  478,  111, 
1,056,  and  803/m2  for  P.  halterale. 

During  the  past  two  decades,  chironomid  larval  and  adult  populations  and  the 
physical  and  chemical  limnology  of  a  few  lakes  in  the  St.  Johns  River  system,  Florida, 
have  been  investigated  and  reported  in  the  literature.  Specifically,  Lakes  Monroe 
(Seminole  and  Volusia  Counties)  and  Jesup  (Seminole  County)  have  been  the  focus  of 
these  investigations  among  lentic  ecosystems,  and  in  lotic  habitats,  tributaries  of  the 
Wekiva  River.  Ali  and  Baggs  (1982)  reported  seasonal  changes  of  chironomid 
populations  in  Lake  Monroe.  Glyptotendipes  paripes  and  C.  crassicaudatus  were 
quantitatively  important,  with  larval  densities  ranging  from  <100  to  6,000/m2  and  <10  to 
1,800/m2,  respectively.  Adults  of  these  species  dispersed  in  significant  numbers  up  to 
200  m  from  the  lake,  with  C.  crassicaudatus  showing  a  relatively  greater  inland  dispersal. 
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At  lakefront,  up  to  350,000  adult  midges  occurred  in  a  New  Jersey  light  trap/night  (Ali 
and  Fowler  1983).  Ali  et  al.  (1983)  reported  that  the  two  dominant  midge  species  in 
Lake  Monroe  emerged  on  an  almost  daily  basis  at  water  temperatures  between  1 8-30  °C, 
resulting  in  numerous  asynchronous  generations  of  these  midges  in  a  year.  Patterns  of 
daily  adult  abundance  of  these  midges  were  documented  and  occurrence  of  oviposition 
on  an  almost  daily  basis  was  suggested  (Ali  et  al.  1985).  Ali  and  Fowler  (1985) 
presented  data  on  midge  populations  in  Lake  Monroe  between  1 979  and  1 984.  A  general 
declining  trend  of  G.  paripes  and  C.  crassicaudatus  populations  was  noted  with 
reductions  amounting  up  to  93%  between  years.  Several  prominent  peaks  of  adult 
emergence  occurred  between  May  and  November.  Ali  et  al.  (1996)  presented  further 
adult  chironomid  data  based  on  New  Jersey  light  trap  collections  made  from  waterfront 
land  areas  of  Lake  Monroe  covering  the  1980-1994  period.  Total  number  of  adults 
fluctuated  from  annual  mean  values  of  269  adults/trap/day  in  1994  to  8,009 
adults/trap/day  in  1980.  Glyptotendipes  paripes  composed  69.5%  of  total  chironomids 
collected  over  the  study  period,  followed  by  C.  crassicaudatus  (26.1%).  The  latter 
species  was  most  abundant  in  late  spring  and  early  summer,  while  G.  paripes  was  most 
abundant  in  late  summer.  Overall,  strong  positive  correlation  of  midge  populations  with 
monthly  mean  air  temperatures  was  noted,  whereas  yearly  densities  were  significantly 
influenced  by  water  levels  in  the  lake  as  well  as  by  annual  rainfall  in  the  area.  Ali  et  al. 
(1988)  presented  data  on  water  and  sediment  physico-chemical  parameters  in  Lake 
Monroe.  Specific  conductance  of  up  to  1,832  pS/cm  indicated  high  levels  of  dissolved 
solids.  Mean  values  of  total  N  and  total  P  were  1.82  and  0.21  mg/L,  respectively.  No 
clear  seasonal  trends  were  noted.  Organic  matter  in  the  sediments  ranged  from  1  to  1 82  g 


C/kg,  indicating  wide  range  of  spatial  distribution.  All  measured  sediment  nutrients  were 
significantly  (P  <  0.01)  correlated  with  sediment  organic  carbon  content.  Little  or  no 
significant  spatial  or  temporal  relationships  existed  between  the  studied  water  and 
sediment  parameters,  though  periodic  hydrologic  flushing  of  the  lake  by  the  St.  Johns 
River  may  have  obscured  some  correlations.  Ali  and  Alam  (1996)  provided  quantitative 
estimates  of  nutrients  at  the  sediment- water  interface  in  Lake  Jesup.  Increases  of 
turbidity  due  to  wind  actions  or  other  factors  had  a  positive  influence  on  soluble  reactive 
P  (SRP)  concentrations  in  water  at  the  sediment-water  interface  and,  in  turn,  caused 
increases  of  chlorophyll  a  concentrations  in  water.  Mean  concentrations  of  Kjeldahl 
nitrogen,  extractable  SRP,  and  total  organic  carbon  in  lake  sediments  were  10.9,  0.019 
and  65.5  mg/g,  respectively,  during  the  March  1988  to  December  1993  period  of 
investigation.  Such  levels  of  concentrations  of  these  chemical  characteristics  are 
indicative  of  the  eutrophic  status  of  the  lake.  Ali  and  Lobinske  (personal  communication) 
recently  completed  some  detailed  investigations  of  chironomid  larval  densities  and 
selected  water  and  sediment  physico-chemical  parameters  in  Lakes  Monroe  and  Jesup, 
and  in  Lake  Monroe  and  a  water-cooling  reservoir  (Ali  et  al.  2001).  Preliminary  analysis 
of  these  data  indicates  a  strong  relationship  between  larval  spatial  distributions  and 
sediment  type  as  well  as  water  depth.  Glyptotendipes  paripes  had  the  highest  occurrence 
frequency  in  sand  although  highest  larval  densities  of  this  species  were  found  in  slightly 
looser  substrates,  such  as  muck  or  a  muck/sand  combination.  This  species  was  generally 
negatively  correlated  with  water  depth  in  the  two  lakes.  Xue  and  Ali  (1997)  in  a 
laboratory  study  described  the  case-making  behavior  of  G.  paripes  larvae  and  the  larval 
preference  for  several  ranges  of  substrate  particle  size.  The  most  successful  case  making 
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was  reported  for  two  sand  particle  ranges,  0.35  and  0.58  mm  (100%),  and  0.59  and  0.84 
mm  (98%),  indicating  the  suitability  of  such  substrates  (particle  size  range)  for 
colonization  and  propagation  of  G.  paripes  in  field  situations. 

Recently,  Lobinske  etal.  (1996,  1997)  studied  lotic  chironomids  and  other 
macroinvertebrates  in  relation  to  prevailing  selected  physico-chemical  parameters  in 
Blackwater  Creek  and  Rock  Springs  Run,  tributaries  of  the  Wekiva  River,  central 
Florida.  Chironomidae  were  numerically  the  dominant  macroinvertebrates,  comprising 
40.9%  and  34.8%  of  total  macroinvertebrates  in  Blackwater  Creek  and  Rock  Springs 
Run,  respectively.  Using  the  Hilsenhoff  biotic  index  (Hilsenhoff  1977,  1982,  1987), 
these  streams  were  rated  at  fair  to  good  water  quality  (Lobinske  et  al.  1997).  Lobinske  et 
al.  (1996)  reported  24  genera  of  chironomids  in  Blackwater  Creek  and  26  in  Rock 
Springs  Run,  with  Tanytarsus  spp.  numerically  the  most  common  in  the  former  stream 
(29.7%)  and  Polypedilum  spp.  the  most  common  in  the  latter  (21 .4%).  Chironomid 
productivity  in  both  streams  was  1.1  ±  0.8  g  dry  weight/m2,  indicating  their  oligotrophic 
status.  Among  the  parameters  studied,  water  level  in  the  streams  was  the  dominant 
influence  on  chironomid  population  fluctuations  in  the  two  streams  during  the  study 
period  (Lobinske  et  al.  1996). 

The  existing  literature  did  not  provide  sufficient  specific  information  to  develop  a 
spatio-temporal  larval  distribution  model  of  G.  paripes.  There  were  no  available 
quantified  relationships  of  G  paripes  larvae  with  water  and  sediment  characteristics  in 
Florida  lakes  that  could  be  readily  used.  Previous  results  showing  the  general  trend  of 
larval  G  paripes  associated  mostly  with  firm  sediments  did  not  provide  specific 
information  on  the  relative  strength  of  this  trend  or  what  kind  of  transition  there  was  in 
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population  levels  between  firm  and  soft  sediments.  While  there  were  some  detailed 
bathymetric  maps  available  for  lakes  in  central  Florida,  the  only  maps  available  of 
sediment  characteristics  were  for  soft  sediment  depth  (Danek  et  al.  1991).  Previously, 
development  time  of  immature  G.  paripes  had  been  only  crudely  estimated  and  the  lower 
development  threshold  temperature  had  not  been  determined.  While  there  was  general 
information  on  fish  predation  of  chironomid  larvae,  species  specific  information  on 
chironomids  consumed  was  not  available. 


CHAPTER  3 
DESCRIPTION  OF  STUDY  LAKES 


Lakes  Dora,  Yale  and  Wauburg  are  located  in  the  Central  Valley  region  of  central 
Florida  in  the  lowlands  between  the  Mt.  Dora  Ridge  and  the  Ocala  Uplift  District  (Brooks 
1982).  These  lakes  were  selected  to  obtain  additional  information  on  G.  pctripes 
population  ecology  for  developing  the  model.  The  Ocklawaha  chain  of  lakes  in  Lake 
County  is  a  major  geographic  component  of  the  area,  with  eight  major  lakes  (including 
Dora  and  Yale)  in  the  chain  having  different  water  qualities,  from  hypereutrophic 
(Apopka,  Beauclair)  to  mesotrophic  (Yale,  Weir)  (Fulton  1995a,  b).  These  lakes  share  a 
common  watershed  complex  and  hydrologic  base  (Fulton  1995b).  Several  water  control 
structures  in  the  basin  limit  natural  water  fluctuations,  so  the  systems  are  like  managed 
reservoirs  (Fulton  1995b).  The  northern  portion  of  the  Central  Valley  district  extending 
through  Marion  County  into  Alachua  County  contains  higher  numbers  of  softwater  lakes 
(including  Lake  Wauburg)  compared  to  lakes  in  the  Ocklawaha  system. 

Lake  Dora  covers  an  area  of  1 ,787  ha  with  a  mean  water  depth  of  3.0  m  and 
maximum  depth  of  5.2  m.  Mean  sediment  depth  was  1 .4  m,  with  the  unconsolidated  floe 
averaging  0.33  m;  soft  sediments  covered  85.0%  of  the  lake  bottom  (Danek  et  al.  1991). 
Lake  Dora  was  characterized  as  hypereutrophic  by  the  National  Eutrophication  Survey 
(EPA  1978).  Median  conductivity  value  amounted  to  320  uS/cm,  Secchi  disk 
transparency  0.4  m,  total  phosphorus  (TP)  0.102  mg/L,  ortho-phosphorus  (OP)  0.022 
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mg/L,  inorganic  nitrogen  (IN)  0.240  mg/L,  total  nitrogen  (TN)  3.290  mg/L  and 
chlorophyll  a  60.0  ug/L.  Mean  physico-chemical  values  for  Lake  Dora  presented  by 
Canfield  (1981)  were  as  follows:  pH  8.9,  conductivity  321  uS/cm,  TN  3.100  mg/L,  TP 
0.903  mg/L,  chlorophyll  a  124  ug/L,  and  Secchi  disk  transparency  0.4  m.  Mean  values 
for  the  east  basin  of  Lake  Dora  from  June  1990  to  November  1996,  published  by  the 
Florida  LAKEWATCH  program,  were  pH  8.7,  conductance  416  uS/cm,  TP  0.092  mg/L, 
TN  3.488  mg/L,  chlorophyll  (not  specified  if  chlorophyll  a  or  total)  156  ug/L  and  Secchi 
disk  transparency  0.4  m  (Florida  LAKEWATCH,  1 997).  Aquatic  plants  covered  a  mean 
of  4%  of  the  lake  area,  with  mean  emergent  plant  biomass  of  4.6  kg  wet  wt/m2,  floating 
leafed  biomass  0.0  kg  wet  wt/m2,  and  submerged  plant  biomass  0.3  kg  wet  wt/m2  (July 
1992  data,  Florida  LAKEWATCH  1997).  Mean  values  for  the  west  basin  of  Lake  Dora 
from  June  1990  to  November  1996  were  pH  8.7,  conductance  404  uS/cm,  TP  0.073 
mg/L,  TN  3.464  mg/L,  chlorophyll  144  ug/L  and  Secchi  disk  transparency  0.4  m. 
Aquatic  plant  cover  amounted  to  4%,  with  3.3  kg  wet  wt/m2  emergent,  0.2  kg  wet  wt/m2 
floating  leafed  and  1.5  kg  wet  wt/m2  submersed  plant  biomass  (July  1992  data,  Florida 
LAKEWATCH  1997). 

Lake  Yale  covers  1,632  ha,  with  a  mean  water  depth  of  3.7  m  and  maximum 
depth  of  7.9  m.  Mean  sediment  depth  was  1 .9  m,  with  unconsolidated  floe  depth  of 
0.18  m.  Soft  sediments  covered  88.9%  of  the  lake  bottom  (Danek  e/  al.  1991).  Lake 
Yale  was  classified  as  eutrophic  by  the  EPA  (1978),  with  conductivity  amounting  to  295 
uS/cm,  Secchi  disk  transparency  1.5  m,  TP  0.027  mg/L,  OP  0.014  mg/L,  IN  0.160  mg/L, 
TN  1.270  mg/L  and  chlorophyll  a  25  A  ug/L.  According  to  Canfield  (1981),  Lake  Yale 
had  subjectively  better  water  quality  than  that  reported  by  the  EPA  (1978)  as  indicated  by 
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the  following  values:  pH  8.3,  conductivity  264  uS/cm,  TN  0.683  mg/L,  TP  0.014  mg/L, 
chlorophyll  a  9.6  ug/L,  and  Secchi  disk  transparency  1 .4  m.  Mean  values  reported  by 
Florida  LAKEWATCH  (1997)  for  Lake  Yale  from  May  1990  to  June  1994  were  pH  8.3, 
conductance  264  uS/cm,  TP  0.01 1  mg/L,  TN  0.805  mg/L,  chlorophyll  7  ug/L  and  Secchi 
disk  transparency  2.4  m.  Canfield  (1981)  considered  Lakes  Dora  and  Yale  as  hardwater 
lakes  with  mean  total  hardness  of  1 15-155  mg/L  (as  CaCC^).  Because  both  of  these  lakes 
tend  to  be  below  the  piezometric  surface,  considerable  amounts  of  water  inputs  to  these 
lakes  are  from  mineralized  groundwater  or  as  run-off  from  calcareous,  nutrient-rich  soils. 

Lake  Wauburg  is  a  much  smaller  lake,  covering  1 50  ha  surface  area  with  mean 
water  depth  of  3.0  m.  The  lake  contains  soft  water  (mean  conductance  79  uS/cm,  total 
alkalinity  19.9  mg/L  as  CaCC^)  and  is  eutrophic  as  supported  by  the  following  values: 
mean  TP  0.107  mg/L,  TN  1.726  mg/L,  chlorophyll  83  ug/L  and  Secchi  disk  transparency 
0.56  m.  These  means  are  based  on  the  March  1990  to  December  1996  data  of  Florida 
LAKEWATCH  (1997). 


CHAPTER  4 

MAPPING  BATHYMETRIC  AND  SELECTED  SEDIMENT  PHYSICO-CHEMICAL 

CHARACTERISTICS  OF  STUDY  LAKES 


Introduction 
The  importance  of  sediment  composition  and  lake  bathymetry  in  determining 
distributions  of  G  paripes  has  been  previously  reported  (Provost  1957,  Curry  1962, 
Cowell  et  al.  1975,  McLachlan  1976,  Milleson  1978,  Cowell  and  Vodopich  1981, 
Callahan  and  Morris  1987,  Ali  et  al.  1998b).  Toward  assisting  our  understanding  of  the 
relation  of  G.  paripes  with  sediment  characteristics,  the  objective  of  this  study  was  to 
produce  systematic  maps  of  basic  sediment  physico-chemical  characteristics  and  lake 
bathymetry.  These  maps  serve  an  important  base  component  in  determining  the  spatial 
distributions  of  larval  populations  in  the  computer  model  and  should  provide  an  effective 
tool  for  lake  managers  targeting  control  measures  for  nuisance  midge  populations. 

Materials  and  Methods 
Systematic  sediment  core  samples  using  the  method  of  Ali  (1984)  were  taken 
from  each  sample  lake.  For  sampling  purposes,  Lakes  Dora  and  Yale  were  mapped  with 
a  grid  interval  of  0.2  minutes  latitude/longitude  and  Lake  Wauburg  with  a  grid  interval  of 
0.05  minutes  latitude/longitude.  In  Lake  Dora,  172  sample  sites  were  surveyed  between 
May  and  September  1998,  while  150  sample  sites  in  Lake  Yale  were  examined  between 
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October  and  December  1998,  and  68  in  Lake  Wauburg  during  February  2000.  At  each 
grid  interval,  three  5-cm  deep  sediment  cores  were  collected,  with  overlying  water 
carefully  decanted  without  disturbing  the  surface  sediments.  Presence  of  sand,  muck, 
clay,  peat,  or  detritus  was  recorded.  The  core  samples  were  composited  in  a  sealed 
polyethylene  bag  and  stored  on  ice  until  transported  the  same  day  to  the  laboratory, 
where  they  were  maintained  at  -10  °C  until  analyzed.  At  each  grid  point  in  Lakes  Dora 
and  Yale,  water  depth  was  measured  with  a  boat-mounted  depth  finder  (Humminbird® 
Model  Wide  100,  Techsonic  Industries,  Inc.  Eufaula,  AL),  while  a  graduated  pole  was 
used  in  Lake  Wauburg  to  measure  water  depth. 

For  various  analyses,  samples  were  thawed  and  transferred  to  labeled  1-L  beakers 
and  thoroughly  mixed.  To  estimate  particle  size  distributions  in  samples  from  Lakes 
Dora  and  Yale,  approximately  1  cm3  of  material  was  transferred  to  a  sieve  series  (1,000, 
350,  250  and  125-um  pore  mesh)  and  material  washed  through  with  deionized  water. 
Material  retained  in  each  sieve  was  carefully  washed  into  labeled,  tared  50-mL  beakers, 
dried  for  24  hours  at  90  °C,  and  dry  material  weighed.  This  analysis  was  not  conducted 
with  samples  from  Lake  Wauburg.  One  gram  wet-sediment  from  each  sample  was 
placed  into  a  tared,  labeled  50-mL  beaker  and  dried  for  24  hours  at  90  °C  to  determine 
overall  sediment  dry  weight.  Sediment  pH  of  wet  sediment  was  determined  using  a 
Model  710A  pH/ISE  meter  (ATI  Orion  Laboratory  Products,  Boston,  MA).  The 
remaining  material  was  air  dried  at  room  temperature  and  ground  with  a  mortar  and 
pestle  to  pass  through  a  350-um  pore  sieve  and  partly  utilized  to  determine  total  organic 
carbon  (TOC)  using  the  wet  oxidation  method  of  Nelson  and  Sommers  (1982)  to  assess 
relative  organic  matter  content  of  the  sediment. 
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Bathymetric  maps  were  generated  with  the  software  Slide  Write®  for  Windows 
Version  5.0  (Advanced  Graphics  Software  Inc.,  Encinitas,  CA)  using  the  study  data,  as 
well  as  the  maps  of  Danek  et  al.  (1991)  for  Lakes  Dora  and  Yale,  and  data  provided  by 
Florida  LAKE  WATCH  (personal  communication)  for  Lake  Wauburg.  Maps  of  sediment 
TOC,  sediment  pH,  and  sediment  dry  weight  were  also  generated  using  Slide  Write®  as 
well  as  a  common  set  of  template  files  for  uniform  appearance.  Sediment  dry  weight 
maps  were  augmented  with  additional  data  collected  from  Lakes  Dora,  Yale,  and 
Wauburg  during  the  ecological  study  of  G.  paripes  detailed  in  the  following  chapter. 
This  software  used  the  Inverse  Distance  method  to  interpolate  values  between  data 
points.  These  were  combined  with  scanned  images  of  lake  boundaries  prepared  from 
United  States  Geologic  Survey  7.5-Minute  Topographic  maps. 

Results  and  Discussion 
Mean  values  of  mapped  parameters  are  shown  in  Table  1  for  the  three  study  lakes. 
Sediment  in  all  three  lakes  was  acidic,  indicating  considerable  microbial  activity,  with  pH 
considerably  lower  in  Lake  Wauburg  (mean  value  6.35  pH  units)  compared  to  Lake  Dora 
(mean  value  6.77  pH  units),  and  Lake  Yale  (mean  value  6.68  pH  units).  The  latter  two 
lakes  were  comparable  to  the  mean  sediment  pH  value  of  6.8  for  Lake  Jesup  reported  by 
Ali  and  Alam  (1996),  but  lower  than  the  7.3-7.5  pH  values  reported  for  four  major 
sediment  types  in  Lake  Apopka  (Rees  et  al.  1996).  Sediments  from  Lake  Dora  had  lower 
TOC  (8.9  %)  and  higher  dry  weight  (41.0  %)  values  than  Yale  (10.6  and  30.6  %, 
respectively)  and  Wauburg  (10.9  and  30.1  %,  respectively).  These  TOC  values  are  near 
the  75  mg/g  (7.5%)  reported  by  Ali  et  al.  (1 988)  for  Lake  Monroe  and  less  than  the 
20.4%  reported  for  Lake  Jesup  (Ali  and  Alam  1996).  Rees  et  al.  (1996)  did  not  give 
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overall  means,  but  reported  TOC  values  of  402  g/kg  (40.2%)  for  peat  sediments,  120  g/kg 
(12.0%)  for  mud,  150  g/kg  (15.0%)  for  littoral  (sand/mud),  and  13  g/kg  (13.0%)  for  sand 
sediments  in  Lake  Apopka.  For  muck  sediments,  the  mean  dry  weight  values  of  10.0% 
(Dora),  10.6%  (Yale),  and  17.4%  (Wauburg)  were  all  greater  than  the  2.9%  (derived  from 
reported  97.1%  moisture  content)  for  unconsolidated  flock  (UCF)  from  Lake  Apopka 
reported  by  Gale  and  Reddy  (1994),  though  some  individual  station  values  in  all  three 
study  lakes  were  comparable.  Mean  water  depth  of  3.5  m  in  Lake  Yale  was  greater  than 
that  of  2.7  m  in  Lake  Dora  and  2.8  m  in  Lake  Wauburg. 

The  deeper  (>3  m  depth)  portions  of  Lake  Dora  were  divided  into  two  basins 
separated  by  a  narrow  neck  of  the  lake  near  28°  47.4'  N,  81°  40.4'  W,  with  the  deepest 
portion  of  the  east  basin  offset  from  center  toward  the  southwest  and  the  deepest  portion 
of  the  west  near  the  narrowest  part  of  the  basin  (Figure  1).  Only  a  small  depression  in  the 
west  basin  and  a  narrow  trench  in  the  east  basin  were  >5  m  in  depth.  An  extensive  area 
of  the  southern  bay  of  Lake  Yale  was  comparatively  shallow  (<3  m  deep)  as  was  the 
northeast  bay  near  28°  55.1 '  N,  81°  43'  W  (Figure  2).  The  deepest  part  of  the  lake 
formed  a  rough  "C"  shape  off  the  southwest  shoreline  and  extended  into  the  center  of  the 
large  northern  basin.  Lake  Yale  exhibited  a  >5  m  deep  trench  in  the  southern  part  of  the 
"C",  with  one  area  >6  m  in  depth.  Lake  Wauburg  (Figure  3)  was  overall  the  shallowest 
of  the  study  lakes,  with  some  scattered  locations  between  4  and  5  m  in  depth,  though  the 
nearshore  areas  in  this  lake  were  considerably  steeper,  with  only  a  narrow  band  <2  m 
deep. 

Sediment  percent  dry  weight  distributions  in  Lake  Dora  exhibited  a  similar  trend 
to  lake  depth  (Figure  4).  Highest  sediment  densities  were  generally  found  in  shallower 
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areas  (<3  m  depth)  with  sand  substrates  that  wave  action  probably  kept  clear  of  soft, 
organic  sediment.  Soft,  low-density  sediments  extended  into  shallower  waters  in  the 
northeast  section  of  the  lake,  the  southeast  section  of  the  west  basin,  and  the  westernmost 
bay.  These  areas  likely  are  more  protected  from  wind  and  experience  less  wave  action, 
or  the  organic  matter  accumulation  is  faster  in  these  areas  than  the  ones  exposed  to  the 
scouring  action  of  waves.  A  band  of  firmer  sediments  in  the  western  portion  of  the  lake 
running  southwest  to  northeast  corresponded  with  a  slight  ridge  in  the  lake  bathymetry. 
Of  interest  is  the  rapid  change  to  firm  sediments  in  the  eastern  basin  just  to  the  south  of 
the  deepest  bathymetric  trench.  Similar  general  trends  were  noted  in  Lake  Yale  (Figure 
5),  with  the  primary  exception  of  the  southern  bay,  which  had  low  density,  soft  sediments 
in  the  shallows  except  for  one  small  area.  A  section  of  firm  sediments  extended  from  the 
southwest  shoreline  into  fairly  deep  water  (4-5  m)  just  to  the  southeast  of  the  deep  trench. 
Firm  sediments  also  extended  from  the  center  of  the  northern  shore,  and  from  the 
northeast  corner.  The  proportion  of  firmer  sediments  was  considerably  lower  in  Lake 
Wauburg  than  in  either  Dora  or  Yale  (Figure  6).  The  largest  concentration  of  firm 
sediments  was  along  the  southwestern  shoreline  and  in  scattered  areas  along  the  south 
and  east  shore.  Pebbles  and  rocks  were  mixed  into  the  sand  sediments  along  the  east 
shore.  A  wide  band  of  soft  sediments  was  encountered  up  the  shoreline  along  the 
northwest  margin  of  the  lake  north  of  29°  81.85'  N.  The  trend  of  shallow  areas  having 
firm/sandy  substrates  and  the  deeper  areas  of  lakes  having  soft  sediments  was  also 
reported  for  Lake  Monroe  (Ali  et  al.  1988)  and  for  Lake  Jesup  (Ali  and  Alam  1996). 
Low  density,  high  water  content  sediments  were  also  reported  for  large  areas  of 
hypereutrophic  Lake  Apopka,  mostly  in  the  form  of  UCF  (Gale  and  Reddy  1994). 
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Sediment  pH  distributions  for  Lakes  Dora,  Yale,  and  Wauburg  are  presented  in 
Figures  7,  8,  and  9,  respectively.  In  the  two  former  lakes,  the  general  trend  was  for  lower 
pH  values  from  soft  sediments  (as  defined  by  lower  percent  dry  weight)  in  the  range  of 
6.2  to  7,  while  firmer  substrates  had  pH  values  generally  between  6.8  and  7.4.  Lake 
Wauburg  had  lower  sediment  pH  values,  with  no  value  over  6.8.  Total  organic  carbon 
distribution  in  Lake  Dora  clearly  demonstrated  an  increase  associated  with  depth,  with 
the  highest  values  noted  in  the  deepest  parts  of  both  basins  (Figure  10).  A  similar  trend 
was  noted  for  Lake  Yale  (Figure  11),  except  for  the  southern  bay,  which  had  elevated 
TOC  values  in  the  shallows,  as  well  as  at  the  tip  of  the  northeastern  bay.  This  southern 
bay  was  flanked  by  wetlands  and  the  sediments  recorded  as  containing  plant  detritus. 
Lake  Wauburg  had  some  exceptionally  high  TOC  values  at  the  center  of  the  west 
shoreline,  northern  tip,  and  southeastern  corner  (Figure  12)  near  observed  wetlands. 
These  areas  were  recorded  as  containing  peat  and  detritus  that  were  likely  of 
allochthonous  origin  from  the  adjacent  wetland  vegetation. 

Sediment  particle  sizes  for  muck  and  sand  sediments  from  Lakes  Dora  and  Yale 
are  shown  in  Table  2.  Particle  sizes  were  not  determined  for  Lake  Wauburg  due  to 
processing  errors.  For  both  muck  and  sand  substrates  from  both  lakes,  the  majority 
(67.1-84.7%)  of  particles  was  retained  by  the  250  and  350-um  pore  sieves.  Significantly 
more  particles  (P  <  0.05)  were  retained  at  250-um  from  both  muck  and  sand  substrates  of 
Lake  Yale,  while  significantly  more  particles  >  1 ,000-um  were  present  in  muck 
sediments,  and  retained  at  125-um  from  sand  sediments  from  Lake  Dora. 

Spatial  distributions  of  particles  are  presented  in  Figure  13  (Lake  Dora)  and 
Figure  14  (Lake  Yale).  In  both  lakes,  the  distributions  of  250  to  350-|am  particles  closely 
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followed  the  areas  of  higher  sediment  dry  weight,  corresponding  with  the  firm,  sand 
substrates.  Particles  between  350  and  1,000-um  were  closely  associated  with  the 
distributions  of  soft  sediments.  Lake  Dora  had  a  wider  distribution  of  large  (>1000-um) 
particles,  especially  in  the  east  basin.  Also,  there  was  an  unusual  concentration  of  fine 
(<125-um)  particles  in  a  small  eastern  bay  of  Lake  Dora. 

All  three  lakes,  with  local  variations,  were  found  to  exhibit  a  similar  trend  in 
bottom  morphology.  Firm,  sandy  substrates  dominated  most  of  the  nearshore  areas  of 
these  lakes,  except  those  adjacent  to  observed  wetlands  that  contained  plant  detritus. 
With  increasing  depth,  soft,  organic  carbon  enriched  sediments  become  dominant,  with 
slightly  lower  pH  values.  These  trends  are  similar  to  those  reported  for  other  central 
Florida  lakes  (Ali  and  Alam  1996,  Ali  et  al.  1988). 
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Table  1 .  Mean  ±  SD  of  selected  sediment  physico-chemical  parameters  and  water  depth 
from  mapping  survey  of  Lakes  Dora  (May-September  1998),  Yale  (October- 
December  1998),  and  Wauburg  (February  2000). 

Lake  Sediment  pH        Total  Organic  Sediment  Dry         Water  Depth 

Carbon  (%)  Weight  (%) (m) 

Dora  6.77  ±0.20  8.9  ±8.3  41.0  ±33.5  2.7  ±1.1 

Yale  6.68  ±0.21  10.6  ±7.9  30.6  ±28.7  3.5  ±1.5 

Wauburg 6.35  ±0.13  10.9  ±8.8  30.1  ±24.0  2.8  ±1.1 


Table  2.  Comparison  of  mean  percent  of  sample  ±  SD  of  sediment  particles  retained  by 
sieve  series  (1,000,  350,  250  and  125-um  pore  mesh)  from  muck  and  sand 
sediments  collected  from  Lake  Dora  (May-September  1998)  and  Lake  Yale 
(October-December  1998). 


Sieve  Pore  Size 

Lake 

1,000-um 

350-um            250-um 

125-um 

>125-um 

Muck 

Dora 

13.6  ±8.4* 

37.4  ±8.1        29.7  ±9.7* 

12.9  ±8.0 

6.5  ±3.9 

Yale 

4.9  ±3.9* 

33.8±  12.2      43.0  ±  11.1* 

Sand 

12.9  ±5.6 

5.4  ±2.9 

Dora 

3.4  ±2.4 

23.7  ±10.0       55.6  ±9.2* 

16.3  ±  11.7* 

1.0±1.3 

Yale 

2.9  ±2.2 

22.3  ±11.1       62.4  ±9.0* 

10.5  ±6.8* 

2.0  ±4.3 

♦Percent  within  particle  size  significantly  (P  <  0.05)  different  between  lakes  by  Student's 
/-test. 
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Figure  13.  Sediment  particle  size  spatial  distributions  in  Lake  Dora 
(Lake  County,  Florida),  May-September  1998. 


38 


s': 


s' : 


1 1 1 1 1 1 1 ll| llll I II ll|l  II  Mill  I |l  II I II I  III  II 1 1 

47  4«  45  44  43  42 


;?: 

250  -  350  um 

58 : 

55 : 

^A^L 

54  z 

»• 

: 

" 

III 

4 

7 

45 

45                   44                   43                  42 

1  tnniMii|iiiiiinniiiiMiii|iiiiminiimim 

47  46  45  44  43  42 


*: 


350  - 1000  um 


■   in  mini  in  iiuiii  n  in  II 1 1  linn  nun  ii  1 1  inn 

47  tl  45  44  43  42 


"     1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1  M  1 1 1 1 1 1 1 1 1 1  1 1 1 1 1 
47  48  45  44  43  42 

Percent  Compostition 


□  <M0 

|  10-20 

20-30 

30-40 

I  40-50 


50-60 
60-70 
70-80 
80-90 


Figure  14.  Sediment  particle  size  spatial  distributions  in  Lake  Yale 
(Lake  County,  Florida),  October-December  1998. 


CHAPTER  5 

LARVAL  DISTRIBUTIONS  OF  Glyptotendipes  paripes  IN  RELATION  TO  WATER 

AND  SEDIMENT  PHYSICO-CHEMICAL  PARAMETERS  IN  STUDY  LAKES 


Introduction 
Availability  of  several  years  of  quantitative  larval  population  data  from  two 
eutrophic  central  Florida  lakes  (Lakes  Jesup  and  Monroe)  (Ali  and  Baggs  1982,  Ali  et  al. 
1998a,  b,  Ali  and  Lobinske,  unpublished  data)  provided  the  initial  basis  for  development 
of  a  distribution  model  of  G  paripes  larvae  in  central  Florida  lakes  and  encouraged 
further  research  into  another  geographic  area  of  the  state,  the  Central  Valley  lake  region 
(Canfield  and  Hoyer  1988).  Previously,  distributions  of  G.  paripes  larvae  have  been 
studied  by  several  authors:  Callahan  and  Morris  (1987),  Cowell  et  al.  (1975),  Cowell  and 
Vodopich  (1981),  Fuller  and  Cowell  (1985),  Milleson  (1978)  and  Provost  (1957).  For 
model  calibration,  two  lakes  in  the  Ocklawaha  chain  of  lakes,  Lake  Dora 
(hypereutrophic)  and  Lake  Yale  (mesotrophic/eutrophic),  and  Lake  Wauburg  (eutrophic) 
in  the  northern  area  of  the  Central  Valley  lake  region  were  included  in  the  investigation 
for  expanding  the  larval  distribution  information  on  this  species  into  different  geographic 
locations.  Chironomid  larvae  and  selected  water  and  sediment  physico-chemical 
parameters  in  these  lakes  were  studied  over  a  two-year  period  using  sampling  plans 
optimized  for  G  paripes  larval  collection  (Cochran  1963).  With  these  data,  the  objective 
of  this  study  was  to  quantify  the  influence  of  water  and  sediment  characteristics  on  spatial 
and  temporal  distributions  of  G  paripes  larvae  in  central  Florida  lakes. 
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Materials  and  Methods 
Preliminary  survey  of  G  paripes  larval  distributions  and  selected  water  and 
sediment  physico-chemical  parameters  in  Lakes  Dora  and  Yale  was  made  in  September 
1997;  Lake  Wauburg  was  surveyed  in  September  1999.  During  the  preliminary  survey  in 
Lakes  Dora  and  Yale,  thirty  stations  were  sampled,  while  1 2  stations  were  sampled  in 
Lake  Wauburg.  Using  a  0.2  precision  level  and  data  from  the  preliminary  surveys,  the 
calculated  sample  size  was  120  (Cochran  1963)  for  the  first  phase  (year)  of  the  study  in 
Lakes  Dora  and  Yale.  For  Lake  Wauburg,  the  calculated  sample  size  was  30.  To 
enhance  sample  statistical  precision,  these  stations  were  randomly  assigned  over  two 
strata  based  on  sediment  dry  weight  (DW)  (stratum  1  <40%  DW  and  stratum  2  >40% 
DW)  and  allocation  optimized  using  the  method  of  Cochran  (1963).  These  strata  were 
chosen  based  on  previous  G  paripes  distributions  from  Lakes  Jesup  and  Monroe  (Ali  et 
al.  1998a,  Ali  and  Lobinske,  unpublished  data).  For  Lake  Dora,  this  allocation  was  22 
stations  in  stratum  1  and  98  in  stratum  2,  while  39  and  81  stations,  respectively,  were 
determined  for  Lake  Yale.  For  Lake  Wauburg,  1 1  and  19  stations,  respectively,  were 
determined.  This  optimized  allocation  increased  the  number  of  samples  collected  from 
the  stratum  with  the  higher  sample  variance  in  the  survey  study.  For  the  first  phase  in 
Lakes  Dora  and  Yale,  samples  were  collected  bimonthly  over  four  consecutive  days  per 
sample  effort,  with  30  stations  sampled  per  day.  After  the  first  year,  allocation  of 
sampling  efforts  was  re-evaluated  using  the  first  year's  data.  This  evaluation  revealed 
that  80  stations  were  required  for  the  second  phase  (year)  of  work,  with  15  stations  in 
stratum  1  and  65  in  stratum  2  for  Lake  Dora,  and  26  and  54  (respectively)  for  Lake  Yale. 
During  the  second  phase,  samples  were  collected  monthly  over  two  days  with  40  stations 
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sampled  per  day.  The  one-year  study  of  Lake  Wauburg  was  conducted  concurrently  with 
the  phase  2  study  of  Lakes  Dora  and  Yale,  with  samples  collected  on  a  monthly  basis. 
For  each  sampling  effort,  paired  latitude  and  longitude  values  for  each  station  were 
randomly  determined  by  a  simple  computer  protocol  operated  within  the  software 
Matlab®,  Version  6.0,  Release  12  (The  Mathworks  Inc.,  Natik,  MA)  (Appendix).  A  7.3- 
m  pontoon  boat  (Lakes  Dora  and  Yale)  or  a  3-m  skiff  (Lake  Wauburg)  was  used  for 
sampling,  with  a  Global  Positioning  System  (GPS)  receiver  without  differential 
correction  (model  GPS  12,  Garmin  International  Inc.,  Olathe,  KS)  used  to  determine  and 
record  each  sampling  station. 

Water  depth  was  measured  at  each  station  as  previously  described  in  Chapter  4. 
In  situ  measurements  of  dissolved  oxygen  (Yellow  Springs  Instruments  Inc.,  Yellow 
Springs,  OH,  model  54A  dissolved  oxygen  meter),  water  temperature  and  conductivity 
(Yellow  Springs  Instruments  Inc.,  Yellow  Springs,  OH,  model  33  salinity-conductivity- 
temperature  meter)  at  the  sediment- water  interface  were  made.  During  phase  1 ,  one 
1 ,000-mL  water  sample  for  chlorophyll  determination  was  collected  at  every  third  station 
from  the  sediment/water  interface  with  a  2.2-L  horizontal  alpha  bottle  (Wildlife  Supply 
Company,  Buffalo,  NY),  stored  on  ice  for  transport  to  the  laboratory  and  at  4  °C  in  a 
refrigerator  until  analyzed  within  24  hours.  Water  samples  for  chlorophyll  determination 
in  Lakes  Dora  and  Yale  were  collected  at  every  fourth  station  during  phase  2  and  at  every 
third  station  in  Lake  Wauburg.  To  determine  chlorophyll  concentration,  250-mL  water 
subsamples  were  vacuum  filtered  through  47-mm  diameter,  0.45-um  pore  nylon 
membranes  (Whatman  International  Ltd,  Maidstone,  UK,  No.  7404-004)  in  a  dark  room 
(with  only  a  single  green  light)  to  reduce  photodegradation  of  chlorophyll.  All 
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subsequent  analyses  of  chlorophyll  were  also  conducted  under  green  light.  Filters  with 
phytoplankton  residue  were  stored  in  covered  dishes  wrapped  in  aluminum  foil  to 
exclude  light  and  maintained  at  -10  °C  until  processed  within  30  days.  The  filters  were 
thawed,  cut  into  5-mm  strips  and  placed  into  capped  test  tubes  containing  20-mL 
(sufficient  to  completely  cover  the  filter  strips)  N,N-dimethylformamide.  The  tubes  were 
wrapped  in  aluminum  foil  and  shaken  for  24  hours.  The  extract  was  filtered  through  a 
0.45-um  pore  nylon  filter  and  absorbance  determined  with  a  spectrophotometer  at  647 
and  664.5  nm  (Inskeep  and  Bloom  1985,  Moran  and  Porath  1980).  Calculations  of 
chlorophyll  content  were  made  as  follows: 

chlorophyll  a  -  12.7ABS664  5  -  2.79ABSM7 

chlorophyll  b  =  20.7ABS647  -  4.62ABS664  5 

total  chlorophyll  =  17.6ABS647+  8.08ABS664.5 

For  all  benthic  samples,  one  1 5  x  1 5  x  1 5  cm  Ekman  dredge  sample  was  collected 
from  each  station  and  washed  through  a  350-um  mesh  sieve  bucket.  Washed  samples 
were  processed  and  chironomid  larvae  counted  in  a  30  x  40  cm  gridded  white  pan  using 
standard  methods  (Ali  et  al.  1977)  and  identified  using  the  keys  of  Epler  (1995).  Before 
washing,  approximately  1  cm3  of  sediment  from  each  benthic  sample  was  transferred  to 
labeled  polyethylene  scintillation  vials  for  sediment  DW  determination  as  previously 
described  (Chapter  4).  These  DW  data  were  also  used  to  augment  the  final  sediment 
maps  prepared  in  Chapter  4.  From  each  benthic  sample,  G.  paripes  larvae  were  separated 
and  length  of  each  measured  (tip  of  head  capsule  to  tip  of  terminal  segment)  to  the 
nearest  mm.  Each  1  -mm  size  class  was  composited  and  dried  at  60  °C  for  24  hours  to 
determine  dry  biomass  (Dermott  and  Paterson  1974);  overall  weight  of  each  size  class 
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was  divided  by  total  number  of  larvae  to  estimate  individual  biomass.  Larval  length  and 
dry  biomass  were  fitted  to  a  regression  function  to  estimate  a  length-biomass  relationship. 
The  result  of  this  analysis  was  used  to  estimate  productivity  using  the  size  frequency 
method  (Hynes  and  Coleman  1968,  Hamilton  1969)  with  the  generation  time 
modification  of  Benke  (1979)  using  the  formula: 

Productivity  =  [^n,-n,+1)(^)]^ 

Where  i  =  number  of  size  classes,  n,  =  mean  number  of  individuals  of  size  classy,  B7  = 
individual  biomass  for  size  classy,  and  CPI  =  cohort  production  interval  in  days.  For  this 
study,  CPI  was  determined  using  overall  mean  water  temperature  observed  in  the  lakes 
and  development  time  calculated  in  Chapter  7  for  that  temperature. 

All  chironomid  larvae  other  than  G.  paripes  were  counted  and  identified  to  lowest 
possible  taxonomic  level;  these  larvae  collected  from  each  station  were  composited  and 
dried  as  above  to  determine  biomass.  During  phase  1  of  the  study,  specific  records  were 
maintained  for  C.  crassicaudatus,  G.  paripes,  Tanytarsini,  Tanypodinae,  other 
Chironomidae  and  total  Chironomidae.  With  the  sampling  plan  re-evaluation  after  phase 
1,  additional  specific  records  were  also  maintained  for  Cryptochironomus  spp.  and 
Psuedochironomus  spp. 

Results  and  Discussion 
Monthly  mean  values  of  selected  water  physico-chemical  parameters  and 
sediment  DW  are  presented  in  Table  3  for  Lake  Dora.  Mean  water  depth  varied  between 
2.1  and  2.8  m,  with  generally  low  water  levels  occurring  during  the  last  four  months  of 
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the  study  because  of  a  prolonged  drought  in  central  Florida.  Water  clarity  in  this  lake 
was  relatively  poor  during  the  entire  study,  with  Secchi  disk  transparency  ranging  from 
19  to  31  cm.  A  range  of  12.8-28.7  °C  daytime  water  temperatures  was  noted,  with  water 
remaining  below  20  °C  between  November  and  March.  Concentration  of  dissolved 
oxygen  varied  between  5.3  and  9.2  ppm,  with  lowest  values  occurring  during  summer 
and  highest  values  during  winter.  Specific  conductance  values  in  the  lake  ranged  from 
338  to  473  uS/cm.  Primary  productivity  indicated  by  phytoplankton  concentration  in 
Lake  Dora  was  high;  chlorophyll  a,  chlorophyll  b,  and  total  chlorophyll  amounted  to 
163.3-316.5,  3.3-46.9,  and  177.7-329.0  ug/m3,  respectively.  Sediment  DW  was  in  the 
range  of  4 1.5  to  59.4%. 

Monthly  mean  values  of  relevant  water  and  sediment  D  W  data  from  Lake  Yale 
are  presented  in  Table  4.  Mean  water  depth  varied  between  2.7  and  3.2  m  during  the 
study  period.  No  samples  could  be  collected  from  Lake  Yale  during  January  2001 
because  of  inaccessibility  due  to  receded  water  levels  making  boat  ramps  unusable. 
Water  clarity  in  Lake  Yale  was  greater  than  in  Lake  Dora,  as  indicated  by  the  45-127  cm 
range  of  monthly  mean  Secchi  disk  transparency.  Daytime  water  temperatures  ranged 
between  12.9  and  29.3  °C,  dissolved  oxygen  4.6-9.8  ppm,  and  conductivity  293-425 
uS/cm.  Phytoplankton  concentration  in  Lake  Yale  was  much  lower  compared  to  Lake 
Dora;  chlorophyll  a  varied  from  13.0  to  94.2  ug/m3,  chlorophyll  b  from  0.2  to  42.3 
ug/m3,  and  total  chlorophyll  from  14.1  to  136.5  ug/m3.  Sediment  DW  ranged  between 
35.6  and  50.0%. 

Water  depth  in  Lake  Wauburg  ranged  between  2.1  and  2.8  m  (Table  5).  Secchi 
disk  transparency  values  ranged  between  59  and  85  cm.  Lake  Wauburg  exhibited  a  wider 
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range  of  daytime  mean  water  temperatures  (10.5-30.2  °C).  Dissolved  oxygen  values 
ranging  between  6.1  and  10.7  ppm  were  generally  higher  compared  to  Lakes  Dora  or 
Yale.  Specific  conductance  (75-1 18  uS/cm)  was  considerably  lower  than  Lakes  Dora 
and  Yale,  representing  a  lower  concentration  of  dissolved  solids.  Chlorophyll  a  varied 
from  16.9  to  73.2  ug/m3,  chlorophyll  b  0.9  to  59.2  ug/m3,  and  total  chlorophyll  from  17.9 
to  134.5  ug/m3.  Dry  weight  of  sampled  sediments  ranged  from  21.3  to  37.7%. 

Monthly  mean  densities  of  chironomid  larvae  from  Lake  Dora  are  summarized  in 
Table  6.  The  target  species,  G.  paripes,  comprised  only  a  small  proportion  of  the  total 
larval  samples,  forming  7.6%  of  total  chironomids  over  both  phases  of  the  study,  and 
only  0.9%  of  the  total  in  the  second  phase  of  the  study.  Initial  values  of  G  paripes 
density  were  rather  high  (2,197/m2  in  March  1999),  but  this  species  declined  quickly 
through  September  1999,  with  numbers  remaining  low  throughout  the  rest  of  the  study, 
except  for  modest  increases  in  density  during  May  to  July  2000.  This  rapid  decline  is 
similar  to  the  reported  decline  in  G  paripes  adult  capture  by  New  Jersey  light  traps 
around  Lake  Monroe  during  1981  and  1982  (Ali  and  Fowler  1985).  It  would  be  of 
interest  to  see  if  the  G  paripes  population  in  Lake  Dora  increases  again  in  the  next  1-2 
years  in  a  manner  similar  to  the  Lake  Monroe  population  in  1983  and  1984.  Tanytarsini 
(mostly  Cladotany tarsus  spp.)  larvae  were  numerically  the  most  abundant,  forming 
74.3%  of  total  midge  larvae  collected  in  this  lake.  Peak  larval  densities  occurred  during 
late  autumn  and  early  winter  (November-January).  The  pest  species  Chironomus 
crassicaudatus  constituted  6.5%  and  Tanypodinae  1.7%  of  total  larvae  collected  during 
the  study.  Cryptochironomus  spp.  amounted  to  9.4%  of  total  larvae  collected  during 
phase  2  (the  only  phase  they  were  specifically  recorded).  Peak  densities  of 
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C.  crassicaudatus  larvae  during  spring  and  early  summer  were  consistent  with  the  trends 
reported  for  this  species  in  Lake  Monroe  (Ali  et  al.  1996),  though  the  early  population 
increase  during  January  and  February  2001  was  unusual  and  may  have  been  a  result  of 
the  unusually  low  water  levels  in  the  lake  after  an  extended  drought.  Relatively  low 
numbers  of  Tanypodinae  were  collected  from  Lake  Dora,  with  monthly  mean  numbers 
ranging  between  3  and  116  larvae/m2,  and  forming  1.7%  of  total  collection.  Small 
numbers  of  other  chironomids  were  also  identified  in  Lake  Dora,  including  species  of 
Orthocladiinae,  Polypedilum  spp.,  Pseudochironomus  spp.,  and  Geoldichironomus  spp. 
Larval  G.  paripes  density  in  Lake  Yale  (Table  7)  was  also  high  during  March 

1999  (1,782/m2),  declining  to  a  minimum  of  3/m2  in  September  1999  and  remaining 
below  100/m2  until  June  2000.  Population  trends  of  G.  paripes  during  summer  and  fall  of 

2000  were  similar  to  those  previously  reported  for  this  species  in  Lake  Monroe  (Ali  et  al. 
1996).  Glyptotendipes  paripes  formed  a  greater  proportion  of  total  chironomids  (25.4% 
of  total  collection)  in  Lake  Yale  than  in  Lake  Dora.  Chironomus  crassicaudatus 
population  peaked  during  May  1 999,  December  1 999  and  March  2000,  and  this  trend  was 
similar  to  that  reported  by  Ali  et  al.  (1996)  in  Lake  Monroe;  however,  C.  crassicaudatus 
population  in  Lake  Yale  fell  below  detectable  levels  after  June  2000.  Chironomus 
crassicaudatus  larvae  were  absent  from  the  samples  during  the  last  8  months  of  the  study; 
this  species  comprised  only  2.6%  of  total  collection.  Numerically,  Tanytarsini  (mostly 
Cladotany tarsus  spp.)  were  the  most  common  midge  in  this  lake,  though  densities  and 
overall  proportion  (53.9%)  were  below  that  of  Lake  Dora.  Population  peak  was  observed 
between  February  and  May  2000,  and  a  general  increase  in  the  following  autumn  and 
winter,  with  a  notable  increase  in  February  2001.  Density  of  Cryptochironomus  spp. 


47 
larvae  in  Lake  Yale  remained  generally  between  30  and  60/m2  during  the  second  phase  of 
the  study,  with  a  maximum  of  78  larvae/m2  during  May  2000  and  a  minimum  of  6 
larvae/m2  in  September.  Density  of  Pseudochironomus  spp.  ranged  between  44  larvae/m2 
(September  2000)  and  228  larvae/m2  (May  2000).  Tanypodinae  were  slightly  more 
numerous  in  Lake  Yale  compared  to  Lake  Dora  (27-124  larvae/m2,  4.1%  of  total 
chironomids).  Orthocladiinae  and  Polypedilum  spp.  larvae  were  found  in  low  numbers  in 
Lake  Yale.  No  samples  were  collected  during  January  2001 . 

Glyptotendipes  paripes  larvae  clearly  dominated  the  chironomid  fauna  of  Lake 
Wauburg,  comprising  97.8%  of  collected  larvae  during  the  study  period  (Table  8).  Mean 
density  varied  from  5,206  larvae/m2  (August  2000)  to  19,442  larvae/m2  (March  2000). 
Any  obvious  seasonal  trend  was  lacking,  though  there  was  some  tendency  for  higher 
densities  to  occur  closer  to  cooler  months  (mean  density  exceeded  10,000  larvae/m2  in 
March,  April,  November  2000  and  January  2001).  Tanypodinae  in  Lake  Wauburg  were 
encountered  at  higher  densities  (monthly  mean  densities  26-361  larvae/m2)  compared  to 
Lakes  Dora  and  Yale,  though  their  percent  composition  of  total  chironomids  remained 
relatively  low  (1.1%).  Chironomus  crassicaudatus,  Tanytarsini,  Polypedilum  spp.  and 
species  of  Orthocladiinae  were  found  in  low  numbers  in  Lake  Wauburg  during  the  study. 
Density  distributions  of  larval  G.  paripes  were  plotted  using  the  same  techniques 
as  the  mapping  study  in  Chapter  4.  A  log(n+l)  transformation  was  used  to  enhance  data 
resolution  of  the  highly  variable  densities  recorded.  Phase  1  results  for  G.  paripes  larvae 
in  Lake  Dora  are  presented  in  Figures  15-17,  while  the  phase  2  results  are  shown  in 
Figures  18-23.  During  March  1999,  the  highest  densities  were  found  in  the  vicinity  of 
the  narrows  between  the  east  and  west  basins  of  the  lake,  and  peripherally  around  much 
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of  the  west  basin;  this  trend  was  repeated  in  May  1 999,  though  the  densities  were  lower. 
After  May  1999,  G  paripes  larvae  were  found  only  in  small,  low  density  clusters,  mostly 
in  the  western  basin  and  usually  associated  with  areas  of  sandy  substrates. 

Larval  G  paripes  distributions  in  Lake  Yale  are  shown  in  Figures  24-26  (phase  1) 
and  27-32  (phase  2).  From  March  to  December  1999,  G.  paripes  distributions  were 
similar  to  those  found  in  Lake  Dora,  primarily  around  the  lake  margin  and  associated 
with  sandy  substrates.  Starting  in  February  2000,  clusters  of  G.  paripes  larvae  were 
uncharacteristically  found  in  deeper  water  areas  with  soft  sediments.  High  densities 
(>  1,000  larvae/m2)  were  found  in  these  deep-water,  soft  sediment  areas,  especially  in  the 
center  of  the  large  northern  basin,  with  low  or  zero  numbers  occurring  in  the  peripheral 
areas.  Starting  in  November  2000,  while  still  predominant  in  the  lake  center,  higher 
densities  of  G.  paripes  were  found  along  the  sandy  nearshore  areas.  Only  the  northern 
portion  of  Lake  Yale  was  sampled  in  December  2000,  and  no  samples  were  collected 
during  January  2001. 

As  noted  in  Table  8,  Lake  Wauburg  had  dramatically  higher  densities  of 
G.  paripes  larvae,  as  is  also  evident  in  Figures  33-38.  Many  areas  of  this  lake  sustained 
very  high  densities  (>  10,000  larvae/m2)  throughout  the  sampling  period.  These  high- 
density  areas  were  found  around  the  lake  margin  in  association  with  the  firmer  substrates, 
with  the  lowest  densities  found  in  the  lake  center  and  along  the  northwest  shoreline, 
where  the  soft  sediments  existed  up  to  the  shoreline.  Due  to  crowding  from 
overabundance,  the  population  of  this  species  extended  to  the  soft  sediments  close  to  the 
sandy  substrates. 
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The  preference  of  G.  paripes  larvae  in  Lakes  Dora  and  Wauburg  for  sand  and 
transitional  sand  to  muck  substrates  is  consistent  with  previously  published  information 
on  this  species.  Ali  and  Baggs  (1982)  reported  a  trend  for  decreased  larval  populations  in 
Lake  Monroe  and  in  a  nearby  water  cooling  reservoir  with  increased  organic  carbon 
content  of  sediments.  This  was  also  reported  by  Ali  et  al.  (1998b)  for  G.  paripes  larvae 
in  Lake  Jesup  (Seminole  County,  Florida).  Provost  (1957)  reported  a  strong  preference 
for  sand  or  sand/peat  bottoms  by  G.  paripes  larvae  in  1 3  lakes  in  the  Winter  Haven, 
Florida  area.  Lakes  in  Polk  County,  Florida  with  large  areas  of  sandy  bottom,  reduced 
emergent  vegetation  and  more  open  shoreline  generally  produced  larger  numbers  of  adult 
G.  paripes  than  those  containing  more  muck  at  the  bottom  (Callahan  and  Morris  1987). 
Curry  (1962)  in  a  review  classified  G.  paripes  larvae  as  preferring  marl,  plants,  sand  and 
peat  bottoms.  While  preferring  sand  bottom,  G.  paripes  larvae  were  noted  in  deeper 
waters  of  Lake  Thonotossassa  (Hillsborough  County,  Florida)  by  Cowell  and  Vodopich 
(1981).  Milleson  (1978)  reported  that  G.  paripes  larvae  in  the  Istokpoga  lake  chain  were 
predominant  in  sand  or  sand/mud  substrates.  While  the  occurrence  of  G  paripes  larvae 
in  deeper,  soft  sediments  has  not  been  previously  reported  for  Florida  lakes,  a  mixed 
population  of  G  paripes  and  Glyptotendipes  glaucus  in  Lough  Neagh,  Ireland,  was 
reported  to  prefer  mud  sediments  with  a  water  depth  of  4-m  (Carter  1975).  Also, 
G.  paripes  was  reported  to  exist  at  densities  over  1,000  larvae/m2  in  a  shallow  bog  lake 
(Baxter  Lough)  in  England  with  a  peat/mud  bottom  (McLachlan  1976). 

The  influence  of  water  and  sediment  physico-chemical  parameters  over  time  on 
G.  paripes  larval  populations  was  investigated  using  multiple  linear  regression  on 
monthly  mean  values  transformed  using  log(n+l).  Contradictory  results  were  found  for 
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Lakes  Dora  and  Yale:  Chlorophyll  a  concentration  showed  a  positive  weak  correlation 
with  G.  paripes  larval  density  in  Lake  Dora  [G.  paripes  density  =  2.043  +  0.014(Chl  a), 
r2  =  0.378,  P  =  0.009],  while  this  relationship  was  negative  in  Lake  Yale  [G.  paripes 
density  =  3.02  -  0.018(Chl  a),  r2  -  0.292,  P  =  0.031].  No  other  significant  correlations  of 
monthly  mean  G.  paripes  larval  densities  with  monthly  mean  water  physico-chemical 
parameters  were  detected  in  these  two  lakes.  Glyptotendipes  paripes  larval  density  was 
positively  correlated  with  water  depth  in  Lake  Wauburg  (larval  density  =  2.35  + 
0.651  depth,  r2  =  0.593,  P  =  0.15).  Performing  multiple  regression  on  monthly  mean 
water  and  sediment  physico-chemical  parameters  and  G.  paripes  larval  density  from  all 
three  study  lakes  yielded  the  following  equation:  G  paripes  larval  density  =  2.39  + 
0.016(Secchi  disk  transparency)  +  0.047(water  temperature)  -  0.006(specific 
conductance),  indicating  highly  significant  relationship,  (r2  =  0.726  and  P  <  0.0001). 

To  elucidate  the  position  of  G  paripes  in  the  chironomid  community  of  the 
investigated  lakes,  canonical  correspondence  analysis  (CCA)  (Braak  and  Smilauer  1998) 
was  conducted  on  the  data  using  the  software  Canoco  for  Windows®  version  4.0  (Centre 
for  Biometry,  CPRO-DLO,  Wageningen,  Netherlands  and  Microcomputer  Power,  Ithaca, 
NY).  The  utility  of  this  analysis  to  differentiate  habitat  clusters  within  the  chironomid 
community  of  a  lake  was  established  by  Stimac  and  Leong  (1977).  These  authors  were 
able  to  indicate  that  while  there  were  three  distinct  habitats  within  the  aquatic  weeds  of 
Laguna  Lake,  CA,  these  were  not  a  major  limiting  factor  in  chironomid  abundance  and 
that  only  almost  complete  elimination  of  the  weeds  would  have  a  significant  effect  on 
nuisance  midge  populations.  In  the  first  step,  the  environmental  variables  that 
significantly  affect  community  composition  were  selected  by  forward  selection  (P  <  0.05, 
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Monte  Carlo  permutation  test  with  199  permutations).  The  ordination  axes  in  CCA  were 
then  constrained  to  these  environmental  parameters.  Results  of  this  analysis  for  Lake 
Dora  are  shown  in  Figure  39.  Type  of  sediment,  namely  presence  of  sand,  muck  and  DW 
comprised  one  gradient  in  the  community  data  that  corresponded  to  the  first  ordination 
axis.  Water  temperature,  detritus  presence,  chlorophyll  b  and  water  depth  composed 
another  gradient  in  community  data  that  corresponded  with  the  second  ordination  axis.  It 
is  obvious  that  C.  crassicaudatus  increased  with  water  depth  and  muck  presence,  whereas 
Cladotanytarsus  spp.  and  other  Tanytarsini  responded  to  this  trend  negatively  (Figure 
39).  Glyptotendipes  paripes  responded  positively  to  water  temperature,  while 
Polypedilum  spp.  and  Geoldichironomus  spp.  responded  negatively  to  this  parameter. 
Sample  plots  clustered  into  roughly  two  groups,  shallower  sand  bottoms  and  deeper  muck 
bottoms. 

Presence  of  muck  was  the  strongest  predictor  of  G  paripes  density  in  Lake  Yale 
using  CCA  analysis  (Figure  40),  with  water  depth  also  being  important. 
Cryptochironomus  spp.  and  Cladotanytarsus  spp.  density  increased  with  sediment  DW. 
Polypedilum  spp.,  Pseudochironomus  spp.  and  other  Tanytarsini  species  were  influenced 
by  presence  of  rooted  vegetation.  The  chironomid  community  in  this  lake  formed  three 
separate  habitat  clusters:  in  shallower  water  with  sand  bottom,  in  sand  bottoms  with 
rooted  vegetation,  and  in  deep  water  muck  bottoms. 

Tanypodinae  increased  with  muck  presence  and  water  depth  in  Lake  Wauburg 
(Figure  41),  while  Tanytarsini  increased  with  Secchi  disk  transparency,  chlorophyll  b, 
and  total  chlorophyll.  Glyptotendipes  paripes  showed  minimal  influence  by  the 
examined  parameters.  Sample  plots  generally  exhibited  a  similar  clustering  (shallower 
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water/sand  bottom  and  deeper  water/muck  bottom)  to  the  habitat  clustering  in  Lake  Dora 
(Figures  39  and  41). 

The  finding  that  sediment  characteristics  and  water  depth  were  common 
influences  on  chironomid  distributions  in  these  lakes  is  in  agreement  with  the  findings  of 
Verneaux  and  Aleya  (1998),  who  reported  that  chironomid  communities  in  Lake  Abbaye 
(France)  were  primarily  influenced  by  bathymetric  and  sediment  gradients  within  the  lake 
and  seasonal  variations  in  water  characteristics.  The  authors  also  reported  that  seasonal 
variation  increased  with  water  depth,  with  deeper  areas  having  a  more  pronounced 
seasonal  change  in  chironomid  fauna. 

The  size  of  field-collected  G  paripes  larvae  ranged  between  3  and  16  mm  in 
overall  length.  Based  on  the  larval  sizes  recorded  during  the  development  study  in 
Chapter  7,  the  processed  samples  retained  mostly  3  rd  and  4th  instar  larvae,  indicating  that 
actual  larval  densities  were  potentially  much  higher.  The  lack  of  data  on  1st  and  2nd  instar 
population  information  also  complicated  productivity  estimates  of  G.  paripes  in  these 
lakes;  however,  since  3 -mm  long  larvae  had  relatively  small  biomass  (Table  9),  the 
overall  loss  was  probably  not  as  great  to  the  estimate.  Mean  dry  biomass  of  G.  paripes 
larvae  by  1-mm  size  increment  is  shown  in  Table  9.  For  estimating  a  larval 
length/biomass  relationship,  the  dry  biomass  values  were  transformed  by  the  cube  root  to 
correct  for  the  linear  length  to  cubic  biomass  relationship.  This  was  analyzed  by  linear 
regression  and  the  data  fit  was  exceptional  [cube  root  mg  dry  biomass  =  -0.088  + 
0.101(mm  length);  r2  -  0.992,  P  <  0.0001).  Productivity  of  3rd  and  4th  instar  G  paripes 
larvae  between  March  2000  and  February  2001  was  estimated  using  the  size-frequency 
method  (Hynes  and  Coleman  1968,  Hamilton  1969)  with  the  modification  of  Benke 
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(1979).  For  this  modification,  CPI  was  estimated  at  51  days  using  an  approximate  mean 
water  temperature  of  23  °C  and  the  temperature  development  times  estimated  in  Chapter 
7.  Lake  Dora  had  the  lowest  G.  paripes  annual  productivity  of  the  lakes  in  this  study, 
205.1  mg/m2  dry  biomass;  followed  by  Lake  Yale  3,482  mg/m2  dry  biomass,  and  Lake 
Wauburg  1 56,900  mg/m2  dry  biomass.  This  exceptionally  high  productivity  value  for 
G.  paripes  larvae  in  Lake  Wauburg  is  close  to  the  high  productivity  of  162  g/m2  dry 
weight  for  Glyptotendipes  barbipes  in  a  sewage  lagoon  that  was  reported  by  Benke 
(1984)  and  Kimerle  and  Anderson  (1971).  Mean  standing  crop  during  the  same  time 
interval  was  6.8  mg/m2  in  Lake  Dora,  189.7  mg/m2  in  Lake  Yale,  and  6,081.5  mg/m2  in 
Lake  Wauburg.  These  productivity  estimates  indicate  a  general  trend  of  increased 
G.  paripes  populations  with  increased  trophic  level  (such  as  between  Lake  Yale  and  Lake 
Wauburg),  then  decreasing  with  hypereutrophic  conditions  (difference  between  Lake 
Wauburg  and  Lake  Dora).  Biomass  turnover  in  these  lakes  was  measured  by  calculating 
annual  productivity/biomass  (P/B)  ratios  of  30.2  for  Lake  Dora,  18.4  for  Lake  Yale,  and 
25.8  for  Lake  Wauburg.  These  values  are  consistent  with  those  reviewed  by  Benke 
(1984)  for  lakes  with  high  chironomid  productivity. 
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Table  3.  Monthly  mean  ±  SD  of  selected  water  physico-chemical  parameters  (water 
depth,  Secchi  disk  transparency,  water  temperature,  dissolved  oxygen,  specific 
conductance,  chlorophyll  a,  chlorophyll  b  and  total  chlorophyll)  and  sediment  dry 
weight  (DW)  for  Lake  Dora  (Lake  County,  Florida)  from  March  1999  to  February 
2001. 


Date 

Depth 

Secchi 

Temp. 

DO 

Cond. 

DW 

Chi  a 

Chi  b 

Total 

n 

(m) 

Transp. 
(cm) 

(°C) 

(ppm) 

(uS/cm) 

(%) 

(ug/m3) 

(ug/m3) 

Chi 

(ug/m3) 

Mar  99 

2.6  ± 

31  ±3 

17.9  ± 

9.2  ± 

356  ± 

46.7  ± 

250.6  ± 

11.7* 

262.2  ± 

120 

0.8 

0.8 

1.3 

3 

34.7 

18.6 

7.4 

22.1 

May  99 

2.4  ± 

23  ±3 

24.8  ± 

5.9  ± 

369  ± 

53.9  ± 

316.5  ± 

12.6* 

329.0  ± 

120 

0.9 

1.4 

2.0 

22 

31.6 

26.1 

1.5 

27.0 

Jul  99 

2.2  ± 

26  ±2 

28.7  ± 

5.3  ± 

421  ± 

57.3  ± 

210.6± 

9.2  ± 

219.8  ± 

60 

0.8 

1.5 

2.1 

11 

32.4 

31.3 

3.4 

33.2 

Sep  99 

2.6  * 

27  ±2 

25.5  ± 

5.8  ± 

396  ± 

41.5  ± 

214.9  ± 

4.8* 

219.7* 

60 

0.8 

0.3 

0.7 

5 

36.3 

32.7 

3.3 

35.3 

Nov  99 

2.6  ± 

28  ±3 

18.6  ± 

7.9  ± 

351  * 

59.4  ± 

239.0  ± 

11.6* 

250.6* 

90 

0.9 

1.6 

0.5 

11 

30.5 

41.8 

2.7 

44.3 

Jan  00 

2.7  ± 

26  ±2 

16.4  ± 

9.0  ± 

338  ± 

55.2  ± 

204.2  ± 

13.2  ± 

217.4* 

120 

0.9 

1.1 

7.7 

8 

32.2 

41.8 

6.3 

46.4 

Mar  00 

2.7  * 

25  ±2 

20.4  ± 

7.1  ± 

369  ± 

54.7  ± 

256.0  ± 

46.9  ± 

302.8* 

80 

0.9 

0.7 

1.4 

5 

31.7 

15.4 

1.3 

16.3 

Apr  00 

2.4  ± 

23  ±3 

22.0  * 

7.4  ± 

391  ± 

50.3  ± 

259.1  ± 

10.6* 

269.6  * 

80 

0.9 

1.6 

1.5 

10 

34.5 

33.6 

1.8 

35.1 

May  00 

2.5  ± 

23  ±3 

23.8  ± 

6.7  ± 

415  ± 

53.2  ± 

234.3  ± 

15.2  ± 

249.4  * 

80 

0.8 

0.7 

1.2 

6 

33.8 

25.9 

2.7 

31.0 

JunOO 

2.1  ± 

21  ±2 

27.8  ± 

- 

473  ± 

57.5  ± 

216.9* 

35.7  ± 

252.5  * 

80 

0.8 

1.0 

8 

31.8 

30.5 

9.9 

39.2 

Jul  00 

2.3  * 

19*2 

28.1  ± 

5.7  ± 

455  ± 

53.3  ± 

239.4  ± 

8.2  ± 

247.5  * 

80 

0.9 

0.3 

1.5 

13 

33.2 

28.1 

2.0 

29.8 

Aug  00 

2.3  ± 

19±2 

28.1  ± 

5.4  ± 

446  ± 

57.5  ± 

218.3  ± 

3.3  ± 

221.5* 

40 

0.6 

0.4 

1.1 

5 

32.7 

32.7 

0.9 

20.5 

Sep  00 

2.3  ± 

21  ±2 

27.9  ± 

5.4  ± 

450  ± 

49.9  * 

206.1  * 

8.6  ± 

214.7* 

80 

0.8 

0.4 

0.9 

6 

35 

30.5 

1.3 

31.6 

Oct  00 

2.3  ± 

21  ±4 

22.3  ± 

7.4  ± 

410* 

54.8  ± 

222.7  * 

12.7  ± 

235.3  * 

80 

0.9 

2.6 

0.7 

20 

33.4 

14.3 

0.9 

14.6 

Nov  00 

2.2  * 

22  ±3 

21.9± 

6.4  ± 

428  ± 

56.8  ± 

- 

. 

_ 

80 

0.7 

0.5 

1.1 

3 

31.9 

Dec  00 

2.1  ± 

22  ±3 

17.3  ± 

8.6  ± 

395  ± 

56.6  ± 

218.8± 

5.0* 

223.8* 

80 

0.8 

1.1 

1.6 

6 

32.2 

12.1 

1.0 

12.5 

Jan  01 

2.2  ± 

28  ±3 

12.8  * 

9.1  ± 

404  ± 

47.6  ± 

163.3  ± 

45.5  ± 

208.7* 

80 

0.9 

1.1 

0.9 

23 

34.9 

19.5 

2.8 

21.8 

Feb  01 

2.1  * 

25  ±2 

14.4  ± 

8.5  ± 

- 

51.1* 

170.6  * 

7.1* 

177.7* 

80 

0.9 

0.5 

1.4 

17.9 

17.9 

1.8 

19.4 

55 


Table  4.  Monthly  mean  ±  SD  of  selected  water  physico-chemical  parameters  (water 
depth,  Secchi  disk  transparency,  water  temperature,  dissolved  oxygen,  specific 
conductance,  chlorophyll  a,  chlorophyll  b  and  total  chlorophyll)  and  sediment  dry 
weight  (DW)  for  Lake  Yale  (Lake  County,  Florida)  from  March  1 999  to  February 
2001. 


Date 

Depth 

Secchi 

Temp. 

DO 

Cond. 

DW 

Chi  a 

Chi  b 

Total 

n 

(m) 

Transp. 
(cm) 

(°C) 

(ppm) 

(US/cm) 

(%) 

(ug/m3) 

(Hg/m3) 

Chi 

(ug/m3) 

Mar  99 

3.0  ± 

87  ±  10 

23.5  ± 

7.8  ± 

326  ± 

44.5  ± 

29.2  ± 

3.2  ± 

32.3  ± 

120 

1.2 

2.1 

0.8 

3 

31.6 

7.3 

5.2 

10.9 

May  99 

2.9  ± 

62  ±7 

27.1  ± 

6.3  ± 

388  ± 

39.6  ± 

45.8  ± 

4.4  ± 

50.2  ± 

80 

1.5 

1.1 

1.5 

11 

29.2 

11.4 

1.5 

11.7 

Jul  99 

2.7  ± 

65  ±5 

29.3± 

5.7  ± 

408  ± 

42.1  ± 

57.9  ± 

15.5  ± 

73.4  ± 

60 

1.2 

0.6 

1.3 

7 

31.6 

21.5 

30.8 

51.1 

Sep  99 

2.9  ± 

65  ±5 

26.5  ± 

5.6  ± 

368  ± 

47.7  ± 

53.6  ± 

0.2  ± 

53.8  ± 

30 

1.3 

0.8 

1.7 

5 

28.9 

14.3 

0.3 

14.6 

Dec  99 

2.9  ± 

57  ±5 

18.2  ± 

7.5  ± 

313  ± 

50.0± 

59.7  ± 

4.1  ± 

63.8  ± 

30 

1.2 

0.6 

1.3 

3 

31.2 

7.9 

0.7 

7.3 

Feb  00 

3.2  ± 

101  ±27 

14.4  ± 

8.8  ± 

293  ± 

44.3  ± 

22.0  ± 

6.6  ± 

28.6  ± 

120 

1.4 

1.5 

6.5 

12 

31.3 

7.6 

2.0 

9.0 

Mar  00 

3.0  ± 

50  ±4 

20.3  ± 

9.1  ± 

336  ± 

45.0  ± 

94.2  ± 

42.3  ± 

136.5  ± 

80 

1.6 

0.3 

0.8 

5 

30.7 

6.7 

1.0 

7.0 

Apr  00 

2.7  ± 

46  ±3 

20.5  ± 

7.8  ± 

340  ± 

40.7  ± 

71.7  ± 

5.9  ± 

77.6  ± 

60 

1.3 

0.9 

0.6 

5 

30.3 

11.7 

2.3 

12.6 

May  00 

2.7  ± 

45  ±3 

24.9  ± 

7.3  ± 

388  ± 

43.6  ± 

66.2  ± 

8.5  ± 

74.7  ± 

80 

1.1 

1.1 

2.6 

10 

31.4 

4.6 

0.7 

5.0 

JunOO 

2.8  ± 

49  ±3 

28.2  ± 

6.3  ± 

425  ± 

42.9  ± 

84.3  ± 

41.5  ± 

125.7  ± 

80 

1.3 

0.6 

1.1 

5 

31.9 

7.6 

1.8 

8.0 

Jul  00 

3.2  ± 

58±3 

28.0  ± 

4.6  ± 

425  ± 

44.1  ± 

47.6  ± 

1.1± 

48.7  ± 

80 

1.4 

0.7 

1.7 

12 

33.1 

7.2 

0.6 

7.5 

Aug  00 

2.8  ± 

61  ±5 

28.6  ± 

6.1  ± 

423  ± 

43.1  ± 

42.7  ± 

1.1  ± 

43.8  ± 

40 

1.3 

0.7 

1.6 

6 

32.5 

5.9 

0.6 

6.1 

Sep  00 

3.0  ± 

66  ±3 

27.0  ± 

7.0  ± 

404  ± 

42.5  ± 

50.3  ± 

3.3  ± 

53.6  ± 

80 

1.4 

0.4 

0.7 

5 

32.4 

5.6 

0.6 

5.7 

Oct  00 

2.7  ± 

67  ±3 

20.7  ± 

8.0  ± 

365  ± 

38.8  ± 

42.0  ± 

7.0  ± 

49.0  ± 

80 

1.3 

0.5 

0.7 

4 

31.5 

6.5 

1.0 

6.7 

Nov  00 

2.7  ± 

70  ±7 

16.7  ± 

9.3  ± 

340  ± 

39.5  ± 

- 

. 

. 

80 

1.1 

0.5 

0.5 

7 

32.6 

Dec  00 

2.9  ± 

77  ±4 

12.9  ± 

9.8  ± 

317  ± 

35.6  ± 

19.1  ± 

1.4± 

20.6  ± 

24 

0.9 

0.1 

0.9 

5 

34.4 

2.8 

1.0 

2.8 

Jan  01 
0 
Feb  01 

- 

- 

- 

- 

- 

- 

- 

- 

- 

2.7  ± 

127  ±5 

18.0  ± 

8.6  ± 

_ 

42.1  ± 

13.0± 

1.1  ± 

14.1  ± 

60 

1.2 

0.6 

0.7 

32.9 

2.4 

0.8 

2.9 

56 


Table  5.  Monthly  mean  ±  SD  of  selected  water  physico-chemical  parameters  (water 
depth,  Secchi  disk  transparency,  water  temperature,  dissolved  oxygen,  specific 
conductance,  chlorophyll  a,  chlorophyll  b  and  total  chlorophyll)  and  sediment  dry 
weight  (DW)  for  Lake  Wauburg  (Alachua  County,  Florida)  from  March  2000  to 
February  2001. 


Date 

Depth 

Secchi 

Temp. 

DO 

Cond. 

DW 

Chi  a 

Chi  b 

Total 

n 

(m) 

Transp. 
(cm) 

(°C) 

(ppm) 

(uS/cm) 

(%) 

(ug/m3) 

(ug/m3) 

Chi 

(ug/m3) 

Mar  00 

2.8  * 

79  ±6 

21.3  ± 

8.6  ± 

95  ±17 

31.0± 

73.2  * 

47.6  * 

120.8* 

30 

1.1 

0.5 

2.1 

26.7 

5.4 

1.7 

6.1 

Apr  00 

2.7  ± 

85  ±3 

22.7  ± 

7.0  ± 

91  ±  1 

32.9  ± 

38.5  ± 

11.7* 

50.2* 

30 

1.2 

0.6 

1.2 

27.2 

5.6 

9.0 

14.5 

May  00 

- 

67  ±4 

27.4  ± 

6.9  ± 

109  ±1 

21.9* 

58.7  ± 

15.8* 

74.5* 

30 

0.7 

1.0 

26.0 

4.5 

2.5 

5.9 

JunOO 

- 

72  ±3 

30.2  ± 

6.2  ± 

118*5 

35.4± 

73.5  ± 

49.6  * 

118.1* 

21 

1.2 

2.7 

29.6 

6.2 

1.1 

7.0 

Jul  00 

- 

63  ±  1 

29.2  ± 

6.3  ± 

115*5 

37.2  ± 

49.4  ± 

4.2  ± 

53.6* 

30 

1.6 

1.6 

30.4 

11.3 

1.2 

12.3 

Aug  00 

2.3  ± 

62  ±3 

29.5  ± 

7.0  ± 

118±5 

21.3  ± 

16.9  ± 

0.9  ± 

17.9* 

18 

0.7 

0.8 

1.8 

12.4 

13.8 

0.5 

13.5 

Sep  00 

2.5  ± 

63  ±4 

25.2  ± 

6.1  ± 

103  ±2 

33.6  ± 

71.1  * 

11.2* 

82.2* 

30 

0.9 

0.2 

1.2 

29.4 

6.5 

1.1 

7.4 

Oct  00 

2.5  * 

62  ±3 

22.5* 

7.8  ± 

99  ±3 

30.2  ± 

52.4  ± 

10.9* 

63.2* 

30 

0.9 

0.8 

1.5 

25.9 

3.8 

1.4 

4.6 

Nov  00 

2.6  ± 

59  ±2 

20.9* 

8.4  ± 

98  ±1 

31.5* 

- 

- 

- 

30 

0.8 

0.2 

0.7 

28.3 

Dec  00 

2.1  ± 

61  ±  3 

13.7  ± 

10.0  ± 

81  ±2 

37.7  ± 

50.4  ± 

6.2* 

56.6* 

30 

1.1 

0.7 

0.5 

28.7 

5.4 

1.4 

6.6 

Jan  01 

2.5  ± 

80  ±3 

10.5  ± 

10.7  ± 

75  ±2 

27.8  ± 

75.4  ± 

59.2* 

134.5* 

30 

0.9 

0.5 

0.7 

26.5 

24.3 

37.1 

60.7 

Feb  01 

2.3  ± 

68  ±4 

19.0  ± 

7.4  ± 

- 

34.5  ± 

72.0  ± 

17.1  * 

89.1  * 

30 

0.9 

0.6 

1.6 

28.8 

7.5 

10.5 

17.3 

57 


Table  6.  Mean  monthly  density  (No./m2)  ±  stratified  SD  of  selected  Chironomidae  larvae 
[Chironomus  crassicaudatus,  Glyptotendipes  paripes,  Cryptochironomus  spp. 
{Crypt.),  Tanytarsini  (Tanyt.,  mostly  Cladotany tarsus  spp.),  Tanypodinae 
(Tanyp.),  total  Chironomidae  larvae,  G.  p.  biomass  (mg/m2)  and  total 
Chironomidae  biomass  (mg/m2)]  collected  for  Lake  Dora  (Lake  County,  Florida), 
March  1999-February  2001. 


Date 

C.  crass. 

G.  paripes 

Crypt. 

Tanyt. 

Tanyp. 

Total 

a  p. 

biomass 

Total 
biomass 

Mar. 

2± 

2,197  ± 

- 

86  ± 

26  ± 

2,395 

98.9  ± 

146.2  ± 

99 

2 

172 

- 

13 

9 

±181 

8.6 

12.9 

May 

770  ± 

477  ± 

- 

1,957 

22  ± 

3,535 

34.4  ± 

670.8  * 

99 

129 

155 

±159 

4 

±241 

4.3 

64.5 

Jul. 

86  ± 

39  ± 

- 

602  ± 

9± 

907  ± 

25.8  ± 

172.0  ± 

99 

22 

9 

108 

4 

129 

4.3 

17.2 

Sep. 

13* 

9± 

- 

2,499 

3± 

2,589 

2.2  ± 

313.9± 

99 

9 

4 

±748 

2 

±744 

1.3 

73.1 

Nov. 

254  ± 

3± 

- 

4,123 

22  ± 

4,945 

1.7* 

1,333.9 

99 

52 

4 

±206 

9 

±219 

3.0 

±68.8 

Jan. 

73  ± 

13± 

- 

2,864 

52  ± 

3,307 

8.6  ± 

602.0  ± 

00 

22 

4 

±176 

9 

±181 

4.3 

38.7 

Mar. 

168  ± 

9* 

69  ± 

1,686 

22  ± 

1,965 

8.6  ± 

245.1  ± 

00 

34 

4 

4 

±  112 

9 

±112 

4.3 

25.8 

Apr. 

237  * 

9± 

138  ± 

1,337 

55  ± 

1,750 

3.4  ± 

172.0  ± 

00 

52 

4 

26 

±168 

4 

±181 

2.0 

21.5 

May 

275  ± 

39  ± 

215  ± 

1,969 

13± 

2,533 

12.9  ± 

210.7  ± 

00 

69 

13 

34 

±198 

4 

±237 

4.3 

34.4 

Jun 

82  ± 

52  ± 

129  ± 

882  ± 

17± 

1,264 

21.5* 

137.6  ± 

00 

17 

17 

9 

82 

4 

±90 

8.6 

8.6 

Jul. 

43  ± 

22  ± 

69  ± 

254  ± 

17± 

417± 

12.9  ± 

64.5  ± 

00 

13 

4 

4 

22 

4 

30 

4.3 

8.6 

Aug. 

17± 

9* 

77  ± 

400  ± 

9± 

520  ± 

1.7* 

43.0  ± 

00 

13 

4 

13 

60 

9 

69 

0.9 

8.6 

Sep. 

17± 

4± 

77  ± 

297  ± 

17± 

426  ± 

1.7* 

38.7  ± 

00 

9 

2 

9 

43 

4 

52 

0.9 

8.6 

Oct. 

9± 

4± 

112  ± 

774  ± 

34  ± 

946  ± 

3.9  ± 

51.6* 

00 

4 

3 

13 

86 

9 

95 

6.5 

8.6 

Nov. 

43  ± 

30  ± 

103  ± 

1,281 

56  ± 

1,531 

8.6  ± 

94.6  ± 

00 

17 

13 

9 

±  138 

9 

±142 

4.3 

8.6 

Dec. 

86  ± 

9± 

559  ± 

4,330 

116± 

5,207 

1.3* 

270.9  ± 

00 

34 

4 

30 

±245 

13 

±262 

0.4 

21.5 

Jan. 

219± 

4± 

194  ± 

1,974 

86  ± 

2,541 

4.3  ± 

399.9  ± 

01 

39 

4 

17 

±133 

26 

±142 

4.3 

43.0 

Feb. 

120  ± 

1± 

215  ± 

1,264 

77  ± 

1,707 

0.4  ± 

301.0± 

01 

22 

1 

17 

±90 

17 

±95 

0.4 

21.5 

58 


Table  7.  Mean  monthly  density  (No./m2)  ±  stratified  SD  of  selected  Chironomidae  larvae 
[Chironomus  crassicaudatus,  Glyptotendipes  paripes,  Cryptochironomus  spp. 
{Crypt.),  Pseudochironomus  spp.  {Pseud.),  Tanytarsini  (Tanyt.,  mostly 
Cladotany  tarsus  spp.),  Tanypodinae  (Tanyp.),  total  Chironomidae  larvae,  G.  p. 
biomass  (mg/m2)  and  total  Chironomidae  biomass  (mg/m2)]  collected  for  Lake 
Yale  (Lake  County,  Florida),  March  1 999-February  2001. 


Date 

C. 
crass. 

G. 
paripes 

Crypt. 

Pseud. 

Tanyt. 

Tanyp. 

Total 

G.p. 
biomass 

Total 
biomass 

Mar. 

1± 

1,782  ± 

- 

- 

193  ± 

45  ± 

2,038 

52.0  ± 

140.2  ± 

99 

1 

163 

7 

6 

±  168 

6.2 

10.8 

May 

173  ± 

182  ± 

- 

- 

243  ± 

34  ± 

1,141 

156.5  ± 

482.9  ± 

99 

29 

34 

44 

7 

±  143 

28.4 

56.8 

Jul. 

10± 

123  ± 

- 

- 

1104 

47  ± 

452  ± 

108.8  ± 

246.0  ± 

99 

5 

34 

14 

12 

67 

34.0 

64.5 

Sep. 

1  ± 

3± 

- 

- 

30  ± 

43  ± 

188  ± 

9.5  ±2.2 

83.0  ± 

99 

3 

1 

3 

21 

26 

14.2 

Dec. 

168  ± 

94 

- 

- 

4184 

60  ± 

890  ± 

3.4  ±  1.3 

500.9  ± 

99 

223 

7 

48 

23 

233 

373.7 

Feb. 

31± 

75  ± 

- 

- 

1,535  ± 

124  ± 

2,202 

64.0  ± 

595.6  ± 

00 

11 

15 

116 

16 

±123 

14.2 

41.7 

Mar. 

168  ± 

10± 

67  ±6 

- 

1,684  ± 

22  ± 

1,964 

9.5  ±4.3 

245.5  ± 

00 

34 

4 

111 

8 

±114 

25.4 

Apr. 

8± 

52  ± 

47  ±8 

161  ± 

1,008  ± 

60  ± 

1,350 

52.9  ± 

194.4  ± 

00 

4 

16 

24 

155 

10 

±162 

21.8 

33.2 

May 

3± 

814 

78  ±9 

228  ± 

2,327  ± 

97  ± 

2,779 

20.4  ± 

182.3  ± 

00 

4 

52 

28 

309 

11 

±314 

8.0 

20.6 

Jun. 

24 

1714 

37  ±4 

714  11 

550  ± 

27  ± 

864  ± 

47.3  ± 

94.6  ± 

00 

1 

58 

67 

7 

101 

18.6 

24.6 

Jul. 

0± 

377  ± 

48  ±13 

1304 

1534 

43  ± 

767  ± 

318.2± 

399.9  ± 

00 

0 

130 

27 

16 

13 

141 

120.0 

124.6 

Aug. 

0± 

342  ± 

28  ±3 

80  4  14 

1704 

43  ± 

669  ± 

313.9± 

356.9  ± 

00 

0 

284 

77 

14 

281 

264.5 

265.3 

Sep. 

0± 

461  ± 

64  1 

44  ±7 

64  ± 

38  ± 

625  ± 

172.0  ± 

210.7± 

00 

0 

272 

11 

14 

271 

105.8 

109.7 

Oct. 

0± 

648  ± 

1244 

61  ±20 

225  ± 

61  ± 

1,026 

653.6  ± 

722.4  4 

00 

0 

237 

29 

11 

±241 

247.7 

250.7 

Nov. 

0± 

235  ± 

39  ±6 

56  ±7 

1,138  ± 

514 

1,528 

133.3  ± 

202.1  ± 

00 

0 

61 

103 

12 

±  129 

41.7 

44.3 

Dec. 

0± 

691  ± 

484  1] 

72  ±21 

885  ± 

65  ± 

1,765 

167.7  ± 

270.9  ± 

00 

0 

284 

208 

18 

±342 

92.0 

86.4 

Jan. 

01 

Feb. 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0± 

359  ± 

50  ±9 

1774 

1,1714 

54  ± 

1,836 

197.8  ± 

348.3  ± 

01 

0 

126 

24 

123 

14 

±183 

116.9 

123.0 
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Table  8.  Mean  monthly  density  (No./m2)  ±  stratified  SD  of  selected  Chironomidae  larvae 
[Glyptotendipes  paripes,  Tanypodinae  (Tanyp.),  total  Chironomidae  larvae,  G.  p. 
biomass  (mg/m2)  and  total  Chironomidae  biomass  (mg/m2)]  collected  for  Lake 
Wauburg  (Alachua  County,  Florida),  March  2000-February  2001. 


Date 

G.  paripes 

Tanyp. 

Total 

G  p.  biomass 

Total  biomass 

Mar.  00 

19,442  ± 

361  ± 

20,085  ± 

15,721.2  ± 

15,863.1  ± 

2,262 

79 

2,275 

1,660.2 

1,662.0 

Apr.  00 

10,295  ± 

181  ± 

10,508  ± 

8,942.3  ± 

9,182.7  ± 

2,187 

43 

2,228 

1,479.6 

1,465.9 

May  00 

6,143  ± 

149  ± 

6,317  ± 

3,302.4  ± 

3,392.3  ± 

1,303 

30 

1,302 

462.3 

461.0 

Jun.  00 

8,565  ± 

72  ± 

8,670  ± 

7,142.3  ± 

7,1 84.9  ± 

1,982 

27 

2,009 

1,594.9 

1,385.5 

Jul.  00 

7,749  ± 

76  ± 

7,898  ± 

5,639.5  ± 

5,672.6  ± 

991 

21 

1,018 

793.8 

797.7 

Aug.  00 

5,206  ± 

100  ± 

5,306  ± 

2,673.7  ± 

2,729.6  ± 

1,342 

46 

1,355 

619.6 

622.6 

Sep.  00 

8,747  ± 

34  ± 

8,897  ± 

6,01 1.4  ± 

6,040.2  ± 

2,282 

26 

2,292 

1,387.6 

1,386.3 

Oct.  00 

6,808  ± 

70  ± 

6,884  ± 

3,799.9  ± 

3,830.9  ± 

1,022 

28 

1,026 

565.0 

567.2 

Nov.  00 

14,507  ± 

61  ± 

14,568  ± 

5,096.8  ± 

5,132.9  ± 

3,652 

13 

3,646 

1,130.9 

1,127.9 

Dec.  00 

7,588  ± 

56  ± 

8,091  ± 

3,077.5  ± 

3,182.0  ± 

2,382 

19 

2,460 

844.1 

857.9 

Jan.  01 

10,816± 

26  ± 

11,149  ± 

6,530.0  ± 

6,601.4  ± 

2,005 

13 

2,030 

1,178.0 

1,182.5 

Feb.  01 

6,167  ± 

34  ± 

6,210  ± 

5,040.5  ± 

5,072.7  ± 

1,174 

16 

1,172 

860.9 

859.1 
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Table  9.  Mean  ±  SD  of  Glyptotendipes  paripes  larval  dry  biomass  (mg/larva)  by  1-mm 
size  class  and  cube  root  of  mean  value  used  for  calculating  a  size/biomass  curve. 

Size  class  (mm)  Mean  dry  biomass                 Cube  root  of  mean  dry  biomass 
(mg/larva)  ±  SD 

3  0.015  ±0.008  0.24 

4  0.032±  0.012  0.32 

5  0.055  ±0.021  0.38 

6  0.120  ±0.048  0.49 

7  0.241  ±0.096  0.62 

8  0.344  ±0.109  0.70 

9  0.479  ±  0.207  0.78 

10  0.833  ±0.409  0.94 

11  1.270  ±0.419  1.08 

12  1.658±0.330  1.18 

13  1.729  ±0.866  1.20 

14  2.586  ±0.449  1.37 

15  2.861  ±1.396  1.42 

16  3.133  ±0.503  1.46 

Pupa 2.147  ±0.905 £ 

Not  required  for  length/biomass  curve  but  used  for  productivity  estimates 
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Figure  15.  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes paripes 
in  Lake  Dora  (Lake  County,  Florida),  March,  May  1999. 


62 


50 


49    - 


48  - 


47   - 


<~?  46   - 

O 

Z 

o 

oo 

<N  45 


July  1999 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

$       45       44       43       42      41       40       39       38 
3    50 


4> 

H    49  - 

I 

48  - 


47  - 


46  - 


45 


September  1999 


nun  in  |ii  n  n  ii  i|i  mi  ii  ii|  n  ii  i ii 1 1  |ii  n  ii  ii  1 1 1  ii  mi  ii|  i  mi  n  n 

46  44  43  42  41  40  39  38 


Longitude  (Minutes  81°  West) 
0-0.5  1.0-1.5  2.0-2.5  3.0-3.5 


0.5-1.0 


1.5-2.0 


2.5-3.0 


3.5-4.0 


4.0-4.5 


4.5-5.0 


Figure  16.  Spatial  larval  distributions  log(No.  /  m2  +1)  of  Glyptotendipes  paripes 
in  Lake  Dora  (Lake  County,  Florida),  July,  September  1999. 
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Figure  17.  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes paripes 
in  Lake  Dora  (Lake  County,  Florida),  November  1999,  January  2000. 
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Figure  18.  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes paripes 
in  Lake  Dora  (Lake  County,  Florida),  March,  April  2000. 
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Figure  19.  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes paripes 
in  Lake  Dora  (Lake  County,  Florida),  May,  June  2000. 
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Figure  20.  Spatial  larval  distributions  log(No.  /  m2  +1)  of  Glyptotendipes paripes 
in  Lake  Dora  (Lake  County,  Florida),  July,  August  2000. 
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Figure  2 1 .  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Dora  (Lake  County,  Florida),  September,  October  2000. 
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Figure  22.  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Dora  (Lake  County,  Florida),  November,  December  2000. 
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Figure  23.  Spatial  larval  distributions  log(No.  /  m2  +  1)  of  Glyptotendipes paripes 
in  Lake  Dora  (Lake  County,  Florida),  January,  February  2001 . 
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Figure  24.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Yale  (Lake  County,  Florida),  March,  May  1999. 
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Figure  25.  Spatial  larval  distributions  log(No./m2  +  1))  of  Glyptotendipes paripes 
in  Lake  Yale  (Lake  County,  Florida),  July,  September  1999. 
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Figure  26.  Spatial  larval  distributions  log(No./m2  +1)  of  Glyptotendipes paripes 
in  Lake  Yale  (Lake  County,  Florida),  December  1999,  February  2000. 
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Figure  27.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Yale  (Lake  County,  Florida),  March,  April  2000. 
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Figure  28.  Spatial  larval  distributions  log(No./m2+  l)of  Glyptotendipes paripes 
in  Lake  Yale  (Lake  County,  Florida),  May,  June  2000. 
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Figure  29.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipes paripes 
in  Lake  Yale  (Lake  County,  Florida),  July,  August  2000. 


76 


o 
% 

o 
00 

rs 

CO 

o 


i 


bl  - 

October  2000 

56  - 

-^^^k. 

55  - 

•5^ 

54- 

^^i^fc 

53  - 

^m 

i          i          i          i 

47 


46 


45 


43 


0.0-0.5 
0.5-1.0 


Longitude  (Minutes  81°  West) 
1.0-1.5         |  2.0-2.5 
1.5-2.0         I  2.5-3.0 


42 


3.0-3.5 
3.5-4.0 


4.0-4.5 
4.5-5.0 


Figure  30.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Yale  (Lake  County,  Florida),  September,  October  2000. 
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Figure  31.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Yale  (Lake  County,  Florida),  November,  December  2000. 
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Figure  32.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipes  paripes 
in  Lake  Yale  (Lake  County,  Florida),  January,  February  2001 . 
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Figure  33.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipesparipes 
in  Lake  Wauburg  (Alachua  County,  Florida),  March,  April  2000. 
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Figure  34.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipesparipes 
in  Lake  Wauburg  (Alachua  County,  Florida),  May,  June  2000. 
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Figure  35.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipesparipes 
in  Lake  Wauburg  (Alachua  County,  Florida),  July,  August  2000. 
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Figure  36.  Spatial  larval  distributions  log(No./m2  +1)  of  Glyptotendipesparipes 

in  Lake  Wauburg  (Alachua  County,  Florida),  September,  October  2000. 
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Figure  37.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipesparipes 

in  Lake  Wauburg  (Alachua  County,  Florida),  November,  December  2000. 
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Figure  38.  Spatial  larval  distributions  log(No./m2  +  1)  of  Glyptotendipesparipes 
in  Lake  Wauburg  (Alachua  County,  Florida),  January,  February  2001 . 
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Figure  39.  Canonical  correspondence  analysis  (CCA)  of  selected  water  and  sediment 
physico-chemical  parameters  with  Chironomidae  larval  community  in  Lake 
Dora  (Lake  County,  FL),  March  2000  -  February  2001 
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Figure  40.  Canonical  correspondence  analysis  (CCA)  of  selected  water  and  sediment 
physico-chemical  parameters  with  Chironomidae  larval  community  in  Lake 
Yale  (Lake  County,  FL),  March  2000  -  February  2001 
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Figure  41.  Canonical  correspondence  analysis  (CCA)  of  selected  water  and  sediment 
physico-chemical  parameters  with  Chironomidae  larval  community  in  Lake 
Wauburg  (Alachua  County,  FL),  March  2000  -  February  2001 


CHAPTER  6 
FISH  PREDATION  ON  CHIRONOMID  MIDGES  IN  STUDY  LAKES 


Introduction 
A  significant  source  of  G.  paripes  larval  reduction  is  likely  through  predation  by 
fish  within  the  habitat.  Previous  work  on  food  habits  offish  in  central  Florida  habitats 
includes  the  work  of  Huish  (1950)  concerning  bluegill  (Lepomis  macrochirus),  black 
crappie  (Pomoxis  nigromaculatus),  and  redear  sunfish  {Lepomis  microlopus)  from  Lake 
George.  Consumption  of  all  life  stages  of  chironomids  by  bluegill  comprised  24%  of 
total  volume  of  food  items  ingested.  Overall  midge  consumption  by  black  crappie  was 
29%  of  total  volume,  while  Chironomidae  formed  only  5.1%  of  total  volume  of  food 
items  ingested  by  redear  sunfish.  However,  Wilbur  (1969)  found  midges  to  be  of  greater 
importance  in  the  diet  of  redear  sunfish  in  Lake  Griffin,  with  the  genera  Clinotanypus, 
Glyptotendipes,  Harnischia,  and  Tanytarsus  constituting  81%  of  food  volume  consumed 
by  these  fish  in  the  lake.  Conversely,  chironomids  comprised  only  4%  of  food  volume  of 
redear  sunfish  from  Lake  Weir.  Wilbur  (1969)  indicated  that  redear  sunfish  feed  upon  a 
variety  of  prey  items,  depending  upon  availability.  Chew  (1974)  studied  the  early  life 
history  of  the  Florida  largemouth  bass  {Micropterous  salmoides  floridanus),  including  a 
diet  analysis.  Chironomids  were  of  minor  importance  to  young  largemouth  bass  in  Lake 
Weir,  being  8.54%  of  total  food  volume  for  fish  1.3-3.0  cm  total  length,  3.23%  for  fish 
3.1-7.5  cm,  and  0.02%  for  fish  7.6-29.6  cm  and  not  present  in  the  diet  of  larger  fish. 
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Midges  were  more  important  to  small  largemouth  bass  in  Lake  Hollingsworth,  forming 
20.5%  of  total  food  volume  for  fish  1.3-3.0  cm  long,  but  utilization  of  midges  as  food 
source  reduced  with  larger  size  classes,  forming  0.15%  for  fish  3.1-7.5  cm  and  0.0%  for 
fish  7.6-29.6  cm. 

The  objective  of  this  study  was  to  determine  species-specific  fish  predation  on 
chironomid  midge  larvae.  For  this  purpose,  midge  predatory  fish  were  collected  from 
Lakes  Dora  and  Yale  on  four  occasions.  Chironomid  larvae  recovered  from  gut  contents 
were  identified  and  the  composition  of  midge  taxa  in  the  fish  diet  was  compared  with  the 
composition  of  standing  crop  of  midge  larvae  in  the  study  lakes  during  the  same  month  to 
elucidate  any  relationships  between  consumption  of  midge  larvae  by  predatory  fish  and 
the  larval  composition  and  distributions  in  the  lakes  at  time  of  feeding. 

Materials  and  Methods 
Fish  were  collected  on  four  occasions  (May  1999,  July,  September  and  December 
2000)  from  Lakes  Dora  and  Yale  by  electrofishing  using  a  Coffelt  Manufacturing 
(Flagstaff,  AZ)  model  Mark  22  system  operated  from  a  7.3-m  pontoon  boat.  Collections 
were  made  in  the  morning  between  0830  and  1200.  These  collections  were  made  with 
the  permission  of  the  Florida  Fish  and  Wildlife  Conservation  Commission  under 
Fisheries  Scientific  Collectors  Permit  #s  FC-99-14  and  FC-2000-14.  Up  to  twenty  fish 
were  obtained  from  each  lake  per  session  and  fish  were  collected  from  near-shore  areas, 
usually  with  obvious  submerged  or  emergent  aquatic  vegetation.  Collected  fish  were 
identified  and  killed  immediately  by  a  blow  to  the  head  and  stored  on  ice  until 
transported  to  the  laboratory,  where  they  were  stored  at  - 1 0  °C  until  examined.  After 
thawing,  maximum  total  length  of  each  fish  was  measured,  foregut  dissected  (Bowen 
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1996),  and  the  gut  contents  transferred  to  4-dram  vials  containing  70%  ethanol  for 
storage.  Gut  contents  were  examined  under  variable  magnification  of  a  dissecting 
microscope  and  all  contents  enumerated.  Chironomidae  head  capsules  and  associated 
fragments  were  wet  mounted  on  slides  and  examined  at  400x  magnification  using  a 
phase-contrast  microscope  and  identified  to  lowest  possible  taxonomic  level  using  the 
keys  of  Epler  (1995). 


Results  and  Discussion 
All  fish  collected  during  this  study  were  bluegill  (Lepomis  macrochirus).  Five  to 
20  fish  successfully  were  collected  per  sampling  occasion  (Table  10),  with  more  fish 
collected  from  Lake  Dora  than  from  Lake  Yale.  During  May  1999  to  September  2000, 
captured  fish  from  Lake  Dora  were  larger  compared  to  Lake  Yale:  1 1.7-16.1  cm  and  8.9- 
10.8  cm,  respectively.  Midge  larvae  of  the  tribe  Tanytarsini  (mostly  Cladotany tarsus 
spp.)  were  the  most  numerous  in  gut  contents  offish  from  both  lakes  (Table  10),  ranging 
between  55.9-62.8%  of  total  consumed  midge  larvae  in  Lake  Dora  and  4.8-48.1%  of  total 
in  Lake  Yale.  Geoldichironomus  spp.  larvae  were  the  next  most  common  in  the  gut 
contents  offish  from  Lake  Dora  (0.0-27.5%  of  total  larval  chironomids).  Other  midge 
larvae  consumed  by  fish  collected  from  Lake  Dora  include  C.  crass icaudatus,  G.  paripes, 
Cryptochironomus  spp.,  Pseudochironomus  spp.  and  Tanypodinae.  Mean  number  of 
midge  larvae  in  fish  gut  contents  ranged  from  4.7  (December  2000)  to  44.0  (May  1999). 
In  all  samples  except  for  those  in  December  2000,  midge  larvae  were  present  in  the  gut 
contents  of  every  fish.  In  December  2000,  40%  of  collected  fish  had  empty  guts.  This 
could  have  been  due  to  low  water  temperatures  reducing  feeding  activity;  reduced  growth 
and  feeding  activity  of  bluegill  during  winter  months  has  been  previously  reported  by 
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Gilinsky  (1984).  Pseudochironomus  spp.  larvae  were  the  second  most  prevalent  midge 
larvae  in  the  fish  gut  contents  of  Lake  Yale,  forming  up  to  46.2%  of  total  midge  larvae, 
these  were  followed  by  C.  crassicaudatus  (collected  only  during  May  1999),  G.  paripes, 
Cryptochironomus  spp.,  and  Tanypodinae.  Fewer  midge  larvae  were  consumed  by  fish 
in  Lake  Yale,  with  mean  number  between  1.0  (December  2000)  and  18.9  (July  2000) 
larvae  per  fish.  Bluegill  feeding  on  midge  larvae  in  Lake  Yale  was  also  reduced  during 
December  2000,  though  only  one  fish  had  an  empty  gut.  Other  food  items  identified 
from  fish  in  these  two  lakes  included  immatures  of  Insecta  (Odonata,  Ephemeroptera  and 
Trichoptera),  Crustacea  (Decapoda,  Amphipoda  and  Ostracoda),  Nematoda,  Oligochaeta, 
Gastropoda,  and  some  unidentifiable  material. 

On  three  of  four  sample  occasions  (May  1 999,  September  and  December  2000), 
bluegill  in  Lake  Dora  had  consumed  significantly  more  midge  larvae  than  in  Lake  Yale 
(Mest;  t  =  3.438,  P  <  0.01,  /  =  4.478,  P  <  0.01,  t  =  3.669,  P  <  0.01,  respectively).  No 
clear  relationship  was  detected  between  consumed  number  of  larvae  and  fish  size. 

To  determine  qualitative  and  quantitative  fish  predation  on  chironomid  larvae, 
percent  composition  of  chironomid  larvae  in  gut  contents  of  collected  fish  was  examined 
in  relation  to  overall  percent  composition  of  chironomid  larvae  in  study  lakes,  as  well  as 
percent  composition  of  midge  larvae  in  the  nearshore  areas  with  firm  sediment  (DW  > 
40%)  representative  of  the  areas  from  where  the  fish  were  caught  (Figures  42  and  43).  In 
Lake  Dora  (Figure  42),  Tanytarsini  were  the  most  common  midge,  exhibiting  similar 
percent  compositions  in  fish  gut  contents  in  the  entire  lake  as  well  as  in  the  nearshore 
areas  of  the  lake.  During  July,  September  and  December  2000,  Geoldichironomus  spp. 
larvae  in  gut  contents  represented  a  much  larger  percentage  of  total  chironomids  when 
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compared  to  their  percentage  composition  in  the  entire  lake  and  nearshore  populations. 
In  Lake  Dora,  no  significant  differences  by  Mest  were  noted  between  midge  larval 
percent  composition  in  gut  contents  and  the  prevailing  populations  in  the  nearshore  or 
entire  lake  area.  In  Lake  Yale,  there  was  less  consistency;  C.  crassicaudatus  was  the 
most  common  chironomid  in  gut  contents  during  May  1999,  but  was  only  a  minor 
component  (<20%)  of  overall  midge  community  in  the  lake  during  that  time.  During  July 
2000,  Pseudochironomus  spp.  and  Tanytarsini  comprised  a  similar  percentage  of 
nearshore  chironomids  and  in  the  gut  contents  of  collected  fish.  In  the  remaining  months 
of  the  study  in  Lake  Yale,  Tanytarsini  were  the  most  common  midge  larvae  in  fish  gut 
contents,  as  well  as  in  the  firm  sediments.  Glyptotendipes  paripes  was  the  most  common 
midge  in  Lake  Yale  during  July  and  September  2000  and  composed  a  large  component 
during  December  2000,  but  comprised  a  smaller  component  of  the  gut  contents.  This  was 
probably  because  G.  paripes  larvae  were  concentrated  in  deeper  areas  of  the  lake  and 
thus  were  not  in  the  immediate  grazing  area  of  the  collected  predatory  fish.  Percent 
composition  of  G.  paripes  larvae  in  gut  contents  and  overall  percent  composition  in  lake 
larval  population  was  significantly  different  (/  =  2.938,  P  -  0.026,  n  =  4),  but  there  was 
no  significant  difference  noted  with  percent  composition  of  firm  sediment  areas.  No 
significant  differences  were  noted  among  any  other  chironomid  taxa  collected  from  Lake 
Yale. 

Bluegill  (L.  macrochirus)  appear  to  be  indiscriminate  feeders  on  midge  larvae 
inhabiting  the  sandy  nearshore  areas  of  Lakes  Dora  and  Yale,  with  the  percent 
composition  of  midge  larvae  in  the  gut  contents  of  fish  being  similar  to  the  composition 
of  chironomid  communities  in  nearshore  areas  of  the  lakes.  In  Lake  Dora,  there  might 
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have  been  some  preference  indicated  by  the  data  for  consumption  of  Geoldichironomus 
spp.  larvae,  but  since  the  composition  differences  were  statistically  not  significant,  this 
result  could  be  a  sampling  artifact.  However,  both  Gilinsky  (1984)  and  Rieradevall  et  al. 
(1995)  reported  selectivity  offish  predation  on  chironomids;  primarily  larger  or  more 
accessible  chironomids  tended  to  be  preyed  upon  more  often,  while  Wolfram- Wais  et  al. 

(1999)  reported  selectivity  among  fish  species  for  sediment  dwelling  or  epiphytic 
chironomids.  The  significant  difference  between  percent  composition  of  G.  paripes  in 
benthos  of  Lake  Yale  and  in  the  gut  contents  offish  indicates  that  bluegill  collected 
during  this  study  were  not  feeding  upon  midge  larvae  in  the  deeper,  soft  sediment  areas 
of  the  lake,  possibly  concentrating  their  feeding  activity  in  the  shallower,  nearshore  areas. 
Butler  et  al.  (1992)  reported  that  midge  predatory  fish  were  preferentially  foraging  in 
restored,  sand  bottom  areas  of  Lake  Tohopekaliga.  Rieradevall  et  al.  (1995)  reported  that 
various  fish  species  in  Lake  Banyoles  (Spain)  demonstrated  spatial  preferences  for  food 
grazing  activity.  Lepomis  gibbosus  was  reported  to  preferentially  forage  for  food  in  an 
emergent  reed  belt  region  of  Neusiedler  See  (Austria)  (Wofram- Wais  et  al.  1999). 
Considering  the  presented  data,  bluegills  perhaps  do  not  have  a  major  influence  on 
overall  chironomid  midge  larval  numbers  in  these  lakes.  Further  work  would  be  needed 
to  enumerate  actual  standing  crop  of  bluegill  in  these  lakes  to  determine  the  total 
consumption  of  midge  larvae  in  relation  to  the  chironomid  standing  crop.  Rieradevall  et 
al.  (1995),  Butler  et  al.  (1992)  and  Gilinsky  (1984)  reported  fish  predation  significantly 
lowered  chironomid  larval  density  in  various  studied  habitats.  However,  Batzer  et  al. 

(2000)  reported  that  fish  predation  in  a  marsh  weedbed  did  not  significantly  reduce 
chironomid  densities  when  compared  to  enclosures  that  excluded  predatory  fish.  The 
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authors  explained  that  the  fish  predation  reduced  other  macroinvertebrates  predacious  on 
chironomids,  resulting  in  little  net  change  in  predation  on  midge  larvae.  Rasmussen 
(1990)  reported  that,  while  constituting  the  major  part  of  the  diet  of  whitefish  in  Hjarbaek 
Fjord,  Denmark,  only  2-5%  of  total  chironomid  productivity  was  consumed  by  the  fish. 
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Table  10.  Mean  length  and  gut  contents  of  bluegill  (Lepomis  macrochirus)  collected  on 
four  occasions  between  May  1999  and  December  2000  by  electrofishing  in  Lakes 
Dora  and  Yale  (Lake  County,  Florida).  Mean  length  and  number  of  chironomid 
larvae  per  fish  ±  SD  and  percent  composition  (shown  in  parentheses). 


Taxa 

May 

July 

September 

December 

1999 

2000 

2000 

2000 

Lake  Dora 

No.  fish  collected 

10 

20 

20 

20 

Fish  length  (cm) 

11.7  ±0.9 

15.1  ±3.4 

16.1  ±  1.6 

11,2±3,1 

Chironomus 

0±0 

0.1  ±0.3 

0±0 

0.8  ±2.6 

crassicaudatus 

(0.0) 

(0.3) 

(0.0) 

(16.1) 

Glyptotendipes 

0.4  ±0.7 

1.8  ±2.5 

0.2  ±0.7 

0.3  ±1.1 

paripes 

(0.1) 

(6.0) 

(0.5) 

(5.4) 

Cryptochironomus 

0±0 

2.5  ±3.4 

4.4  ±7.3 

0.3  ±  1.1 

spp. 

(0.0) 

(8.2) 

(10.4) 

(6.5) 

Pseudochironomus 

0±0 

0±0 

0.5  ±  1.5 

0±0 

spp. 

(0.0) 

(0-0) 

(1.2) 

(0.0) 

Geoldichironomus 

0±0 

8.3  ±11.0 

10.5  ±13.1 

0.7  ±1.7 

spp. 

(0.0) 

(27.5) 

(24.7) 

(14.0) 

Tanytarsini 

25.1  ±18.2 

16.9±17.5 

26.7  ±22.1 

2.6  ±6.3 

(57.0) 

(56.3) 

(62.8) 

(55.9) 

Tanypodinae 

0.7  ±  1.3 

0.5  ±  1.1 

0.1  ±0.3 

0±0 

(1.6) 

(1.5) 

(0.2) 

(0.0) 

Other 

17.8  ±7.8 

0±0 

0.1  ±0.4 

0.1  ±0.4 

chironomids/No  ID 

(40.5) 

(0.0) 

(0.2) 

(2.1) 

Total 

44.0  ±21.0 

30.0  ±31.6 

42.5  ±34.2 

4.7  ±  12.0 

Lake  Yale 

No.  fish  collected 

5 

7 

15 

10 

Fish  length  (cm) 

10.8  ±3.5 

9.2  ±2.0 

8.9  ±  1.6 

11.8±1.3 

Chironomus 

6.2  ±13.9 

0±0 

0±0 

0±0 

crassicaudatus 

(75.6) 

(0.0) 

(0.0) 

(0.0) 

Glyptotendipes 

0.2  ±0.4 

1.7*4.1 

0±0 

0.2  ±0.4 

paripes 

(2.4) 

(9.0) 

(0.0) 

(20.0) 

Cryptochironomus 

0±0 

1.7  ±  3.1 

0.5  ±1.3 

0.2  ±  0.4 

spp. 

(0.0) 

(9.0) 

(18.5) 

(20.0) 

Pseudochironomus 

0±0 

8.6  ±9.6 

0.6  ±1.6 

0.2  ±0.4 

spp. 

(0.0) 

(46.2) 

(22.2) 

(20.0) 

Tanytarsini 

0.4  ±0.9 

6.0  ±6.4 

1.3  ±  1.3 

0.3  ±  0.5 

(4.8) 

(32.2) 

(48.1) 

(30.0) 

Tanypodinae 

1.4  ±  3.1 

0.4  ±  1.1 

0±0 

0.1  ±0.3 

(17.1) 

(2.1) 

(0.0) 

(10.0) 

Other 

0±0 

0.5  ±  1.1 

0.3±  0.8 

0±0 

chironomids/No  ID 

(0.0) 

(2.6) 

(11.1) 

(0.0) 

Total 

■_             ■ 

8.2  ±  13.5 

18.9  ±  17.2 

2.7  ±3.0 

1.0  ±0.9 

Larvae  could  not  be  identified  due  to  damage  or  loss  of  head  capsules 
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Figure  42.  Composite  percent  composition  of  midge  larvae  in  bluegill  (Lepomis 

macrochirus)  gut  contents  and  field  populations  in  Lake  Dora,  May  1999- 
December  2000.  C.  crassicaudatus  (C.  c),  G.  paripes  (G.  p.), 
Cryptochironomus  (Cryp),  Pseudochironomus  (Pseu),  Geoldichironomus 
(Geol),  Tanytarsini  (Tant),  Tanypodinae  (Tanp),  and  other  (Oth). 
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Figure  43.  Composite  percent  composition  of  midge  larvae  in  bluegill  (Lepomis 

macrochirus)  gut  contents  and  field  populations  in  Lake  Yale,  May  1999- 
December  2000.  C.  crassicaudatus  (C.  a),  G.  paripes  (G.  p.), 
Cryptochironomus  (Cryp),  Pseudochironomus  (Pseu),  Geoldichironomus 
(Geol),  Tanytarsini  (Tant),  Tanypodinae  (Tanp),  and  other  (Oth). 


CHAPTER  7 
DETERMINATION  OF  TEMPERATURE-RELATED  DEVELOPMENT  OF 

Glyptotendipes  paripes 


Introduction 
To  understand  the  development  dynamics  of  an  insect  species,  development  time 
in  relation  to  temperature  must  be  investigated  for  that  species.  An  important  and  pivotal 
concept  of  this  type  of  work  is  degree-days  (DD),  representing  a  heat  over  time  unit, 
usually  expressed  as  elapsed  days  at  degrees  Celsius  over  a  lower  threshold  value  to 
complete  given  development  stages.  This  concept  has  been  derived  from  work  with 
agricultural  plants,  investigating  temperature  requirements  for  crop  growth  to  harvest 
(Arnold  1959,  1960).  Pruess  (1983)  reviewed  uses  of  degree-day  methods  for 
entomology  and  pest  management.  Techniques  for  estimating  DD  requirements  include 
direct  laboratory  measurements  and  regression  (Campbell  et  al.  1974,  Luckmann  et  al. 
1976),  use  of  sine- wave  estimates  from  maximum-minimum  temperature  field 
measurements  (Allen  1976),  and  fluctuating  field  temperatures  (Messenger  1964,  Eubank 
et  al.  1973).  For  this  study,  modifications  of  the  methods  used  by  Stevens  (1998), 
Campbell  et  al.  (1974)  and  Luckmann  et  al.  (1976)  were  employed  to  estimate  the  DD 
requirements  of  G.  paripes  collected  in  central  Florida.  The  objective  of  this  study  was 
to  provide  information  needed  for  simulating  temperature  effects  on  development  of 
larval  G.  paripes  populations  in  the  computer  model. 
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Materials  and  Methods 

Adult  G.  paripes  were  collected  by  an  aspirator  during  early  evening  on  several 
occasions  in  March  2001  from  Ferran  Park,  Eustis,  central  Florida  on  the  eastern  shore  of 
Lake  Eustis,  a  eutrophic  lake  in  the  Ocklawaha  chain  of  lakes.  Collected  adults  were 
released  into  a  30  x  30  x  30  cm  screen  cage  provided  with  a  10-cm  diameter  dish 
containing  150-250  mL  preaerated  IX  Martin's  rearing  solution  (Martin  et  al.  1980)  for 
oviposition.  On  the  following  morning,  observed  egg  masses  were  collected  by  pipette  to 
provide  material  for  rearing  experiments. 

For  the  rearing  experiment,  sections  of  an  egg  mass  containing  30-60  eggs  were 
placed  into  tissue  culture  plate  wells  (Costar  model  3513,  12  well  cluster,  Corning  Inc., 
Corning,  NY)  with  4-mL  preaerated  IX  Martin's  solution,  one  mass  section  per  plate  in 
well  #1  using  a  modification  of  Stevens  (1998).  Plates  were  incubated  at  a  temperature 
range  of  10.0,  12.5,  15.0,  17.5,  20.0,  22.5,  25.0,  27.5,  30.0,  32.5,  and  35.0  °C,  with  a 
14:10  hour  light:dark  photoperiod  maintained  in  separate  growth  chambers  (Model  E- 
30B,  Percival  Mfg.  Co.,  Boone,  IA).  Three  plates  (replicates)  per  temperature  in  a 
growth  chamber  were  used,  with  each  replicate  from  a  separate  egg  mass.  Plates  were 
examined  at  12-hour  intervals  and  media  replaced  every  24  hours  with  fresh  preaerated 
IX  Martin's  solution.  Once  eggs  hatched  and  larvae  had  exited  the  egg  mass  sheath,  one 
larva  was  transferred  by  pipette  to  each  of  the  12  wells  (containing  4-mL  preaerated  IX 
Martin's  solution)  in  the  same  plate  (replicate),  and  excess  larvae  discarded.  After 
removal  of  excess  larvae  and  sheath  remains,  fresh  preaerated  media  was  placed  in  well 
#1.  The  12-hour  examination  interval  was  maintained  and  any  developmental  changes 
were  recorded  (a  molt  was  considered  to  occur  after  an  associated  exuvia  was  observed). 
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Upon  emergence,  sex  of  the  imago  was  recorded.  Each  developmental  event  was 
considered  to  occur  halfway  between  examinations.  Overall  lengths  of  larvae  after  each 
instar  molt  were  estimated  with  a  dissecting  microscope  equipped  with  an  eyepiece- 
mounted  micrometer.  Media  replacement  was  continued  daily  with  4  mL  fresh 
preaerated  IX  Martin's  solution  and  35.0  uL  food  suspension.  The  food  suspension  was 
prepared  by  mixing  1 00  mg  K.9    fish  food  (Friskies,  Victoria,  Australia,  provided  by 
M.  M.  Stevens)  and  25  mg  baker's  yeast,  then  macerating  with  a  tissue  grinder  in  12  mL 
IX  Martin's  solution.  Stevens  (1993,  1998)  recommended  the  use  of  one  whole  grain  of 
K9®  fish  food  per  well  (replaced  daily),  but  preliminary  work  with  G  paripes  indicated 
that  the  food  needed  to  be  finer  grained  and  that  additional  nutrients  were  needed.  After 
molting  to  4th  instar,  a  small  amount  of  sand,  sterilized  (held  at  100  °C  for  24  hours)  and 
sieved  (250-um  mesh),  was  added  to  each  well  to  provide  substrate  for  tube  building. 
Stevens  et  al.  (2001)  reported  that  G.  paripes  do  not  survive  as  well  in  culture  with  the 
tissue  paper  strips  recommended  for  Chironomus  tepperi  and  that  sand  was  more 
suitable.  Upon  emergence,  sex  was  determined  and  the  imago  discarded.  Due  to  time 
and  mechanical  limitations,  the  experiment  was  terminated  before  full  emergence  of  all 
test  larvae  at  all  temperatures.  Data  were  organized  and  analyzed  using  Excel  97® 
(Microsoft  Inc.,  Redmond,  WA)  and  Instat  Version  3.00  (Graphpad  Software,  San  Diego, 
CA). 

Results  and  Discussion 
Some  difficulties  were  experienced  with  the  mechanical  reliability  and  consistent 
operation  of  the  growth  chambers  and  in  adjusting  the  larval  diet  for  G.  paripes.  As  a 
result,  there  was  an  inconsistent  mortality  and/or  loss  of  replicates  during  the  course  of 
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the  study.  Where  possible,  replacement  sets  of  replicates  were  run.  Since  the  larvae  were 
commonly  active  under  light,  precise  measurements  of  overall  larval  length  were 
difficult,  so  no  statistics  were  calculated  for  larval  length,  but  some  general  observations 
were  made.  At  hatch,  larvae  were  approximately  0.75  mm  in  length,  after  molt  to  2nd 
instar,  about  1 .5  mm,  then  3 -mm  at  3  rd  instar,  and  4th  instar  were  between  6  and  12-mm. 
Some  larvae  were  observed  to  reach  15 -mm  in  length,  but  these  became  sluggish  and 
appeared  to  be  less  flexible,  and  died  before  pupating.  Observed  mean  development 
times  are  summarized  in  Table  1 1 .  Eggs  held  at  10  °C  failed  to  hatch  after  30  days  and 
were  discarded.  Larvae  held  at  35  °C  hatched  normally  but  survived  less  than  one  week 
and  failed  to  molt  to  2n  instar.  At  most  temperatures,  mean  time  to  adult  emergence  was 
slightly  lower  for  males,  though  statistically  not  significant  from  females  by  Student's  t- 
test.  The  sex  ratio  of  emerged  adults  was  close  to  even,  53.1%  female  and  46.9%  male. 

The  linear  regression  technique  described  by  Campbell  et  al.  (1974)  and 
Luckmann  et  al.  (1976)  were  used  to  estimate  the  lower  temperature  development 
threshold  (To).  This  method  plots  percent  development  to  next  stage  per  day  as  a 
function  of  temperature,  with  the  x-intercept  of  the  equation  equal  to  To.  This  was 
modified  using  the  mean  development  time  instead  of  the  median  time  because  the  mean 
values  gave  more  consistent  linear  regression  results  (Table  12).  The  overall  To  value  of 
9.0  °C  was  calculated  from  the  mean  of  To  values  derived  from  significant  (P  <  0.05) 
regressions  (hatch,  2n  ,  3rd  instar  and  adult)  (Table  12).  The  lack  of  significance  in  the 
regressions  for  development  to  4th  instar  and  to  pupa  is  likely  a  result  of  the  wide 
variability  in  development  times  noted  for  these  stages  (Table  1 1).  The  reason  for  this 
high  variability  and  inconsistency  is  not  clear,  though  qualitative  or  quantitative 
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deficiencies  in  the  provided  diet  are  strongly  suspected  in  this,  such  as  those  discussed  by 
Rasmussen  (1984)  or  Sankerperumal  and  Pandian  (1991). 

Once  T0  was  determined,  the  DD  requirements  for  G.  paripes  development  (Table 
13)  were  estimated  by  calculating  the  DD  above  T0  for  each  life  stage  using  the 
respective  mean  values  from  Table  1 1 .  The  total  DD  estimate  of  716  is  greater  than  the 
preliminary  estimate  of  600  DD  >  0  °C  proposed  for  this  species  by  Ali  et  al.  (1985) 
based  on  a  weak  periodicity  of  daily  adult  emergence  patterns  from  Lake  Monroe,  central 
Florida.  While  some  of  this  difference  may  be  a  result  of  variability  in  the  Ali  et  al. 
(1985)  data  (as  indicated  by  the  weak  autocorrelation),  the  remaining  difference  is  likely 
due  to  the  previously  mentioned  dietary  deficiencies.  Rasmussen  (1984)  reported  that 
food  deficits  were  the  likely  reason  for  slow  early-spring  growth  of  both  Chironomus 
riparius  and  G.  paripes  in  a  prairie  pond  in  Alberta,  Canada.  Both  species  were  reported 
to  be  univoltine  in  that  habitat.  The  development  requirement  for  C.  riparius  was 
reported  to  be  475  DD  >  4.1  °C  in  the  laboratory  (requirements  were  not  reported  by 
these  authors  for  G.  paripes).  However,  under  field  conditions  they  reported  that  over 
820  DD  were  required  for  C.  riparius  to  reach  early  4th  instar.  Food  availability  had  a 
considerable  effect  on  the  development  of  Chironomus  circumdatus  at  constant 
temperatures  (Sankarperumal  and  Pandian  1991).  At  37  °C,  larval  development  time 
increased  from  8  days  to  16  days  as  the  concentration  of  live  Chlorella  sp.  cells  provided 
as  food  was  reduced  from  25  x  106  cells/mL  to  13  x  106  cells/mL. 

Other  studies  have  reported  a  fairly  wide  range  of  DD  requirements  for 
chironomid  development.  Mackey  (1977)  reported  estimated  To  and  development  time  at 
15  °C  values  from  laboratory  studies  of  12  lotic  midge  species  from  the  Thames  River, 
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UK.  In  this  study,  the  author  reported  T0  values  between  2.4  and  9.7  °C  and  days  to 
develop  at  15  °C  from  5.4  to  35.8.  Stevens  (1998)  reported  150.5  DD  >  10.5  °C  for  male 
C.  tepperi  and  167.1  DD  >  10.3  °C  for  females  of  the  same  species.  While  not  directly 
determining  DD,  Kimerle  and  Anderson  (1971)  reported  a  considerable  variability  in 
development  times  for  Glyptotendipes  barbipes  reared  at  24  °C,  with  up  to  a  20  day 
range  for  molting  to  4th  instar  in  a  cohort.  Considerable  variations  within  a  given 
temperature  were  also  noted  for  Chironomus  decorus,  Endochironomus  nigricans  and 
Chironomus  staegeri  emergence  (Danks  1978),  with  a  difference  of  up  to  34  days 
between  first  and  last  emergence.  Chironomus  decorus  were  reported  to  have  an  unusual 
growth  rate  in  relation  to  temperature,  forming  a  distinct  bell  curve  with  maximum 
growth  at  20  °C,  a  lower  limit  of  8-10  °C,  and  an  upper  tolerance  of  28-30  °C  (Maier  et 
al.  1990).  Glyptotendipes  paripes  larvae  were  reared  from  egg  to  adult  in  19  days  at 
27  °C  on  a  natural  substrate  supplemented  with  bakers  yeast  (Cowell  and  Vodopich 
1981). 

This  work  provides  a  first  detailed  look  at  the  temperature  requirements  for 
G.  paripes  development;  however,  there  is  still  further  work  needed  to  clarify  these 
requirements  and  to  test  if  these  findings  need  to  be  amended.  In  this  study,  the  simplest 
linear  regression  method  was  used.  Commonly,  the  growth  curve  tends  to  be  logarithmic, 
with  a  slow  response  at  cooler  temperature,  a  linear  section  through  most  of  the 
temperature  range,  then  a  slowing  growth  rate,  and  finally  a  downturn  toward  the  upper 
thermal  limit  of  the  organism  (Campbell  et  al.  1974).  This  study  indicates  the  linear 
portion  of  the  growth  curve  is  between  10  and  32  °C,  with  a  lower  threshold  of  9.0  °C 
and  an  upper  threshold  of  35  °C.  Development  at  32.5  °C  for  4l  instar  to  pupa  was 
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slower,  but  it  is  unclear  from  these  data  as  to  whether  this  is  a  temperature  or  a  diet  effect. 
Future  studies  will  examine  dietary  effects  on  development  and  attempt  to  clarify  the 
nonlinear  portions  of  the  growth  curve. 


Table  1 1 .  Mean  ±  SD  development  time  (hours)  for  each  life  stage  of  Glyptotendipes 
paripes  reared  at  a  range  of  constant  temperatures  ±  1  °C  and  14:10  hour 
lightdark  photoperiod.  Number  of  individuals  (n)  included  in  calculations  are 
shown  in  parentheses. 


Temperature 
°C 

ls,Instar 

2nd  Instar 

3rd  Instar 

4th  Instar 

Pupa 

Adult 

10.0 

naa 

12.5 

(36) 
198  ±30 

632  ±  19 

nac 

nac 

nac 

nac 

15.0 

(72) 
138  ±  0 

(17) 
443  ±  68 

439  ±56 

nac 

nac 

nac 

17.5 

(24) 
70  ±6 

(23) 
234  ±  48 

(22) 
234  ±  42 

343  ±68 

nac 

nac 

20.0 

(36) 
66  ±0 

(35) 
211  ±  35 

(34) 
215  ±  33 

(30) 

247  ±  76 

384  ±0 

72  ±0 

22.5 

(36) 
58  ±6 

(26) 
197  ±25 

(26) 
205  ±  26 

(24) 
271  ±94 

(1) 

474  ±71 

(1) 
48  ±0 

25.0 

(36) 

42  ±0 

(17) 
200  ±21 

(17) 
203  ±18 

(15) 
364 ± 126 

(6) 

355  ±73 

(5) 

42  ±7 

27.5 

(36) 
38  ±6 

(15) 
127  ±20 

(15) 
136  ±  21 

(12) 

233  ±43 

(5) 
363  ±109 

(4) 
39±11 

30.0 

(36) 
34  ±6 

(H) 

123  ±39 

(ID 

126  ±25 

(10) 
229  ±  65 

(8) 

327 ± 102 

(8) 

36±11 

32.5 

(36) 
30  ±0 

(21) 
105  ±14 

(21) 
111  ±14 

(19) 
153  ±43 

(13) 

424  ±  73 

(13) 

25  ±4 

(36) 

(36) 

(36) 

(36) 

(19) 

(18) 

35.0 

30.0  ±0 
(48) 

nab 

''No  egg  hatch  was  observed  up  to  30  days. 

bComplete  mortality  occurred  before  molting  to  2n  instar. 

Experiment  terminated  before  test  organisms  developed  to  stage. 
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Table  12.  Slope  and  x-intercept  from  linear  regression  equations  of  mean  development 
times  per  life  stage  of  Glyptotendipes  paripes  with  temperature  used  to  estimate 
lower  temperature  development  threshold  (T0).  Number  of  data  points  used  for 
each  stage  is  shown  in  parentheses. 


Stage 

Slope 

x-intercept 
(To) 

r2 

P 

Hatch 

3.396 

8.964 

0.985 

<  0.0001 

(9) 
2nd  Instar 

0.910 

8.288 

0.947 

<  0.0001 

(9) 
3rd  Instar 

0.843 

7.609 

0.927 

<  0.0001 

h(8) 
4th  Instar 

0.413 

1.298 

0.558 

0.053 

(7) 
Pupae 

0.040 

-130.8 

0.047 

0.680 

(6) 
Adult 

3.926 

10.945 

0.923 

0.0023 

(6) 

Table  13.  Estimated  mean  ±  SD  degree-days  above  9.0  °C  for  each  development  stage  of 
Glyptotendipes  paripes. 

1st  Instar         2nd  Instar  3rd  Instar  4th  Instar  Pupa  Adult 


29.8  ±2.8      104.5  ±14.0      108.5  ±15.0      166.3  ±45.6      277.6  ±79.1      29.3  ±2.8 


CHAPTER  8 

DEVELOPMENT  OF  A  SPATIO-TEMPORAL  COMPUTER  MODEL  FOR 

Glyptotendipes  paripes  LARVAL  DISTRIBUTIONS  IN  CENTRAL  FLORIDA  LAKES 


Introduction 
The  objective  of  this  study  was  to  develop  a  qualitative  spatio-temporal  model  of 
G.  paripes  larval  populations  in  central  Florida  lakes  using  data  from  the  current  study 
for  calibration  and  to  qualitatively  validate  the  model  using  previous  data  from  Lakes 
Jesup  and  Monroe,  central  Florida.  Limited  larval,  water  depth  and  sediment  DW  data 
from  previous  work  on  G.  paripes  in  Lakes  Jesup  and  Monroe  (Ali  et  al.  1998a,  Ali  and 
Lobinske,  unpublished  data)  provided  the  initial  concept  for  the  model.  Data  gathered  in 
the  present  study  were  used  to  develop  and  calibrate  an  exploratory  computer  model  of 
G  paripes  spatial  and  temporal  larval  distributions.  This  model  will  enhance  the  ability 
of  lake  managers  to  plan  sampling  and  management  programs  of  G.  paripes  by  providing 
preliminary  information  on  area(s)  of  likely  larval  distributions  of  G  paripes  in  a  target 
habitat,  allowing  for  more  efficient  use  of  resources.  This  model  will  also  provide  a  tool 
for  theoretical  experimentation  on  G  paripes  larval  populations.  Various  habitat 
scenarios  can  be  tested  by  manipulating  selected  environmental  variables  within  the 
model;  hypothetical  distributions  could  be  determined,  which  may  provide  insight  to 
possible  population  dynamics  and  testable  scenarios  for  field  populations. 

Since  insects  develop  through  multiple  immature  life  stages  before  adulthood,  a 
model  examining  such  a  population  requires  that  age-structure  be  included.  This  allows 


106 


107 

each  life  stage  to  be  individually  considered  during  the  run  of  the  model,  with  separate 
survival  and  development  calculations  for  each  stage.  To  examine  a  heterogeneous 
ecosystem  such  as  a  lake,  the  model  will  also  need  to  include  spatially  specific 
information  on  habitat  suitability  previously  indicated  in  this  study,  such  as  water  depth 
and  sediment  physical  composition  in  the  case  of  G.  paripes  larvae.  Allen  et  al.  (1996), 
Brewster  and  Allen  (1997)  and  Brewster  et  al.  (1997)  developed  specific  techniques  to 
incorporate  specific  habitat  requirements  into  computer  models  of  insect  populations 
using  the  software  package  Matlab®  (The  Math  works  Inc.,  Natik,  MA).  These  models 
use  matrix  functions  to  provide  discrete  spatial  information  on  resources  and  modeled 
populations.  The  current  model  expands  on  these  techniques  to  include  more  detailed 
habitat  information  in  the  model  that  can  be  adjusted  by  initial  conditions. 

Materials  and  Methods 
Spatial  and  temporal  distributions  of  G.  paripes  were  described  by  examining 
each  life  stage  (egg,  larval  instars,  pupa  and  adult)  in  individual  spatial  grids  and  time 
steps.  Each  simulated  lake  was  expressed  as  a  50  x  50  spatial  grid  matrix.  The 
population  is  defined  by  a  50  x  50  x  7  three-dimensional  matrix  (Figure  44).  At  each 
geographic  grid  unit  (x,  y),  the  z  value  represents  the  population  of  a  given  life  stage 
(egg,  4  larval  instars,  pupa  and  adult).  For  example,  the  three-dimensional  coordinates 
(10,  10,  1)  would  represent  the  number  of  eggs  at  geographic  grid  unit  (10,  10),  while 
(5,  10,  6)  would  represent  the  pupal  population  at  geographic  grid  unit  (5,  10).  A 
Lefkovitch  population  matrix  (Lefkovitch  1965)  was  used  to  simulate  survival  and 
development  at  each  time-step  in  individual  spatial  grid  units  (Figure  45).  The  simulated 
habitat  was  defined  by  three  50  x  50  matrices  that  specified  lake  bathymetry,  sediment 
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DW,  and  other  habitat  modifiers  (such  as  vegetation  or  unusual  substrates)  that  were 
combined  to  produce  a  single  habitat  suitability  map  that  influenced  survival  rate.  Other 
modifiers  to  the  population  included  temperature  effects  on  rate  of  development,  food 
availability  influence  on  fecundity,  influence  of  density  on  survival  (density  dependence), 
and  dispersal  range  of  the  adult.  A  basic  flow  chart  of  the  computer  program  is  shown  in 
Figure  46  and  the  text  of  the  Matlab®  program  (m-file)  is  provided  in  the  Appendix. 

Within  the  Matlab®  program  shell,  the  program  (m-file)  is  started  by  entering  the 
file  name  into  the  command  line,  followed  by  arguments  defining  Secchi  disk 
transparency  (cm),  water  temperature  (°C)  and  deviation  of  water  level  from  historic 
mean  (m).  Vectors  (a  one-dimensional  array  of  numbers)  defining  maximum  survival 
rate  (Smax),  maximum  fecundity  (Fmax)  and  DD  requirements  (Chapter  7)  are  set  within 
the  m-file.  Each  of  these  vectors  contains  seven  numbers,  representing  the  value  for  each 
life  stage.  The  Smax  vector  used  maximum  values  of  0.9  for  egg  and  larval  instar  stages, 
and  0.6  for  pupae  and  adult.  Since  precise  values  were  not  known,  given  the  simple 
assumption  that  barring  other  factors,  most  larvae  would  survive,  0.9  would  be  a 
reasonable  estimate  that  includes  some  background  mortality  not  otherwise  included  in 
the  model.  Since  pupae  and  emerging  adults  are  more  exposed  to  predation,  a  reduced 
value  was  chosen.  Since  only  one  clutch  of  eggs  is  laid  by  an  adult  female  (Xue  et  al. 
1994),  all  but  the  last  value  in  the  Fmax  vector  were  zeros;  the  rate  of  the  reproductive 
stage  was  set  at  1 ,520  eggs/female  (760  per  individual  assuming  a  1 : 1  sex  ratio), 
determined  for  G.  paripes  by  Xue  et.  al  (1994).  The  functional  development  vector  (D) 
used  by  the  program  is  calculated  from  the  DD  vector  using  the  temperature  value  set  in 
the  command  line  and  calibrated  to  a  three-day  time  step. 
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The  relationship  of  G.  paripes  larval  density  with  sediment  DW  and  water  depth 
formed  the  core  of  the  habitat  suitability  maps  used  in  the  model.  Base  habitat  classes 
were  established  using  1-m  water  depth  intervals  and  10%  DW  intervals,  for  example, 
>l-2  m  depth,  >20-30%  sediment  DW.  The  habitat  class  with  the  highest  G.  paripes 
mean  density  (composite  of  all  three  study  lakes)  was  assigned  a  value  of  1 ,  and  all  other 
classes  were  assigned  a  proportional  0-1  value  based  on  mean  larval  density  in  that  class 
and  defined  habitat  suitability.  Maps  of  sediment  DW  and  water  depth  (at  mean  lake 
elevation)  were  prepared  for  each  lake  in  the  form  of  a  50  x  50  matrix.  These  matrix 
maps  were  prepared  by  overlaying  a  50  x  50  grid  onto  the  bathymetric  and  DW  maps 
derived  in  Chapter  4  and  recording  the  value  in  each  grid  unit.  Where  values  within  a 
grid  unit  were  variable,  the  median  value  was  used  in  the  matrix.  These  data  were 
organized  using  an  Excel  97®  (Microsoft  Inc.,  Redmond,  WA)  spreadsheet  and  saved  as 
tab-delimited  ASCII  text  files.  Separate  files  specific  to  each  lake  were  prepared.  These 
ASCII-text  files  were  imported  by  the  program  as  a  matrix.  During  simulations,  the 
water  depth  matrix  map  was  adjusted  by  the  user-defined  water  elevation  value. 
Temporary  maps  were  generated  using  true:false  (1:0)  statements  to  generate  maps 
indicating  each  DW  class  (by  10%  DW  increments)  and  water  depth  class  (by  1-m 
intervals).  These  matrix  maps  had  a  "1"  value  in  grid  units  that  matched  the  DW  or  depth 
class,  and  a  "0"  value  where  the  grid  units  did  not  match.  These  1 :0  water  depth  and  DW 
class  matrix  maps  were  combined  by  sequentially  adding  DW  data  class  by  water  depth 
class  to  assemble  a  new  matrix  within  the  model,  forming  resource  suitability  maps  that 
define  the  simulated  habitat.  Element-by-element  multiplication  by  a  correction  factor 
map  (also  prepared  as  a  tab-delimited  ASCII-text  file)  modified  the  above  map  for  the 
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effects  of  vegetation  and  unusual  substrates  found  in  individual  lakes  and  produced  the 
final  habitat  map. 

A  randomly  generated  50  x  50  x  7  matrix  represents  the  initial  G.  paripes 
population.  The  Matlab®  random  number  generator  was  reset  to  a  new  random  seed  each 
time  the  m-file  was  operated.  To  simulate  reproduction,  a  fecundity  matrix  was  produced 
of  the  same  dimensions  as  the  population  matrix.  This  matrix  was  first  set  up  as  a  50  x 
50  two-dimensional  matrix  with  a  value  of  1  in  each  cell.  An  element-by-element 
multiplication  with  a  modification  of  the  habitat  map  defined  the  confines  of  the  lake  in 
the  fecundity  matrix  (mapf).  Secchi  disk  transparency  (represented  a  rough  estimate  of 
chlorophyll  concentration)  was  used  to  represent  food  availability  and  was  defined  by  the 
user  for  each  simulation.  Since  Xue  et  al.  (1994)  indicated  a  direct  relationship  between 
adult  female  size  and  egg  production,  it  was  reasonable  to  consider  that  the  nutritional 
status  of  the  larvae  is  important  in  determining  adult  size,  hence  overall  fecundity. 
Optimum  phytoplankton  abundance  assumed  the  maximum  fecundity  rate,  while  other 
values  reduced  fecundity.  Based  on  current  and  previous  data  (Ali  and  Lobinske, 
unpublished  data),  Secchi  disk  transparency  of  30  cm  was  assigned  the  optimum  value. 
Secchi  values  above  30  cm  assumed  less  food  available,  whereas  values  under  30  cm 
assumed  adverse  effects  by  excessive  eutrophication.  The  mapf  matrix  was  modified 
using  simple  regression  equations  (Appendix)  based  on  this  trend.  This  matrix  was 
duplicated,  and  the  duplicate  versions  added  onto  the  first  matrix  as  new  cells  in  the  z 
direction  [using  the  "cat"  (concatenate)  command  within  Matlab®],  to  produce  a  three- 
dimensional  array  with  the  same  dimensions  as  the  population  matrix.  Element-by- 
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element  multiplication  of  the  z  values  in  each  geographic  grid  by  the  Fmax  vector 
produced  the  final  fecundity  matrix  (MAPF). 

Since  development  time  of  individual  life  stages  of  G.  paripes  is  not  equal,  a 
Lefkovitch  population  growth  matrix  (Lefkovitch  1965)  (Figure  45)  was  used  to  simulate 
survival  and  development  for  populations  at  each  spatial  grid.  To  accomplish  this,  two 
50x50x7  three-dimensional  arrays  were  prepared  by  concatenating  habitat  maps 
together  in  the  same  manner  as  the  fecundity  matrix  to  produce  the  base  survival  matrix 
(MAPS).  The  survive/develop  matrix  (representing  the  SnDn  entries  of  the  Lefkovitch 
matrix  in  Figure  45)  was  defined  by  multiplying  the  z  values  in  each  cell  of  MAPS  with 
the  Smax  and  D  vectors  to  produce  the  S/D  matrix.  A  non-develop  vector  was  defined  as 
1-D  and  was  used  in  place  of  the  D  vector  in  producing  the  survive/not  develop  matrix 
[S-ND,  representing  the  Sn(l-Dn)  entries  of  the  Lefkovitch  matrix  in  Figure  45]  using  the 
same  method  as  the  S/D  matrix. 

Dispersal  of  adults  was  simulated  by  a  filter  function  that  diffuses  the  simulated 
adult  population  in  each  spatial  cell  into  nearby  cells  based  on  a  normal  distribution 
centered  on  the  origin  cell  using  the  method  described  by  Allen  et  al.  (1996)  and 
Brewster  and  Allen  (1997).  Specifics  of  adult  dispersal  were  not  included  in  the  model; 
the  normal  distribution  was  chosen  to  provide  a  simple  dispersal  function  until  further 
data  are  available. 

The  time  loop  was  started  and  run  for  the  specified  number  of  time-steps  (tmax). 
Both  the  S/D  and  S/ND  matrices  were  multiplied  element  by  element  with  the  population 
matrix  and  this  result  modified  by  a  density-dependent  Ricker  survival  function: 
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Nt+1  =  Nte<-bN<> 
Where  N  is  population,  t  is  time  step  and  b  is  relative  strength  of  density  dependence 
(Ricker  1954).  Life  stages  that  the  Ricker  function  applies  to  (all  larval  instars  and  pupa) 
were  defined  within  the  program.  Carrying  capacity  was  not  determined  during  these 
investigations,  but  from  preliminary  trial  and  error  analysis  of  the  model,  a  value  of  0.001 
was  found  to  be  effective  for  strength  of  density  dependence  on  survival  in  the  equation. 
These  two  resulting  matrices  were  added  together  element-by-element,  the  dispersal 
function  was  applied  to  the  adult  stage,  MAPF  was  element-by-element  multiplied  to  the 
result  and  the  resultant  "eggs"  were  added  to  the  z  =  1  cells  of  the  population  matrix. 

Population  of  all  life  stages  at  each  grid  unit  was  totaled,  and  then  a  log  plot  of  the 
population  distribution  (No./grid  unit,  calibrated  to  represent  No./m2)  was  displayed. 
Mean  density  (No./grid  unit  =  No./m2)  was  determined  for  each  time  step  and  tracked. 
The  time  loop  was  repeated  with  the  final  population  at  the  end  of  the  previous  time  step 
used  as  the  start  population  in  the  new  time  step.  When  the  time  series  was  complete,  the 
final  population  (log  transformed)  density  distribution,  mean  population  trend  (No./m2) 
for  the  simulation,  and  the  habitat  map  were  plotted. 

For  validation  testing,  the  computer  model  was  run  the  for  ten  time  steps 
(representing  30  days)  using  water  level,  mean  Secchi  disk  transparency  and  mean  water 
temperature  recorded  concurrently  with  historical  larval  distributions  from  Lakes  Jesup 
and  Monroe  used.  To  test  with  the  calibration  data  sets,  the  model  was  run  with  the  same 
kind  of  data  recorded  with  larval  G.  paripes  collected  from  Lakes  Dora,  Yale,  and 
Wauburg.  The  observed  G.  paripes  larval  density  distributions  from  field  data  were 
qualitatively  compared  with  simulated  distributions  calculated  by  the  model  by  visual 
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comparison  of  location  and  magnitude  of  high  larval  density  areas.  Simulations  were  run 
for  all  monthly  data  sets  collected  during  this  study,  and  on  eight  data  sets  for  each  of  the 
validation  lakes  (Jesup  and  Monroe). 

To  examine  the  utility  of  the  model  for  generating  theoretical  populations  based 
on  different  environmental  conditions  (such  as  water  level)  or  biological  functions  (such 
as  maximum  fecundity),  the  model  was  run  for  50  time  steps  to  establish  a  baseline  data 
set  using  the  Lake  Monroe  water  depth,  sediment  DW,  and  habitat  modification  maps. 
One  variable  at  a  time  was  sequentially  varied  to  examine  the  influences  that  variable  had 
on  the  simulated  population.  Baseline  values  used  for  these  simulations  were  those 
described  above  for  development  rate,  fecundity,  and  strength  of  density  dependence. 
Dispersal  range  was  set  with  a  standard  deviation  of  2  grid  units  and  calculated  to  include 
an  area  3  SD  in  radius.  The  arbitrary  user  input  base  values  chosen  for  establishing  the 
baseline  output  were:  Secchi  disk  transparency  (60  cm),  water  temperature  (22  °C)  and 
water  level  (deviation  0  m  from  average  value).  Variables  tested  by  this  simulation  series 
were:  Secchi  disk  transparency,  water  temperature,  water  level,  fecundity,  strength  of 
density  dependence  on  survival,  and  dispersal  range. 

Results  and  Discussion 
Representative  validation  results  using  data  for  Lake  Jesup  during  May,  1996  are 
shown  in  Figure  47.  The  model  identified  relatively  high-density  areas  to  the  southeast 
of  the  small  island  in  the  center  of  the  lake,  along  sections  of  the  northern  shore  and  the 
northeastern  tip  near  the  confluence  with  the  St.  Johns  River.  The  model  population, 
however,  had  a  wider  larval  geographic  distribution  than  the  field  data.  There  was  a 
close  match  using  April  1996  validation  data  from  Lake  Monroe  (Figure  48),  with  model 
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and  field  data  in  agreement  with  each  other.  The  model  demonstrates  the  moderate  (100- 
1 ,000  larvae/m2)  population  densities  near  the  northwest  shore  and  in  the  eastern  bays, 
and  the  lower  densities  along  the  northern  and  southern  shorelines.  Similar  validation 
results  were  obtained  for  other  data  sets  from  both  Lake  Jesup  and  Lake  Monroe  (data  not 
shown).  This  type  of  model  has  previously  been  successfully  used  to  model  the  whitefly, 
Bemisia  argentifoli,  populations  in  a  crop-system  model  (Brewster  et  al.  1997). 
However,  these  authors  had  reported  that  some  portions  of  the  model  did  not  work  as 
well  as  others  and  that  further  research  was  required. 

It  is  also  important  to  see  how  the  model  describes  G.  paripes  population  patterns 
in  the  lakes  used  to  calibrate  the  model,  as  shown  in  Figures  49-52.  There  were  good 
qualitative  agreements  in  geographic  distributions  and  densities  between  the  model  and 
May  1 999  data  from  Lake  Dora  (Figure  49).  The  high  densities  along  the  southern  shore 
of  the  east  basin  are  not  expressed  in  the  model,  but  the  concentration  at  the  narrows 
between  the  east  and  middle  basins,  as  well  as  the  western  end  of  the  lake,  are  present. 
The  unexplained  diminished  and  widely  scattered  larval  field  populations  found  in  the 
lake  during  the  latter  part  of  this  study  were  not  simulated  well,  and  the  model  generally 
overestimated  density  and  area  of  distribution. 

The  model  produced  uneven  results  with  Lake  Yale  data  sets  (Figures  50  and  51). 
The  population  center  in  the  northeast  bay  is  shown  in  the  model  and  the  field  data  from 
May  1999,  as  well  as  the  modest  densities  along  the  shoreline  of  the  southern  bay,  though 
the  concentrations  in  the  field  data  for  the  northwest  areas  are  not  expressed  as  well  by 
the  model.  The  model  fit  was  poor  for  data  from  the  latter  part  of  the  study,  such  as  July 
2000,  as  shown  in  Figure  51 .  The  unusual  high  densities  of  G.  paripes  larvae  found  in 
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the  deep  portion  of  the  northern  basin  are  not  predicted  in  the  model,  while  the  peripheral 
concentrations  in  the  model  population  were  not  found  in  the  field  data.  The  factors 
resulting  in  the  shift  of  G.  paripes  larvae  to  the  deep-water  part  of  Lake  Yale  are  still  not 
understood  and  require  further  research  on  their  identification  and  quantification  for 
incorporation  into  the  model. 

The  high  densities  of  G.  paripes  larvae  found  in  the  southwestern  area  of  Lake 
Wauburg  during  July  2000  are  well  represented  in  the  simulation  model  (Figure  52), 
though  the  densities  along  the  eastern  shore  estimated  by  the  model  were  relatively 
reduced  compared  to  the  field  data.  A  similar  pattern  was  noted  with  repeated 
simulations  using  other  data  sets  from  Lake  Wauburg.  The  low  eastern  shore  estimate 
may  be  an  artifact  of  the  narrow  band  of  suitable  habitat  along  that  portion  of  the  lake  and 
the  50  x  50  matrix  used  for  the  simulation.  This  population  may  perhaps  be  better 
simulated  using  a  100  x  100  matrix  that  would  reduce  the  proportion  of  the  simulated 
population  that  dispersed  out  of  the  suitable  habitat  zone. 

The  baseline  output  of  the  program  (using  Lake  Monroe  data)  for  running 
theoretical  manipulations  is  shown  in  Figure  53.  In  all  simulations,  population  trend  is 
shown  as  mean  number  of  larvae/m2  and  density  distribution  as  log(No.  larvae/m2).  The 
population  trend  cycled  every  7-8  time  steps,  with  increasing  amplitude  and  peak  values 
between  400  and  1,700.  These  cycles  were  present  in  all  of  the  theoretical  simulations 
and  are  consistent  with  the  rough  cyclic  trend  in  adult  collection  data  reported  by  AH  et 
al.  (1983,  1985).  The  population  declines  sharply  after  each  peak,  but  the  lowest  value  of 
each  cycle  also  increased  over  time  (from  1  to  375).  The  density  distribution  showed 
highest  densities  along  the  northwestern  shore  and  the  southeastern  bays. 
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Varying  food  availability  (represented  by  Secchi  disk  transparency)  resulted  in 
considerable  differences  in  the  response  of  the  simulated  population  (Figure  54). 
Increased  food  (Secchi  disk  transparency  30  cm  compared  to  baseline  of  60  cm)  resulted 
in  peak  mean  populations  over  5,500,  and  the  cycle  lows  exceeding  1,500.  Variations 
within  each  cycle  increased  over  the  time  run,  creating  a  more  rapid  series  of  peaks 
during  the  last  ten  time  steps  of  the  run.  Larval  distribution  was  over  a  wider  area, 
including  portions  of  the  southern  shoreline  and  the  northeast  shoreline,  with  densities 
over  10,000  covering  a  much  larger  area  than  baseline  simulation.  With  decreased  food 
concentration  (Secchi  disk  transparency  90  cm),  population  peaks  were  considerably 
lower  than  the  baseline  values  (peak  values  between  200  and  550).  Distribution  was 
reduced  to  a  northwestern  band  (with  only  the  ends  reaching  1 ,000)  and  in  the 
southeastern  bays  (with  only  one  area  with  densities  over  1,000).  The  variability  of  the 
chlorophyll  data  in  relation  to  larval  populations  that  was  reported  in  Chapter  5  did  not 
lend  itself  to  conversion  into  the  model  process.  Since  Secchi  disk  transparency  can  also 
be  influenced  by  water  color  and  turbidity  from  other  sources  (such  as  zooplankton  or 
sediments),  further  work  is  needed  to  fully  include  chlorophyll  concentration  in  place  of 
the  less  accurate  Secchi  disk  transparency  value  for  estimating  food  availability  in  the 
model. 

Water  temperature  effects  on  the  model  produced  unusual  results  not 
characteristic  of  the  field  population  (Figure  55).  Increased  development  time  resulting 
from  lower  temperatures  (12  °C)  caused  high  mean  population  peaks  (up  to  7,000)  and 
cycle  minima  reaching  almost  2,000.  Density  distribution  was  similar  to  that  produced 
by  Secchi  disk  transparency  of  30  cm.  This  may  be  a  result  of  the  increased  residence 
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time  within  the  system  while  reproduction  is  still  occurring,  allowing  for  a  considerable 
accumulation  of  simulated  larvae.  In  central  Florida,  temperatures  this  low  are  not 
encountered  for  this  length  of  time.  Elevated  temperature  and  reduced  development 
times  produced  somewhat  larger  mean  population  peaks  than  baseline,  up  to  2,500,  and 
slightly  higher  peak  densities  in  the  northwest  shoreline  and  southeast  bay  areas,  plus 
extending  higher  densities  along  the  southern  shoreline. 

Decreasing  water  level  by  1  meter  from  mean  resulted  in  slightly  higher  mean 
population  peaks  than  baseline  (maximum  2,000)  that  seemed  to  stabilize  during  the  last 
two  cycles  recorded  (Figure  56).  Density  distribution  along  the  northwest  shore  shifted 
toward  the  south,  with  a  high  density  area  appearing  in  the  southwest  bay,  a  similar 
southward  shift  was  also  noted  in  the  eastern  section  of  the  lake,  and  more  of  the  southern 
shore  had  higher  densities.  Increasing  lake  level  by  1  meter  reduced  mean  population 
peaks  from  the  baseline  (maximum  1,300),  while  cycle  minima  increased  to  almost  600. 
Population  distribution  was  reduced  to  a  narrower  band  along  the  northwest  shore,  and 
primarily  limited  to  the  northernmost  portion  of  the  southeastern  bays.  These  results  are 
in  agreement  with  the  findings  of  Ali  et  al.  (1996),  who  reported  increased  populations  of 
G.  paripes  collected  by  New  Jersey  light  traps  along  the  southern  shore  of  Lake  Monroe 
when  the  lake  levels  were  reduced,  though  this  was  also  attributed  to  enhanced  nutrient 
concentrations  (from  the  existing  amount  of  nutrients  being  dissolved  in  a  smaller  volume 
of  water),  leading  to  increased  phytoplankton  populations. 

Altering  the  base  maximum  fecundity  of  the  simulated  population  produced  large 
changes  in  the  simulated  population  densities  (Figure  57).  Decreasing  this  value  to  760 
eggs/female  reduced  cycle  peaks  (after  the  initial  randomly  generated  peak)  to  1 50  or  less 


and  only  slightly  increased  amplitude.  Simulated  larval  distribution  was  limited  to  low 
levels  along  the  northwest  shore  and  the  southern  bay  on  the  east  side  of  the  lake,  with 
only  two  small  areas  approaching  a  density  of  1,000.  Doubling  fecundity  from  the 
baseline  value  resulted  in  much  higher  peaks  (up  to  5,500),  though  the  cycles  start  to  be 
erratic  and  more  frequent,  with  peak  values  slightly  decreasing  after  40  time  steps.  High- 
density  areas  cover  most  of  the  lake  perimeter  except  for  a  small  portion  of  the  south 
shore  and  two  small  areas  of  the  north  shore. 

The  greatest  differences  in  response  by  the  computer  model  were  generated  by 
changes  in  the  strength  of  density  dependence  on  survival  (Figure  58).  Increasing  the 
strength  of  density  influence  on  survival  from  0.001  to  0.01  reduced  population  peaks 
dramatically  from  baseline,  with  a  slight  decline  over  time,  and  an  overall  cyclic 
amplitude  decrease.  Overall,  low  to  modest  densities  of  larvae  were  present  in  the 
simulated  system  in  the  general  area  of  the  baseline  simulation,  with  no  areas  with 
densities  near  1,000.  Reducing  effect  of  density  on  survival  to  0.0001  resulted  in  the 
greatest  mean  population  peak  values,  reaching  13,000  during  the  simulation  run,  with 
broad  areas  of  the  simulated  distribution  with  densities  greater  than  10,000.  This 
dramatic  change  indicates  that  this  function  has  a  very  strong  influence  on  the  overall 
populations  in  lakes  and  the  factors  contributing  to  density-dependent  survival  should  be 
investigated. 

Decreasing  the  dispersal  range  (Figure  59)  of  the  simulated  population  to  a 
standard  deviation  of  1  grid  unit  resulted  in  a  more  aggregated  distribution  with  peak 
areas  approaching  densities  of  10,000.  The  population  trend  quickly  produced  higher 
overall  peaks  from  baseline,  with  more  erratic  cycles  and  peaks  developing  after  time 
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step  40.  Increasing  the  dispersal  area  to  a  standard  deviation  of  4  grid  units  spread  out 
the  distribution  considerably  to  cover  much  of  the  lake  with  lower  densities  and  lower 
mean  population  values.  These  trends  probably  resulted  from  loss  of  population  that 
dispersed  into  unfavorable  areas  of  the  lake  basin  where  survival  rates  were  lower. 

Based  on  the  preceding  information,  the  present  computer  model  was  moderately 
successful  in  predicting  G  paripes  larval  distributions  in  central  Florida  lakes  in  many 
circumstances.  These  successes  indicate  that  the  basic  structure  of  the  computer  model  is 
sound,  but  the  two  cases  of  the  model  deviating  from  field  data  (the  latter  portion  of  the 
field  data  from  Lakes  Dora  and  Yale)  show  that  further  work  is  required  to  elucidate  and 
include  additional  environmental  factors  that  influence  larval  G.  paripes  populations. 
Increased  quantification  of  the  physical  and  chemical  characteristics  of  lake  sediments  in 
relation  to  the  existing  larval  populations  may  also  be  required  to  understand  such 
deviations  between  model  and  field  data  as  the  shift  of  G.  paripes  populations  to  deeper 
water  and  soft  organic  sediments  during  summer  and  autumn  2000.  More  specific 
information  on  phytoplankton  community  that  the  larvae  feed  on  is  likely  needed,  since 
variability  in  food  quality  between  different  phytoplankton  species  may  play  a  role  in 
larval  survival  and  nutritional  status  (influencing  fecundity  rate).  Other  work  may 
include  determining  the  actual  carrying  capacity  of  different  sediments,  allowing  for  a 
more  precise  definition  of  density  dependence  on  survival.  This  model  used  simple 
assumptions  on  larval  G.  paripes  mortality;  detailed  development  of  a  life  table  would  be 
of  important  use  in  the  further  development  of  population  models  for  this  species.  Future 
work  will  also  include  operating  the  present  model  using  larger  matrix  maps  to  simulate 
the  lake  systems,  such  as  100  x  100  or  200  x  200  grid  units.  This  may  further  improve 
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the  operation  of  the  model  and  possibly  diminish  or  eliminate  any  scale-induced  effects, 
such  as  noted  for  the  eastern  shore  of  Lake  Wauburg. 

The  theoretical  simulation  runs  presented  in  this  study  show  the  utility  of  the 
model  for  doing  preliminary  examinations  of  larval  G  paripes  field  population 
distributions  under  various  conditions.  This  kind  of  theoretical  testing  could  be  useful  in 
designing  new  investigations  of  larval  populations,  and  examining  possible  manipulation 
effects  of  habitats  which  support  nuisance  levels  of  G.  paripes.  The  utility  of  this  type  of 
modeling  on  planning  experiments  has  already  been  demonstrated  by  Brewster  et.  al 
(1997)  in  a  study  of  organic  crop  arrangement  for  management  of  the  whitefly,  B. 
argentifoli. 
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Input  Secchi  disk  transparency,  water  temperature  and  water  level 


Set  maximum  fecundity  and  survival  vectors 
Calculate  development  vector 

V 

Import  sediment  dry  weight  map 
Import  water  depth  map  and  adjust  to  water  level 

V 

Sequential  setup  of  habitat  map  by  depth  and  sediment  dry  weight 

V 


Apply  habitat  modification  map 


V 


Set  up  dispersal  function 


V 


Set  up  fecundity  matrix 


V 


Set  up  survival/development  matrix  with  density  dependence 


V 


Begin  time  loop 

Survival  and  development 

Dispersal 

Reproduction 

Track  population 

End  loop 


V 


Plot  final  population  distribution 

Plot  population  trend 

Plot  resource  map 


Figure  46.  Flow  chart  of  Glyptotendipes  paripes  larval  distribution  model. 
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Figure  47.  Comparison  of  model  output  (bottom)  of  Glyptotendipes  paripes  larval 

distribution  in  Lake  Jesup  under  matching  conditions  with  May  1996  field 
collected  data  (top).  Scale  on  right  is  log(  No./m2). 
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Figure  48.  Comparison  of  model  output  (bottom)  of  Glyptotendipes paripes  larval 

distribution  in  Lake  Monroe  under  matching  conditions  with  April  1996  field 
collected  data  (top).  Scale  on  right  is  log(No./m2). 
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Figure  49.  Comparison  of  model  output  (bottom)  of  Glyptotendipes  paripes  larval 
distribution  in  Lake  Dora  under  matching  conditions  with  May  1999  field 
collected  data  (top).  Scale  on  right  is  log(No./m2). 
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10  20  30  40  50 

Figure  50.  Comparison  of  model  output  (bottom)  of  Glyptotendipes  paripes  larval 
distribution  in  Lake  Yale  under  matching  conditions  with  May  1999  field 
collected  data  (top).  Scale  on  right  is  log(No./m2). 
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Figure  51.  Comparison  of  model  output  (bottom)  of  Glyptotendipes  paripes  larval 
distribution  in  Lake  Yale  under  matching  conditions  with  July  2000  field 
collected  data  (top).  Scale  on  right  is  log(No./m2). 
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Figure  52.  Comparison  of  model  output  (bottom)  of  Glyptotendipes paripes  larval 
distribution  in  Lake  Wauburg  under  matching  conditions  with  July  2000 
field  collected  data  (top).  Scale  on  right  is  log(No./m2). 
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Figure  53.  Baseline  Glyptotendipes paripes  mean  population  trend  (No./m2) 

(top)  from  computer  model  run  for  50  time  steps  and  final  population 
log(No./m2)  distribution  (bottom). 
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Figure  54.  Glyptotendipes paripes  mean  No./m2  population  trends  from  computer 
model  run  for  50  time  steps  (left)  and  log(No./m2)  final  population 
distributions  (right).  Secchi  disk  transparency  30  (top)  and  90  cm  (bottom). 
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Figure  55.  Glyptotendipes paripes  mean  No./m2  population  trends  from  computer 
model  run  for  50  time  steps  (left)  and  log(No./m2)  final  population 
distributions  (right).  Water  temperature  12  °C  (top)  and  32  °C  (bottom). 
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Figure  56.  Glyptotendipes paripes  mean  No./m2  population  trends  from  computer 
model  run  for  50  time  steps  (left)  and  log(No./m2)  final  population 
distributions  (right).  Water  elevation  -1  m  from  mean  (top)  and  +1  m 
(bottom). 
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Figure  57.  Glyptotendipes paripes  mean  No./m2  population  trends  from  computer 
model  run  for  50  time  steps  (left)  and  log(No./m2)  final  population 
distributions  (right).  Fecundity  760  eggs  /  female  (top)  and  3,040 
eggs  /  female(bottom). 
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Figure  58.  Glyptotendipes paripes  mean  No./m2  population  trends  from  computer 
model  run  for  50  time  steps  (left)  and  log(No./m2)  final  population 
distributions  (right).  Strength  of  density  dependence  0.01  (top)  and 
0.0001  (bottom). 
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Figure  59.  Glyptotendipes paripes  mean  No./m2  population  trends  from  computer 
model  run  for  50  time  steps  (left)  and  log(No./m2)  final  population 
distributions  (right).  Dispersal  standard  deviation  1  grid  unit  (top)  and 
4  grid  units  (bottom). 


CHAPTER  9 
SUMMARY  AND  CONCLUSIONS 


The  current  study  has  expanded  the  existing  usable  database  on  chironomid  midge 
larval  populations  and  associated  water  and  sediment  physico-chemical  parameters  in 
central  Florida  lakes.  Larval  development  time  has  been  estimated  and  relative  predation 
pressure  by  bluegill  (Lepomis  macrochirus)  determined.  This  has  resulted  in  the 
development  of  a  qualitative  spatio-temporal  population  model  of  G.  paripes  larvae  in 
central  Florida  lakes. 

Detailed  sediment  and  bathymetric  maps  necessary  for  defining  larval  G.  paripes 
habitat  were  made  of  the  study  lakes.  Mean  water  depths  in  these  lakes  varied  from  2.7 
to  3.5  m.  All  of  the  study  lakes  had  slightly  acidic  sediments  (mean  6.35  to  6.77  pH 
units),  mean  total  organic  carbon  concentration  ranged  between  8.9  and  10.9%,  and  mean 
sediment  dry  weight  varied  from  30.6  to  41 .0%.  With  a  few  exceptions,  the  perimeter 
areas  of  these  lakes  had  sandy  substrates,  while  the  deeper  parts  had  soft,  organic 
sediments.  This  is  a  physical  effect  of  wave  action  preventing  accumulations  of  organic 
materials  in  the  nearshore  areas.  With  less  wave  action,  organic  material,  most  likely 
from  dead  and  decomposing  phytoplankton,  can  accumulate.  Since  decomposition 
bacteria  quickly  deplete  available  oxygen  in  these  sediments,  much  slower  anaerobic 
decomposition  occurs  in  the  sediments  deeper  than  only  a  few  millimeters.  Since 
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138 
deposition  from  the  dead  phytoplankton  in  these  lakes  is  faster  than  decomposition,  there 
is  a  net  increase  in  organic  sediments. 

Previously  reported  trophic  status  of  the  study  lakes  (Lake  Dora  as 
hypereutrophic,  Lake  Wauburg  as  eutrophic  and  Lake  Yale  as  mesotrophic  to  eutrophic) 
was  confirmed.  Glyptotendipes  paripes  abundance  in  these  lakes  was  influenced  by 
trophic  state,  with  increasing  abundance  from  mesotrophic  to  eutrophic,  then  declining 
abundance  with  hypereutrophic  conditions.  Species  of  Cladotany tarsus  were  numerically 
the  dominant  midge  larvae  in  Lakes  Dora  and  Yale,  while  G  paripes  was  dominant  in 
Lake  Wauburg.  Reasons  for  the  prevalence  of  Cladotanytarsus  spp.  in  Lakes  Dora  and 
Yale  were  not  clear  from  the  present  data.  A  possible  reason  is  the  generally  smaller  size 
that  may  allow  for  better  utilization  or  survival  in  nearshore  habitats. 

In  Lakes  Dora  and  Wauburg,  G.  paripes  were  distributed  primarily  on  firm 
sediments,  as  reported  previously  in  the  literature.  This  is  likely  a  result  of  more  suitable 
substrates  for  tube  building  and  less  oxygen  stress  when  compared  to  organic  substrates 
where  microbial  activity  rapidly  depletes  available  oxygen  in  the  interstitial  water  of 
sediments  and  particle  size  is  often  smaller.  However,  in  Lake  Wauburg,  relatively  high 
densities  extended  into  soft  organic  sediments  adjacent  to  the  normally  preferred  firm 
sediments  because  of  overcrowding.  While  G.  paripes  larval  distributions  in  Lake  Yale 
were  similar  to  Lakes  Dora  and  Yale  during  the  early  part  of  the  study,  an  unusual  shift  to 
deeper  water  and  soft,  organic  sediments  was  recorded  during  the  latter  part  of  the 
survey.  The  reason  for  this  shift  presently  remains  unexplained  and  requires  further 
research.  Initially,  moderately  high  populations  of  G.  paripes  were  recorded  in  Lake 
Dora,  declining  dramatically  during  the  first  year  of  the  study  and  remaining  low  during 
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the  remainder  of  the  investigation.  The  high  nutrient  loading  may  have  contributed  to 
this  lake  being  less  suitable  to  G.  paripes  populations.  The  cause  of  this  decline  remains 
unexplained.  Lake  Wauburg  supported  high  densities  (mean  up  to  19,442  larvae/m2)  of 
G.  paripes.  Favorable,  eutrophic  conditions  supplying  abundant  phytoplankton  as  food 
and  possibly  fewer  competing  chironomids  (as  shown  by  the  limited  midge  diversity  in 
this  lake)  contributed  to  these  high  densities.  Glyptotendipes  paripes  productivity  ranged 
from  rather  low  in  Lake  Dora  to  high  in  Lake  Wauburg,  where  the  annual  productivity 
estimate  of  1 56.9  g/m2  dry  weight  was  close  to  the  levels  reported  for  midge  larval 
habitats  such  as  sewage  treatment  ponds.  Canonical  correspondence  analysis  indicated 
that  sediment  parameters  including  dry  weight,  sand  presence,  muck  presence  and 
vegetation  were  the  primary  factors  influencing  larval  chironomid  distributions  in  the 
study  lakes.  These  detailed  relationships  between  larval  G.  paripes  populations  and 
environmental  conditions  allowed  for  the  inclusion  of  a  detailed  habitat  suitability  map 
within  the  model  that  determined  the  base  survival  rate  of  simulated  populations. 

Fish  gut  analysis  revealed  that  bluegill  were  nonselective  consumers  of  midge 
larvae  in  these  lakes  and  possibly  not  a  major  cause  of  G.  paripes  larval  reductions. 
Further  research  would  be  required  to  quantify  overall  bluegill  populations  in  the  lakes 
and  determine  the  overall  effect  of  predation  on  midge  larval  populations. 

Examination  of  the  temperature  influences  on  larval  development  revealed  a 
lower  development  threshold  of  9.0  °C  and  an  upper  threshold  between  32.5  and  35.0  °C, 
with  an  apparent  optimum  temperature  between  27.5  and  30.0  °C.  The  calculated 
development  time  in  the  laboratory  was  716  degree-days  above  9.0  °C  for  adult 
emergence  after  egg  laying.  While  this  value  is  higher  than  previous  indirect 
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measurements,  the  information  in  this  study  provides  the  first  detailed  examination  of 
temperature-dependent  development  of  this  species.  Additional  studies  should  include 
revision  and  possible  improvement  in  the  larval  diet  to  refine  the  degree-day  estimates. 
Determination  of  the  development  rate  of  larval  G.  paripes  enables  the  use  of  the 
Lefkovitch  population  matrix  in  determining  development  rate  of  the  simulated 
population  within  the  model.  Comparison  of  larval  length  measurements  in  this  study 
with  those  from  the  field  study  indicated  that  most  1st  and  2nd  instar  larvae  (3-mm  or  less 
in  length)  were  not  being  detected  by  the  field  collection  method  and  mostly  3rd  and  4th 
instar  larvae,  3-mm  or  greater  in  length,  were  retained.  Further  work  on  larval  G.  paripes 
populations  should  include  modifications  to  sampling  procedures  to  include  these  earlier 
instars. 

The  computer  model  generated  during  this  study  successfully  estimated  relative 
qualitative  G.  paripes  larval  densities  and  distributions  in  two  validation  lakes  (Lakes 
Jesup  and  Monroe)  and  in  the  three  study  lakes  under  most  conditions.  The  most  notable 
deviation  of  the  model  from  field  data  was  the  occurrence  of  G.  paripes  in  deep 
water/soft  sediment  observed  during  the  latter  part  of  the  study  in  Lake  Yale  which  was 
not  replicated  in  the  model,  as  well  as  the  overall  reduced  G  paripes  population  in  Lake 
Dora  not  effectively  estimated.  Due  to  these  deficiencies,  further  research  is  necessary  to 
improve  the  efficiency  of  the  model.  Improved  quantification  of  sediment  physical  and 
chemical  conditions  and  inclusion  of  more  parameters  for  better  understanding  the 
ecology  of  the  habitat  are  suggested.  Also,  qualitative  and  quantitative  compositions  of 
phytoplankton  community,  which  may  influence  midge  larval  diet,  need  to  be 
investigated. 
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While  work  is  still  needed  on  the  model,  it  could  be  recommended  for  use  by  lake 
managers  to  assist  in  sampling  and  treating  habitats  producing  nuisance  G.  paripes 
populations.  Once  proper  bathymetric  and  sediment  dry  weight  maps  are  generated  for  a 
given  lake,  quick  estimates  of  the  likely  areas  of  a  lake  supporting  large  or  nuisance 
levels  of  G.  paripes  populations  are  possible  which  should  enable  more  time-  and 
resource-efficient  sampling  of  the  populations  for  management  and  research  purposes. 
With  better  estimates  of  spatial  distributions,  more  time-  and  cost-efficient  control 
measures  can  be  used  to  reduce  midge  nuisance.  Selective  targeting  of  control  measures 
can  reduce  material  and  labor  costs,  as  well  as  minimize  nontarget  effects  within  a 
habitat.. 

This  research  provides  a  useful  new  tool  for  postulating  and  testing  new 
hypotheses  and  concepts  on  the  larval  ecology  of  G  paripes  and  possible  insights  for 
reducing  midge  nuisance.  The  model  can  allow  quick  examination  of  different 
environmental  conditions  or  situations  and  the  development  of  possible  hypotheses. 
These  could  then  be  tested  either  in  the  laboratory  or  in  the  field  and  could  be  aimed 
primarily  at  understanding  larval  ecology  or  developing  new  control  strategies. 
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APPENDIX 
MATLAB®  COMPUTER  PROGRAMS  USED  IN  STUDY 


Random  Generator  for  Sample  Locations 

function  dorasite(tmax) 

%  Format  dorasite(tmax) 

%  calculates  random  latitude  (minutes  at  28°N)  and  longitude  (81°W) 

%  tmax=number  of  sites  to  be  calculated 

rand('  state'  ,sum(  1 00*clock)); 
fort=l:tmax 

N=(rand*1.9)+46.4; 

W=(rand*6.1)+38.5); 
N_W=[N,W] 
end 


Simulation  Model  of  Glyptotendipes  paripes  Larval  Distribution 

function  lakegppop(secc,temp,el) 

%  Age  and  Space-Structured  Population  Projection  Model  Using  3-d  Arrays 

%  Density  Dependent  Survival  (Ricker  Function) 

%  The  3  dimensions  are  used  as: 

%      DIM1  &  DIM2  =  2-d  spatial  domain  (i  =  'North'  =  DIM1 ,  j  =  'East'  =  DIM2.) 

%      k  =  DIM3  =  Age  structure  at  location  (i,j). 

%    So  N(ij,k)  =  population  at  location  (i,j)  in  age  group  k. 

%    Parameters  are  similarly  structured  arrays. 

%    So  F(i,j,k)  =  fecundity  at  location(ij)  for  age  group  k. 

% 

%  FORMAT: 

%      lakegppop(secc,  temp,  el) 

%    where: 

%       secc  ■  Secchi  disk  transparency  (cm) 

%       temp  =  mean  water  temp  C 

%       el  =  lake  elevation  deviation  from  mean  (m) 

%       n  =  size  of  spatial  grid 

%        Fmax  =  Maximum  Fecundity  Vector  (reduced  by  habitat  map) 
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%  length(Fmax)  =  #  of  stages 

%  Smax  =  Maximum  Survival  Vector  (reduced  by  habitat  map) 

%  length(Smax)  =  #  of  stages 

%  D     =  Probability  of  development  vector, 

%  D(i)  =  P(developing  to  stage  i+1  in  a  time  step). 

%  length(D)  =  #  of  stages 

%  map  =  0-1  habitat  favorability  map 

%  m  =  age  group  that  disperses. 

%  s  =  std.  dev.  of  dispersal  pattern  (in  grid  cell  units)  (Try  2) 

%  ns  =  number  of  standard  deviations  to  make  kernel  (Try  3) 

%  tmax  °  run  time 

%         First  we  assign  (maximum)  fecundity  &  survival  vectors  by  age. 

Fmax  =  [000000  760];       %  Fecundity  vector  by  age 

Smax  =  [0.9  0.9  0.9  0.9  0.9  0.6  0.6];  %  base  survival  vector 

Dbase  =  [29.8  104.5  108.5  166.3  277.6  29.3  1];  %  Degree-Day  vector  for  development 

%       final  "1"  removes  all  holdover  adults  from  simulated  population 

Ddiv  =  [3333333];  %  conversion  factor  for  3  day  time-steps 

D2  =  Dbase/(temp-8.95);  %  degree-day  calculation,  8.95  =  lower  development  threshold 

D2  =  (Ddiv./D2);  %  adjusts  to  3-day  time  step 

D3  =  (D2<=1);  %  next  four  steps  adjust  D  so  that  no  value  >1 

D4  =  (D2>1); 

D5  =  D3.*D2; 

D  =  D5+D4;  %  final  development  vector 

m=7;  %  adult  dispersal 

s=2; 

ns=3; 

n  =  50;  %  size  of  spatial  domain 

tmax  =  50;       %  length  of  time  to  run 

area  =  346;  %  number  of  cells  defining  lake 

close  all  %  Close  any  open  figure  windows 
set(0,'Defaultaxesfontsize',14)  %  Set  ALL  figure  fontsizes  larger. 
rand('state',sum(100*clock));  %  Initialize  the  random  generators  with  a  random  seed. 

% Import  sediment  dry  weight  map,  set  up  range  matrices 

load  lakedw.txt  %  sediment  dry  weight  data  as  a  .txt  file 

dw=lakedw;  %  sets  0,1  matrices  for  each  dw  class 

wl=(dw==l); 

w2=(dw==2); 

w3=(dw=3); 

w4=(dw==4); 

w5=(dw==5); 

w6=(dw==6); 

w7=(dw=7); 


154 


w8=(dw==8); 

% Import  depth  map 

load  lakedep.txt  %  mean  water  depth  as  a  .txt  file 

% adjust  depth  map  for  lake  elevation — 

dep=lakedep; 
dep=dep+el; 
dc=(dep~=el); 
dep=dep.*dc; 

%  - set  up  larval  density  output  matrix- 

gpout=zeros(n); 


%-- set  up  depth  range  matrices 

%    sets  0,1  matrices  for  each  depth  range 

dl=(dep>0)&(dep<=1.00) 
d2=(dep>l)&(dep<=2.00) 
d3=(dep>2)&(dep<=3 .00) 
d4=(dep>3)&(dep<=4.00) 
d5=(dep>4)&(dep<=5 .00) 
d6=(dep>5)&(dep<=6.00) 
d7=(dep>6)&(dep<=7.00) 
d8=(dep>7.00); 

%  set  up  density  by  depth  and  dry  weight 

%  element  by  element  multiplication  of  depth  maps  with  dw  maps,  *  habitat  suitability 

% to  1  meter 

d=dl.*wl*0.043; 

gpout=gpout+d; 

d=dl.*w2*0.003; 

gpout=gpout+d; 

d=dl.*w3*0.01; 

gpout=gpout+d; 

d=dl.*w4*0.15; 

gpout=gpout+d; 

d=dl.*w5*0.876; 

gpout=gpout+d; 

d=dl.*w6*0.864; 

gpout=gpout+d; 

d=dl.*w7*0.600; 

gpout=gpout+d; 

d=dl.*w8*0.443; 

gpout=gpout+d; 
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% 1  to  2  meters 

d=d2.*wl*0.063; 

gpout=gpout+d; 
d=d2.*w2*0.034; 

gpout=gpout+d; 
d=d2.*w3*0.354; 

gpout=gpout+d; 
d=d2.*w4*0.428; 

gpout=gpout+d; 
d=d2.*w5*0.480; 

gpout=gpout+d; 
d=d2.*w6*0.813; 

gpout=gpout+d; 
d=d2.*w7*0.807; 

gpout=gpout+d; 
d=d2.*w8*  1.000; 

gpout=gpout+d; 

d=d3.*wl*0.192; 

gpout=gpout+d; 
d=d3.*w2*0.111; 

gpout=gpout+d; 
d=d3.*w3*0.181; 

gpout=gpout+d; 
d=d3.*w4*0.333; 

gpout=gpout+d; 
d=d3.*w5*0.296; 

gpout=gpout+d; 
d=d3.*w6*0.65; 

gpout=gpout+d; 
d=d3.*w7*0.533; 

gpout=gpout+d; 
d=d3.*w8*0.877; 

gpout=gpout+d; 

°/n ~X  fn  4  meters 

d=d4.*wl*0.310; 

gpout=gpout+d; 
d=d4.*w2*0.077; 

gpout=gpout+d; 
d=d4.*w3*0.091; 

gpout=gpout+d; 
d=d4.*w4*0.061; 

gpout=gpout+d; 
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d=d4.*w5*0.078; 

gpout=gpout+d; 

d=d4.*w6*0.227; 

gpout=gpout+d; 

d=d4.*w7*0.493; 

gpout=gpout+d; 

d=d4.*w8*0.922; 

gpout=gpout+d; 


% 4  to  5  meters- 

d=d5.*wl*0.321; 

gpout=gpout+d; 

d=d5.*w2*0.018; 

gpout=gpout+d; 

d=d5.*w3*0.002; 

gpout=gpout+d; 

d=d5.*w4*0.002; 

gpout=gpout+d; 

d=d5.*w5*0.046; 

gpout=gpout+d; 

d=d5.*w6*0.002; 

gpout=gpout+d; 

d=d5.*w7*0.007; 

gpout=gpout+d; 

d=d5.*w8*0.09; 

gpout=gpout+d; 


% 5  to  6  meters- 

d=d6.*wl*0.400; 

gpout=gpout+d; 

d=d6.*w2*0.002; 

gpout=gpout+d; 

d=d6.*w3*0.001; 

gpout=gpout+d; 

d=d6.*w4*0.001; 

gpout=gpout+d; 

d=d6.*w5*0.002; 

gpout=gpout+d; 

d=d6.*w6*0.001; 

gpout=gpout+d; 

d=d6.*w7*0.002; 

gpout=gpout+d; 

d=d6.*w8*0.05; 

gpout=gpout+d; 


% 6  to  7  meters- 
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d=d7.*wl*0.367; 

gpout=gpout+d; 

d=d7.*w2*0.001; 

gpout=gpout+d; 

d=d7.*w3*0.004; 

gpout=gpout+d; 

d=d7.*w4*0.005; 

gpout=gpout+d; 

d=d7.*w5*0.008; 

gpout=gpout+d; 

d=d7.*w6*0.002; 

gpout=gpout+d; 

d=d7.*w7*0.001; 

gpout=gpout+d; 

d=d7.*w8*0.001; 

gpout=gpout+d; 

% below  8  meters- 

d=d8*0.001; 
gpout=gpout+d; 


% final  habitat  modification- 
load  lakemod.txt; 
gpout=gpout.  *  lakemod; 
map=gpout; 


% Construct  Spatial  Filter  K  as  2-d  N(0,0,s)— 

% Calculate  the  dispersal  kernel  from  Normal  dist. 

Ksize  =  round(s*ns);    %  Size  of  K  is  2*round(s*ns)+l 
inc=0.01;    %  Spatial  step  for  integration 
if(Ksize==0) 

Ksize  =  1    %  Set  minimum  size  of  K  to  3x3  (-1  ->  1) 
end 

k=[]; 

for  i=-Ksize:Ksize 
Krow=[]; 
forj=-Ksize:Ksize 

ymin=i-l/2;  ymax=i+l/2; 

xmin=j-l/2;    xmax=j+l/2; 

vol=normalpiece(0,s,xmin,xmax,ymin,ymax,inc); 

Krow=[Krow  vol]; 
end 

K=[K;  Krow]; 
end 
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Kmiddle  =  K(:,(Ksize+l)); 

size(K); 

Kvol=sum(sum(K)); 

K  =  K/Kvol;  %  normalize  K  to  sum  to  one. 
close  all 

% random  seed  population  matrix 

N  =  rand(n,n,length(f));  %  Random  U(0,1)  initial  population. 

%  First  we  assign  (maximum)  fecundity  &  survival  vectors  by  age. 

f=Fmax;     %  fecundity  vector  by  age 
s  =  Smax;  %  survival  vector  by  age 
%  set  up  fecundity  and  survival  matrices 
mapf=ones(50); 


if(secc>=30) 

smod=22.8/(0.76*secc)  %  adjusts  for  phytoplankton  availability  measured  by  Secchi 
disk 
end 

if(secc<30) 

smod=(0.053*secc)-0.6  %  adjust  down  representing  excessive  eutrophication 
end 
mapf=smod*mapf; 

fmod=(lakedw~=0);  %  removes  any  population  outside  of  lake 
mapf=mapf .  *  fmod ; 
MAPF  =  mapf; 

fork=l:length(f)-l 

MAPF  =  cat(3,MAPF,  mapf);  %  Stack  length(f)  'maps'  on  top  of  each  other  end  along 
DIM3  dimension. 

%  Then  we  construct  the  3-d  arrays  for  Population,  Fecundity  &  Survival 


fork=l:length(f) 

F(:,:,k)  =  f(k).*MAPF(:,:,k);  %  Spatial  Fecundity  array, 
end 


figure 
colormap(jet) 
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SD  =  zeros(n,n,length(s));  %  Set-up  survival-develop  array  SD  with  initial  zeros. 

SND  =  zeros(n,n,length(s));  %  Set-up  survival-not-develop  array  SND  with  initial  zeros. 


b  =  0.001;  %  Strength  of  density  dependence 

p  =  2;     %  1  st  D.D.  Age  Group 

q  =  6;     %  Last  D.D.  Age  Group 

SD1  =  SD;    %  Initiate  SD1  same  size  as  SD 

SND1  -  SND; 

TotPop  =  sum(sum(sum(N)))/nA2; 

MAPS  =  map;%  S  vary  together  according  to  your  habitat  map. 

fork=l:length(s)-l 
MAPS  =  cat(3,MAPS,  map);  %  Stack  length(s)  'maps'  on  top  of  each  other 
end  %         along  DIM3  dimension. 


for  k=l:length(s) 
SD(:,:,k)  -  s(k).*D(k).*MAPS(:,:,k);  %  Spatial  Survival-Develop  array, 
end 

fork=l:length(s) 
SND(:,:,k)  -  s(k).*(l-D(k)).*MAPS(:,:,k);  %  Spatial  Survival-NOT-Develop  array, 
end 

for  t=l  :tmax    %  Begin  time  loop 


SDl(:,:,p:q)  =  SD(:,:,p:q).*exp(-b*N(:,:,p:q));  %  Ricker  D.D.  Survival-Stage  p:q 
SNDl(:,:,p:q)  =  SND(:,:,p:q).*exp(-b*N(:,:,p:q));  %  Ricker  D.D.  Survival-Stage  p:q 
N(:,:,m)  =  filter2(K,N(:,:,m));  %  Dispersal  by  kernel  K  @  stage  m. 
Nl  =  sum(F.*N,3)+SNDl(:,:,l).*N(:,:,l);  %  Nl  is  Reproduction+Surviving  Holdover 
N2  =  SDl(:,:,l:size(N,3)-l).*N(:,:,l:size(N,3)- 

l)+SNDl(:,:,2:size(N,3)).*N(:,:,l:size(N,3)-l); 

%  N2  is  survive-develop+survive-NOT-develop  holdover  at  each  stage. 
N  =  cat(3,  Nl,  N2);  %  Concatenate  reproduction  in  front  of  surviving  stages. 

TotN  =  sum(N,3);  %  Total  Population  at  each  location 

LogTotN  =  loglO(TotN+l);  %log  transform  of  population  data 

TotPop  =  [TotPop  sum(sum(TotN))/area]; 

imagesc(LogTotN)  %  View  log  transformed  total  population  map. 

caxis([0,5]); 

colorbar; 
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drawnow 
end  %  End  time  loop 

imagesc(LogTotN); 

caxis([0,5]); 

colorbar 

title(['Log  Total  Population  Map']); 

% Time  Plot  of  Total  Population  over  Space  &  Age 

figure 

plot(TotPop) 

Title(['Mean  population  density ']); 

%  end 

%  Look  at  resource  map. 

figure 

colormap(jet) 

imagesc(map) 

colorbar 

title('Habitat  Suitability  Map') 


Normal  Distribution  Generator  for  Model 

function  vol  =  normalpiece(mean,SD,xmin,xmax,ymin,ymax,inc); 

% 

%  FORMAT: 

%  vol  =  normalpiece(mean,SD,min,max,inc);    SEMICOLON  IMPORTANT!! 

%     WHERE: 

%         mean  =  mean 

%         SD  =  Standard  Deviation 

%         min  =  lower  bound  for  integration 

%         max  =  upper  bound  for  integration 

%         inc  =  integration  step  size  (try  0.01) 

%         vol  =  volume  under  the  normal  dist.  from  min  to  max  in  both  directions 

[X,Y]  =  meshgrid(xmin:inc:xmax,ymin:inc:ymax); 

var  =  SDA2; 

Z  =  (l/(2*pi*var))*exp(-((x-mean).A2+(y-mean).A2./(2*var)); 

vol  =  sum(sum(Z))*incA2; 

This  program  courtesy  of  Dr.  Jon  C.  Allen 
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